Next Article in Journal
The Protective Role of Melatonin in Sperm Cryopreservation of Farm Animals and Human: Lessons for Male Fish Cryopreservation
Previous Article in Journal
Animal Tourism: Thai Caregivers’ Perspectives on Their Relationships with Elephants and Tigers
Previous Article in Special Issue
Selection Response Due to Different Combination of Antagonistic Milk, Beef, and Morphological Traits in the Alpine Grey Cattle Breed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes

1
School of Life Sciences, Inner Mongolia University, Hohhot 010021, China
2
School of Veterinary Medicine, Inner Mongolia Agricultural University, Hohhot 010010, China
3
Institute of Animal Husbandry, Inner Mongolia Academy of Agricultural and Animal Husbandry Sciences, Hohhot 010031, China
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(6), 792; https://doi.org/10.3390/ani12060792
Submission received: 13 January 2022 / Revised: 17 March 2022 / Accepted: 17 March 2022 / Published: 21 March 2022

Abstract

:

Simple Summary

TBXT (c.333G > C; c.334G > T) has been identified as a molecular genetic marker in short-tailed sheep. This paper describes a high-resolution melting (HRM) analysis using unlabeled probes and asymmetric PCR for the detection of genetic variants of TBXT in short-tailed sheep populations. The detection results of this method are consistent with those of Sanger sequencing and can help farmers with marker-assisted breeding.

Abstract

The short-tailed phenotype has long been considered one of the best traits for population genetic improvement in sheep breeding. In short-tailed sheep, not only is tail fat eliminated but also the pubic area is exposed due to the lack of a tail covering, giving them an advantage in reproduction. Recent studies have shown that two linked mutations in sheep TBXT at nucleotides 333 and 334 are associated with the short-tailed phenotype. In the population of short-tailed sheep, several heterozygous mutants of this gene are found. In our research, we used high-resolution melting (HRM) to identify homozygous and heterozygous genotypes in a flock of short-tailed sheep and compared the results with those of Sanger sequencing, which were identical. This demonstrates that our established HRM method, a rapid and inexpensive genotyping method, can be used to identify homozygous and heterozygous individuals in short-tailed sheep flocks.

1. Introduction

The tails of most animals serve multiple functions in nature [1], including balance, movement, information exchange with peers [1,2], and swatting flies to remove them from the animals. After a long period of domestication, sheep with different tail length phenotypes have been developed [3,4], including the fat-long-tail, fat-short-tail, fat-waist-tail, slender tail, and thin-short-tail phenotypes. There is evidence [5,6] that the ancestors of sheep had a narrow tail phenotype that evolved into a fat-tail phenotype as they adapted to their frigid living environment. Sheep in China can be roughly divided according to their tail phenotypes into fat-waist-tail and thin-short-tail [3,7]. In several areas [8], long-tailed sheep are the main source of the fat intake of the local population but this phenotype presents many disadvantages [9], including that the overly long tail covers the pubic area, which makes mating inconvenient, and sheep with fat tails require more feed [10], which increases the cost of rearing. Although the tail docking of fat-tailed sheep is an option [11], this not only increases the labor of the farmer but is also exceedingly painful for the sheep; therefore, this is detrimental in terms of animal welfare [12,13] and, as a result, tail docking is controversial [11,13,14]. Short-tailed sheep accumulate more intramuscular fat than fat-tailed sheep, leading to more tender meat [14]. As a result, short-tailed sheep are more popular with buyers and producers are more likely to breed sheep with short tails [9,15,16].
Genes potentially associated with the formation of the short-tail phenotype in sheep have been found in many studies such as bone morphogenetic protein 2 (BMP2) [7,17], platelet-derived growth factor D (PDGFD) [16,17] and T-box transcription factor T (TBXT) [18,19], among which the mutation of TBXT is presumed to be related to the number of tail vertebrae. In a 1990 study [20], it was discovered that homozygous TBXT mutations could cause early fetal death in mice whereas mice with heterozygous TBXT mutations were short-tailed; this provided proof that this trait was linked to the TBXT gene. However, the embryonic development of brachyury, an encoded protein, also plays an extremely important role [21]. Brachyury, an important nuclear transcription factor throughout embryonic development [22], participates in and regulates numerous signaling pathways that influence critical processes such as mesoderm differentiation and somite segmentation [23,24]. Nucleotides 333 and 334 in sheep TBXT correspond with residues 111 and 112 in the T-box domain of the brachyury protein. Interestingly, most mice and cats [25,26] with the short-tailed phenotype are heterozygous for TBXT mutations and homozygous TBXT mutations have a negative influence on early embryonic development. In a 2013 study of Manx cats with a short-tailed phenotype [26], the researchers found that a heterozygous mutation in TBXT could cause a phenomenon known as haploinsufficiency. When one allele is mutated, the other allele can still be normally expressed but in cases of haploinsufficiency, having only half of the usual level of protein is not sufficient to maintain normal cell function; this is the phenomenon responsible for Manx cats having a short-tailed phenotype. The same hypothesis was proposed in a 2016 study [27] of cats with the short-tailed phenotype in Southeast Asia.
Hulunbuir sheep are a unique breed of Mongolian sheep in the Xini River Basin in the Hulunbuir region of the Inner Mongolia Autonomous Region of China [19,28], which can be divided according to the two main strains [29]: Barag sheep with the fat-long-tail phenotype; and short-tailed sheep with the short-tail phenotype. The very short tails of short-tailed sheep do not cover the anus whereas Barag sheep have tails that reach to their hocks or lower. With the exception of the tail phenotype, Barag sheep and short-tailed sheep are otherwise nearly identical [19], making these two sheep strains with different tail morphologies good biological models to research the mechanisms underlying the formation of short tails.
In the study of Zhi [19], the results of genome resequencing suggested that mutations in two linked nucleotides in TBXT, 333 and 334, might affect the number of caudal vertebrae in sheep. In the hybridization experiments of Han [18], several sheep with the short-tailed phenotype were found to have heterozygous mutations and individuals with heterozygous TBXT mutations could retain the short-tailed phenotype. This was different from the discoveries in mice and other mammals. TBXT homozygous mutations are not lethal during the embryonic development of short-tailed sheep. This suggests that the mechanism by which the TBXT mutation causes the short-tailed phenotype in sheep is different from that in mice and cats. Due to the widespread breeding of Barag sheep and short-tailed sheep in the Hulunbuir area, it is not uncommon for the offspring to be short-tailed with TBXT heterozygosity.
Farmers require a fast and low-cost genotyping method to ensure the breeding of homozygotes and to avoid progeny trait segregation during short-tailed sheep breeding, thereby reducing the time it takes to detect the genotype of progeny [30]. It can be difficult for farmers raising sheep with a short-tailed phenotype to distinguish between homozygotes and heterozygotes for the TBXT mutation. High-resolution melting (HRM) has been used to genotype various diseases and economic traits [31,32,33] because it is less expensive and time-consuming than Sanger sequencing and the TaqMan method [34]. The basic principle of HRM is that the heterozygous sample is amplified by PCR to generate DNA fragments and a mixture of two types of homoduplex DNA and two types of heteroduplex DNA will form upon the denaturation and renaturation of these fragments. A unique melting curve is formed during the melting of the resulting DNA fragments that can be distinguished from those of the homozygous DNA fragments [35]. Many mutated genes for economic traits in sheep and cattle have been identified by HRM analyses. These mutated genes include bone morphogenetic protein 15 (BMP15) [30], growth differentiation factor 9 (GDF9) [36], and major histocompatibility complex class II DRB3 (BoLA-DRB3) [37]. These well-established assay protocols demonstrate the benefits and viability of applying an HRM analysis in sheep breeding. In our research, we propose an HRM analysis using unlabeled probes and asymmetric PCR. The unlabeled probes ensure that the Tm difference between the homozygous and heterozygous samples is 0.3 °C, which allows the two samples to be clearly distinguished. The HRM method can be used to effectively select homozygous sheep, thus helping farmers to produce more sheep with the short-tailed phenotype.

2. Materials and Methods

2.1. Ethical Statement

All sheep experimental procedures and protocols were approved and authorized by the animal care and use committee of the Inner Mongolia Agricultural University (Inner Mongolia Autonomous Region, China) in this study (License No. 2020008).

2.2. Materials

Short-tailed sheep and Barag sheep were grown in Ewenki Autonomous County (Inner Mongolia Autonomous Region, China). We randomly chose 60 short-tailed sheep and 60 Barag sheep of the same age from 3 farms (40 on each farm). An equal number of rams and ewes was used to eliminate the effect of sex. A 2 mL sample of venous blood from each sheep was collected in a vacuum container containing EDTA and stored at −20 °C. Among the selected sheep flocks, 16 Barag sheep and 14 short-tailed sheep with known genotypes were selected to obtain caudal vertebrae X-ray images.

2.3. Methods

2.3.1. Genomic DNA Extraction

Genomic DNA was extracted using a FastPure Blood DNA Isolation Reagent kit (Vazyme, Nanjing, China) according to the manufacturer’s instructions. DNA concentration and purity were evaluated using a NanoDropTM 2000 ultraviolet-visible spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The extracted DNA was diluted to roughly 50 to 70 ng/μL.

2.3.2. PCR Amplification and Sequence Alignment

The primers and probes for the HRM and Sanger sequencing analyses were designed using Primer Express 3.0.1 software (Thermo Fisher Scientific, Waltham, MA, USA) (Table 1). The 3′ end of the probe was blocked with a C3 spacer and the primers and probes were synthesized by Sangon Biotech Co. Ltd. (Shanghai, China). A total reaction volume of 50 μL was used for the PCR amplification procedure, which included 25 μL of 2 × Phanta Max Master Mix (Vazyme, Nanjing, China), 2 μL each of upstream and downstream primers, 20 μL ddH2O, and 1 μL of a DNA isolate (concentration ranging from 50 to 70 ng/μL). The PCR reaction program was as follows: pre-denaturation at 95 °C for 3 min (1 cycle); 95 °C for 15 s, 60 °C for 15 s, and 72 °C for 3 s (35 cycles); and an extension at 72 °C for 3 min (1 cycle). The PCR products were purified using an EZ-10 Column DNA Purification kit (Sangon Biotech, Shanghai, China) following the manufacturer’s instructions, and the purified DNA fragments were sent to Sangon Biotech Co. Ltd. (Shanghai, China) for Sanger sequencing. The amplified DNA fragments were aligned with the sheep TBXT reference sequence (NCBI accession number: 101114280) using Chromas 2.6.5 software (Technelysium Pty Ltd., SouthBrisbane, QLD, Australia).

2.3.3. HRM Reaction

Using an HRM analysis kit (Tiangen, Beijing, China), the total reaction volume was 10 μL, including 5 μL of 2 × HRM Analysis PreMix (with EvaGreen) with 0.2 μL each of upstream and downstream primers, 10 μM of downstream primers, and 1 μM of upstream primers as well as 0.2 μL of an unlabeled probe (concentration 10 μM), 3.9 μL of ddH2O, and 0.5 μL of DNA amplification using PikoRealTM real-time fluorescence quantitative PCR equipment (Thermo Fisher Scientific, Waltham, MA, USA). The PCR reaction program included pre-denaturation at 95 °C for 2 min and 50 cycles at 95 °C for 10 s and 60 °C for 30 s. The subsequent melting curve analysis consisted of 2 steps: denaturation at 95 °C for 30 s and renaturation at 60 °C for 1 min followed by a continuous fluorescence reading mode from 60 to 95 °C with a ramp rate of 0.02 °C/s. To ensure reliable genotyping results, the reactions were performed in triplicate. PikoRealTM 2.2.250.602 software (Thermo Fisher Scientific, Waltham, MA, USA) was used to evaluate the melting curves, normalize the plots of fluorescence versus temperature, and conduct automated genotyping based on the analysis of the HRM characteristics.

2.3.4. Acquisition of X-ray Images of Sheep Caudal Vertebrae

To observe the number of caudal vertebrae in the Barag sheep and short-tailed sheep, we obtained X-ray images of the caudal vertebrae of sheep of a known genotype using a portable X-ray machine (General Electric Company, Fairfield, CT, USA).

2.3.5. Statistical Analysis

All statistical analyses were performed using SPSS.25 software (IBM Corp., Armonk, NY, USA). Chi-squared tests were used to investigate the genotypic differences between Barag and short-tailed sheep as well as the differences in the number of caudal vertebrae between the homozygous and heterozygous short-tailed sheep.

3. Results

A statistical analysis and a chi-squared analysis were performed on the Sanger sequencing results of 60 Barag sheep and 60 short-tailed sheep (Table 2). In both groups of sheep, we detected a total of three genotypes (C-T/C-T, C-T/G-G, and G-G/G-G; Figure 1). In the Barag sheep, nucleotides 333 and 334 in TBXT were both G. No individuals with c.333G > C or c.334G > T were found. We discovered that the primary genotype in the short-tailed sheep population was c.333G > C; c.334G > T homozygotes (C-T/C-T) accounted for 86.7% of the sheep and the heterozygous genotype (C-T/G-G) accounted for 13.3% of the sheep. Similarly, we found no individuals with G at TBXT nucleotides 333 and 334 among the short-tailed sheep. This result showed that TBXT (c.333G > C or c.334G > T) was strongly associated with the sheep tail traits (p < 0.05) (Table 2).
To investigate if the homozygous short-tailed sheep (C-T/C-T) had fewer or shorter caudal vertebrae than the heterozygous short-tailed sheep (C-T/G-G), we obtained X-ray image phenotypes of 14 short-tailed sheep and 16 Barag sheep. In Figure 2, we show the representative X-ray images of the caudal vertebrae of the short-tailed sheep and Barag sheep. It is worth noting that the ends of the caudal vertebrae of the short-tailed sheep were deformed. The number of caudal vertebrae of the Barag sheep were 9, 10, or 12 (Figure 2b and Table 3) whereas the number of caudal vertebrae in the short-tailed sheep ranged from 4 to 8 (Figure 2a and Table 3). We recorded the number of caudal vertebrae corresponding with the two genotypes of short-tailed sheep and performed a chi-squared test (Table 4). The results showed that the genotype of short-tailed sheep was not associated with the number of caudal vertebrae (p > 0.05). This showed that it is difficult to distinguish homozygous individuals from heterozygous individuals by observing the tail phenotype alone.
Using the HRM method, we genotyped 60 short-tailed sheep and 60 Barag sheep (Figure 3). We added unlabeled probes complementary to the detected DNA fragments using asymmetric PCR to form a single strand of nucleotides [38]. During the PCR amplification, the unlabeled probes competed with and bound to the template DNA sequence of the genotype G-G/G-G. Therefore, the Tm of sample genotype G-G/G-G was the lowest at only 86.84 °C and that of the heterozygous sample C-T/G-G, for which only half the sequence could be combined with the unlabeled probe, was 87.32 °C. The homozygous genotype of the short-tailed sheep, C-T/C-T, did not bind to the unlabeled probes so the Tm was 87.63 °C, which was the highest. The melting curves of the three genotypes (C-T/C-T, C-T/G-G, and G-G/G-G) can be clearly seen in Figure 3. HRM using unlabeled probes and asymmetric PCR increased the Tm difference between the homozygotes and heterozygotes compared with HRM; with smaller amplicons, the accuracy of the genotyping was improved [38,39,40]. To ensure that the HRM genotyping results were accurate, we compared them with the results of Sanger sequencing, which demonstrated that the same genotype was obtained using both methods.

4. Discussion

One of the purposes of the tail of a sheep is to help the sheep retain energy and fat [41]. When food is scarce, fat from the tail can be utilized to keep the body functioning smoothly [14]. Tail docking not only causes a stress reaction in lambs but also increases the risk of wound infection [42,43,44]. Short-tailed sheep can be actively raised by farmers to avoid the previously mentioned aberrant state and enhance animal welfare [9,45,46,47]. Compared with the danger of tail docking [13], producing sheep with a short-tailed phenotype in their natural condition is more in line with current sheep breeding demands. The number and length of caudal vertebrae determine the overall length of the tail and although short-tailed sheep do not have the same amount of fat in their tails as long-tailed sheep, their bodies compensate by increasing intramuscular fat [48] because tail fat can only be deposited on the bones of the caudal vertebrae [9]. The flesh quality of short-tailed lambs is superior to that of long-tailed lambs. Two linked mutations in the TBXT gene, c.333G > C and c.334G > T, have been linked to the short-tailed phenotype in various animals and can be utilized as genetic markers in breeding operations. As a result, utilizing TBXT (c.333G > C or c.334G > T) for homozygous short-tailed sheep breeding is reliable [9,16,17].
Unlike other genotyping studies [49], most of the mutation sites identified by HRM involve only one base pair. However, there are two linked mutation sites to be identified in sheep with the short-tailed phenotype. The Tm generated by two linked mutations cannot be directly increased by HRM utilizing short amplicons. To improve the Tm difference between the homozygous and heterozygous DNA fragments in this experiment, we used unlabeled probes and asymmetric PCR. As a result, the detected fluorescent signal was enhanced. In screening individuals with genetic mutations in flocks, an HRM analysis has a wide range of applications [35,50]. This technology is less costly than TaqMan and Sanger sequencing but has the same precision and applicability. This approach primarily employs nonspecific binding of saturated fluorescent dyes to double-stranded DNA during PCR [51]. Fluorescent colors are removed from the DNA fragments when they are heated and denatured. Using a fluorescent trap, the binding of double-stranded DNA and fluorescent dyes can be observed in real-time during the heating phase [52]. The melting curve is the result of the detection machine monitoring the fluorescent signal in real-time as it changes. When the bases in the amplified fragments are different, the Tm values of the DNA fragments also differ and distinct melting curves can be used to distinguish different genotypes.
Utilizing unlabeled probes and asymmetric PCR in the HRM system can increase Tm differences across the genotypes [31,53]; therefore, we recommend that this method is used because of the small difference in Tm between the homozygotes and heterozygotes for the short-tailed sheep TBXT mutation. When designing unlabeled probes for HRM, we discovered that the accuracy of genotyping could be improved by placing the mutation in the middle of the probe fragment, which increased the sensitivity of detecting both homozygous and heterozygous individuals. Furthermore, we discovered that the concentration of the DNA template in the HRM system should be investigated. According to previous studies [54], the DNA template concentration is an important determinant for the accuracy of genotyping results; however, when the DNA concentration was less than 5 ng/μL in our tests, the genotyping results were inaccurate because multiple nonspecific amplifications occurred during the HRM reaction, reducing the PCR amplification efficiency. The application of molecular identification methods can facilitate a more rapid and accurate determination than a phenotypic identification; therefore, it is important to establish a rapid detection technology for economically important traits in sheep.

5. Conclusions

In conclusion, we developed an HRM approach based on unlabeled probes and asymmetric PCR that could reliably detect TBXT heterozygotes in short-tailed sheep. This method may be used in short-tailed sheep breeding.

6. Patents

An application for a Chinese invention patent has been submitted for this work (application number: 202110056027.2)

Author Contributions

G.Y. performed the research. G.Y. and C.W. analyzed and visualized the data. H.S. and D.W. contributed to the sample and data acquisition on the farms. G.C. and Y.S. designed the study. Y.Y. and X.X. drafted the manuscript. A.D., L.C., T.M., M.L., J.S., X.X. and G.C. reviewed the manuscript. T.H. and X.L. project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 31860689; the Inner Mongolia Autonomous Region Science and Technology Plan of China, grant number 2019GG241; 2020ZD0007; and the Inner Mongolia Autonomous Region Science and Technology Achievements Transformation Project, grant number 2020CG0078.

Institutional Review Board Statement

The animal study was reviewed and approved by the special committee on scientific research and the academic ethics of the Inner Mongolia Agricultural University. (License No. 2020008).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

G.C., G.Y., Y.S., A.D., T.M., C.W., H.S., M.L., D.W., L.C., X.X., T.H., Y.Y., J.S. and X.L. are co-authors of an issued Chinese invention patent (as described above) covering the high-resolution melting curve used in this study. The efficacy of this patented approach has not been demonstrated nor do the authors believe the patent or any future benefits they may receive because of it inappropriately influence any work described in this manuscript. The other authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Mincer, S.T.; Russo, G.A. Substrate use drives the macroevolution of mammalian tail length diversity. Proc. Boil. Sci. 2020, 287, 20192885. [Google Scholar] [CrossRef]
  2. Hickman, G.C. The mammalian tail: A review of functions. Mammal Rev. 1979, 9, 143–157. [Google Scholar] [CrossRef]
  3. Zhu, C.; Cheng, H.; Li, N.; Liu, T.; Ma, Y. Isobaric tags for relative and absolute quantification-based proteomics reveals candidate proteins of fat deposition in chinese indigenous sheep with morphologically different tails. Front. Genet. 2021, 12, 710449. [Google Scholar] [CrossRef] [PubMed]
  4. Abied, A.; Bagadi, A.; Bordbar, F.; Pu, Y.; Augustino, S.M.; Xue, X.; Xing, F.; Gebreselassie, G.; Han, J.-L.; Mwacharo, J.-L.H.J.M.; et al. Genomic diversity, population structure, and signature of selection in five chinese native sheep breeds adapted to extreme environments. Genes 2020, 11, 494. [Google Scholar] [CrossRef]
  5. Moradi, M.H.; Nejati-Javaremi, A.; Moradi-Shahrbabak, M.; Dodds, K.G.; McEwan, J.C. Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition. BMC Genet. 2012, 13, 10. [Google Scholar] [CrossRef] [Green Version]
  6. Rocha, J.; Chen, S.; Beja-Pereira, A. Molecular evidence for fat-tailed sheep domestication. Trop. Anim. Health Prod. 2011, 43, 1237–1243. [Google Scholar] [CrossRef]
  7. Zhu, C.; Li, N.; Cheng, H.; Ma, Y. Genome wide association study for the identification of genes associated with tail fat deposition in Chinese sheep breeds. Biol. Open 2021, 10, bio054932. [Google Scholar] [CrossRef]
  8. Xu, S.-S.; Ren, X.; Yang, G.-L.; Xie, X.-L.; Zhao, Y.-X.; Zhang, M.; Shen, Z.-Q.; Ren, Y.-L.; Gao, L.; Shen, M.; et al. Genome-wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries). Anim. Genet. 2017, 48, 560–569. [Google Scholar] [CrossRef] [Green Version]
  9. Kalds, P.; Luo, Q.; Sun, K.; Zhou, S.; Chen, Y.; Wang, X. Trends towards revealing the genetic architecture of sheep tail patterning: Promising genes and investigatory pathways. Anim. Genet. 2021, 52, 799–812. [Google Scholar] [CrossRef]
  10. Ben Ettoumia, R.; Vernet, J.; Ortigues-Marty, I.; Kraiem, K.; Majdoub-Mathlouthi, L. Effects of metabolizable energy intake on post weaning lamb growth performance, carcass tissue composition and internal fat depend on animal characteristics: A meta-analysis. Meat Sci. 2021, 185, 108719. [Google Scholar] [CrossRef]
  11. Soriano, V.; Phillips, C.; Taconeli, C.; Fragoso, A.; Molento, C. Mind the Gap: Animal Protection Law and Opinion of Sheep Farmers and Lay Citizens Regarding Animal Maltreatment in Sheep Farming in Southern Brazil. Animals 2021, 11, 1903. [Google Scholar] [CrossRef] [PubMed]
  12. McCracken, L.; Waran, N.; Mitchinson, S.; Johnson, C.B. Effect of age at castration on behavioural response to subsequent tail docking in lambs. Vet. Anaesth. Analg. 2010, 37, 375–381. [Google Scholar] [CrossRef]
  13. Small, A.; Fisher, A.; Lee, C.; Colditz, I. Analgesia for Sheep in Commercial Production: Where to Next? Animals 2021, 11, 1127. [Google Scholar] [CrossRef] [PubMed]
  14. Phillips, C.J.C.; Wojciechowska, J.; Meng, J.; Cross, N. Perceptions of the importance of different welfare issues in livestock production. Animals 2009, 3, 1152–1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Abdullah, A.Y.; Kridli, R.T.; Shaker, M.M.; Obeidat, M.D. Investigation of growth and carcass characteristics of pure and crossbred Awassi lambs. Small Rumin. Res. 2010, 94, 167–175. [Google Scholar] [CrossRef]
  16. Han, J.; Yang, M.; Guo, T.; Niu, C.; Liu, J.; Yue, Y.; Yuan, C.; Yang, B. Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep. Anim. Genet. 2019, 50, 772–777. [Google Scholar] [CrossRef] [Green Version]
  17. Zhi, D.; Da, L.; Liu, M.; Cheng, C.; Zhang, Y.; Wang, X.; Li, X.; Tian, Z.; Yang, Y.; He, T.; et al. Whole genome sequencing of hulunbuir short-tailed sheep for identifying candidate genes related to the short-tail phenotype. G3 Genes Genomes Genet. 2018, 8, 377–383. [Google Scholar] [CrossRef] [Green Version]
  18. Ahbara, A.; Bahbahani, H.; Almathen, F.; Al Abri, M.; Agoub, M.O.; Abeba, A.; Kebede, A.; Musa, H.H.; Mastrangelo, S.; Pilla, F.; et al. Genome-Wide Variation, Candidate Regions and Genes Associated with Fat Deposition and Tail Morphology in Ethiopian Indigenous Sheep. Front. Genet. 2019, 9, 699. [Google Scholar] [CrossRef]
  19. Dong, K.; Yang, M.; Han, J.; Ma, Q.; Han, J.; Song, Z.; Luosang, C.; Gorkhali, N.A.; Yang, B.; He, X.; et al. Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep. BMC Genom. 2020, 21, 800. [Google Scholar] [CrossRef]
  20. Li, X.; Yang, J.; Shen, M.; Xie, X.-L.; Liu, G.-J.; Xu, Y.-X.; Lv, F.-H.; Yang, H.; Yang, Y.-L.; Liu, C.-B.; et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits. Nat. Commun. 2020, 11, 2815. [Google Scholar] [CrossRef]
  21. Herrmann, B.G.; Labeit, S.; Poustka, A.; King, T.R.; Lehrach, H. Cloning of the T gene required in mesoderm formation in the mouse. Nature 1990, 343, 617–622. [Google Scholar] [CrossRef]
  22. Lolas, M.; Valenzuela, P.D.T.; Tjian, R.; Liu, Z. Charting Brachyury-mediated developmental pathways during early mouse embryogenesis. Proc. Natl. Acad. Sci. USA 2014, 111, 4478–4483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Schifferl, D.; Scholze-Wittler, M.; Wittler, L.; Veenvliet, J.V.; Koch, F.; Herrmann, B.G. A 37 kb region upstream of brachyury comprising a notochord enhancer is essential for notochord and tail development. Development 2021, 148, 173–182. [Google Scholar] [CrossRef] [PubMed]
  24. Technau, U. Brachyury, the blastopore and the evolution of the mesoderm. Bioessays 2001, 23, 788–794. [Google Scholar] [CrossRef] [PubMed]
  25. Wu, B.; Shao, Y.; Chen, B.; Liu, C.; Xue, Z.; Wu, P.; Li, H. Identification of a novel mouse brachyury (T) allele causing a short tail mutation in mice. Cell Biochem. Biophys. 2010, 58, 129–135. [Google Scholar] [CrossRef]
  26. Buckingham, K.J.; McMillin, M.J.; Brassil, M.M.; Shively, K.M.; Magnaye, K.M.; Cortes, A.; Weinmann, A.S.; Lyons, L.A.; Bamshad, M.J. Multiple mutant T alleles cause haploinsufficiency of Brachyury and short tails in Manx cats. Mamm. Genome 2013, 24, 400–408. [Google Scholar] [CrossRef] [Green Version]
  27. Xu, X.; Sun, X.; Hu, X.-S.; Zhuang, Y.; Liu, Y.-C.; Meng, H.; Miao, L.; Yu, H.; Luo, S.-J. Whole Genome Sequencing Identifies a Missense Mutation in HES7 Associated with Short Tails in Asian Domestic Cats. Sci. Rep. 2016, 6, 31583. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, T.; Gao, H.; Sahana, G.; Zan, Y.; Fan, H.; Liu, J.; Shi, L.; Wang, H.; Du, L.; Wang, L.; et al. Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep. J. Anim. Breed. Genet. 2019, 136, 362–370. [Google Scholar] [CrossRef]
  29. Fan, H.; Hou, Y.; Sahana, G.; Gao, H.; Zhu, C.; Du, L.; Zhao, F.; Wang, L. A transcriptomic study of the tail fat deposition in two types of hulun buir sheep according to tail size and sex. Animals 2019, 9, 655. [Google Scholar] [CrossRef] [Green Version]
  30. Smołucha, G.; Piórkowska, K.; Ropka-Molik, K.; Sikora, J. Use of the HRM method in quick identification of FecXO mutation in highly prolific olkuska sheep. Animals 2020, 10, 844. [Google Scholar] [CrossRef]
  31. Wittwer, C.T. High-resolution DNA melting analysis: Advancements and limitations. Hum. Mutat. 2009, 30, 857–859. [Google Scholar] [CrossRef] [PubMed]
  32. Vossen, R.H.; Aten, E.; Roos, A.; Den Dunnen, J.T. High-Resolution Melting Analysis (HRMA)-More than just sequence variant screening. Hum. Mutat. 2009, 30, 860–866. [Google Scholar] [CrossRef] [PubMed]
  33. Simko, I. High-resolution DNA melting analysis in plant research. Trends Plant Sci. 2016, 21, 528–537. [Google Scholar] [CrossRef] [PubMed]
  34. England, A.; Kheravii, S.; Musigwa, S.; Kumar, A.; Daneshmand, A.; Sharma, N.; Gharib-Naseri, K.; Wu, S. Sexing chickens (Gallus gallus domesticus) with high-resolution melting analysis using feather crude DNA. Poult. Sci. 2021, 100, 100924. [Google Scholar] [CrossRef] [PubMed]
  35. Keikha, M.; Karbalaei, M. High resolution melting assay as a reliable method for diagnosing drug-resistant TB cases: A systematic review and meta-analysis. BMC Infect. Dis. 2021, 21, 989. [Google Scholar] [CrossRef] [PubMed]
  36. Feng, T.; Geng, C.X.; Lang, X.Z.; Chu, M.X.; Cao, G.L.; Di, R.; Fang, L.; Chen, H.Q.; Liu, X.L.; Li, N. Polymorphisms of caprine GDF9 gene and their association with litter size in Jining Grey goats. Mol. Biol. Rep. 2010, 38, 5189–5197. [Google Scholar] [CrossRef]
  37. Alkafajy, A.; Al-Karagoly, H.; Brujeni, G.N. Comparison of cattle BoLA-DRB3 typing by PCR-RFLP, direct sequencing, and high-resolution DNA melting curve analysis. Vet. Res. Forum 2020, 11, 21–26. [Google Scholar]
  38. Luo, R.; Zhang, X.; Wang, L.; Zhang, L.; Li, G.; Zheng, Z. GLIS1, a potential candidate gene affect fat deposition in sheep tail. Mol. Biol. Rep. 2021, 48, 4925–4931. [Google Scholar] [CrossRef]
  39. Seipp, M.T.; Durtschi, J.D.; Liew, M.A.; Williams, J.; Damjanovich, K.; Pont-Kingdon, G.; Lyon, E.; Voelkerding, K.V.; Wittwer, C.T. Unlabeled oligonucleotides as internal temperature controls for genotyping by amplicon melting. J. Mol. Diagn. 2007, 9, 284–289. [Google Scholar] [CrossRef] [Green Version]
  40. Liew, M.; Seipp, M.; Durtschi, J.; Margraf, R.L.; Dames, S.; Erali, M.; Voelkerding, K.; Wittwer, C. Closed-tube SNP genotyping without labeled probes/a comparison between unlabeled probe and amplicon melting. Am. J. Clin. Pathol. 2007, 127, 341–348. [Google Scholar] [CrossRef]
  41. Erali, M.; Palais, R.; Wittwer, C. SNP Genotyping by unlabeled probe melting analysis. Methods Mol. Biol. 2008, 429, 199–206. [Google Scholar] [CrossRef] [PubMed]
  42. Rocha, J.L.; Godinho, R.; Brito, J.C.; Nielsen, R. Life in deserts: The genetic basis of mammalian desert adaptation. Trends Ecol. Evol. 2021, 36, 637–650. [Google Scholar] [CrossRef] [PubMed]
  43. Zhao, F.; Deng, T.; Shi, L.; Wang, W.; Zhang, Q.; Du, L.; Wang, L. Genomic scan for selection signature reveals fat deposition in chinese indigenous sheep with extreme tail types. Animals 2020, 10, 773. [Google Scholar] [CrossRef] [PubMed]
  44. He, X.; Wu, R.; Yun, Y.; Qin, X.; Chen, L.; Han, Y.; Wu, J.; Sha, L.; Borjigin, G. Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep. Food Sci. Nutr. 2021, 9, 5722–5734. [Google Scholar] [CrossRef]
  45. Small, A.; Fetiveau, M.; Smith, R.; Colditz, I. Three studies evaluating the potential for lidocaine, bupivacaine or procaine to reduce pain-related behaviors following ring castration and/or tail docking in lambs. Animals 2021, 11, 3583. [Google Scholar] [CrossRef]
  46. Yousefi, A.R.; Kohram, H.; Shahneh, A.Z.; Nik-Khah, A.; Campbell, A.W. Comparison of the meat quality and fatty acid composition of traditional fat-tailed (Chall) and tailed (Zel) Iranian sheep breeds. Meat Sci. 2012, 92, 417–422. [Google Scholar] [CrossRef]
  47. Small, A.; Marini, D.; Colditz, I. Local anesthetic delivered with a dual action ring and injection applicator reduces the acute pain response of lambs during tail docking. Animals 2021, 11, 2242. [Google Scholar] [CrossRef]
  48. Han, J.; Guo, T.; Yue, Y.; Lu, Z.; Liu, J.; Yuan, C.; Niu, C.; Yang, M.; Yang, B. Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs. PLoS ONE 2021, 16, e0246279. [Google Scholar] [CrossRef]
  49. Obšteter, J.; Jenko, J.; Gorjanc, G. Genomic selection for any dairy breeding program via optimized investment in phenotyping and genotyping. Front. Genet. 2021, 12, 637017. [Google Scholar] [CrossRef]
  50. Słomka, M.; Sobalska-Kwapis, M.; Wachulec, M.; Bartosz, G.; Strapagiel, D. High Resolution Melting (HRM) for high-throughput genotyping—Limitations and caveats in practical case studies. Int. J. Mol. Sci. 2017, 18, 2316. [Google Scholar] [CrossRef] [Green Version]
  51. Smith, B.L.; Lu, C.; Bremer, J.R.A. High-resolution melting analysis (HRMA): A highly sensitive inexpensive genotyping alternative for population studies. Mol. Ecol. Resour. 2010, 10, 193–196. [Google Scholar] [CrossRef] [PubMed]
  52. Froehlich, T.; Geulen, O. Hybridization Probe Pairs and Single-Labeled Probes: An Alternative Approach for Genotyping and Quantification. Methods Mol. Biol. 2008, 429, 117–133. [Google Scholar] [CrossRef] [PubMed]
  53. Ririe, K.M.; Rasmussen, R.P.; Wittwer, C.T. Product Differentiation by Analysis of DNA Melting Curves during the Polymerase Chain Reaction. Anal. Biochem. 1997, 245, 154–160. [Google Scholar] [CrossRef] [PubMed]
  54. Lu, C.; Saint-Pierre, C.; Gasparutto, D.; Roupioz, Y.; Ravelet, C.; Peyrin, E.; Buhot, A. Melting Curve Analysis of Aptachains: Adenosine Detection with Internal Calibration. Biosensors 2021, 11, 112. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Three genotypes of TBXT in Barag sheep and short-tailed sheep. Blue arrows point to two polymorphic sites. (a) Fragments of chromatograms with genotype C-T/C-T. (b) Fragments of chromatograms for heterozygous genotype C-T/G-G. (c) Fragments of chromatograms with genotype G-G/G-G.
Figure 1. Three genotypes of TBXT in Barag sheep and short-tailed sheep. Blue arrows point to two polymorphic sites. (a) Fragments of chromatograms with genotype C-T/C-T. (b) Fragments of chromatograms for heterozygous genotype C-T/G-G. (c) Fragments of chromatograms with genotype G-G/G-G.
Animals 12 00792 g001
Figure 2. X-ray images of the caudal vertebrae of Barag sheep and short-tailed sheep. (a) Representative X-ray images of the caudal vertebrae of short-tailed sheep. (b) Representative X-ray images of the three types of caudal vertebrae of Barag sheep.
Figure 2. X-ray images of the caudal vertebrae of Barag sheep and short-tailed sheep. (a) Representative X-ray images of the caudal vertebrae of short-tailed sheep. (b) Representative X-ray images of the three types of caudal vertebrae of Barag sheep.
Animals 12 00792 g002
Figure 3. HRM genotyping results. (a) Typical shapes of the melting curves corresponding with the three genotypes (C-T/C-T, C-T/G-G, and G-G/G-G). The color of the curve is automatically assigned by PikoRealTM software. The vertical axis represents the change in fluorescence intensity and the abscissa represents the Tm value. The value of the abscissa corresponding with the main peak of each curve is the Tm value of the amplified fragment. (b) The melting curve image generated during the genotyping of Barag sheep and short-tailed sheep (a total of 32 samples are shown in the image; each sample was run in triplicate).
Figure 3. HRM genotyping results. (a) Typical shapes of the melting curves corresponding with the three genotypes (C-T/C-T, C-T/G-G, and G-G/G-G). The color of the curve is automatically assigned by PikoRealTM software. The vertical axis represents the change in fluorescence intensity and the abscissa represents the Tm value. The value of the abscissa corresponding with the main peak of each curve is the Tm value of the amplified fragment. (b) The melting curve image generated during the genotyping of Barag sheep and short-tailed sheep (a total of 32 samples are shown in the image; each sample was run in triplicate).
Animals 12 00792 g003
Table 1. Primers used in HRM and Sanger sequencing.
Table 1. Primers used in HRM and Sanger sequencing.
NameSequence (5′-3′)Analysis Type PCR Product Length
TBXT-FTGCGCCCCTTCCTTTTCAGHRM/PCR203 bp
TBXT-RGGGGGAGTCGGGGTGGATGTAGHRM/PCR
TBXT-ProbeGCTTGCCCCAGGGCACCCAHRM
Table 2. Genotyping statistics for Barag sheep and short-tailed sheep.
Table 2. Genotyping statistics for Barag sheep and short-tailed sheep.
BreedGenotype (TBXT c.333G > C; c.334G > T)X2p-Value
G-G/G-GC-T/G-GC-T/C-T
Barag sheep100%00154.872.07 × 10−35
Short-tailed sheep013.3%86.7%
Table 3. Statistics on the number of caudal vertebrae in Barag sheep and short-tailed sheep.
Table 3. Statistics on the number of caudal vertebrae in Barag sheep and short-tailed sheep.
Number of Caudal VertebraeTotal
456789101112
Number of short-tailed sheep24161000014
Number of Barag sheep00000680216
Table 4. Statistics on the number and genotype of caudal vertebrae in short-tailed sheep.
Table 4. Statistics on the number and genotype of caudal vertebrae in short-tailed sheep.
Number of Caudal VertebraeTotalX2p-Value
45678
Number of homozygotes1213072.2820451
Number of heterozygotes120317
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, G.; Wang, C.; Su, H.; Wang, D.; Dou, A.; Chen, L.; Ma, T.; Liu, M.; Su, J.; Xu, X.; et al. Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes. Animals 2022, 12, 792. https://doi.org/10.3390/ani12060792

AMA Style

Yang G, Wang C, Su H, Wang D, Dou A, Chen L, Ma T, Liu M, Su J, Xu X, et al. Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes. Animals. 2022; 12(6):792. https://doi.org/10.3390/ani12060792

Chicago/Turabian Style

Yang, Guang, Caiyun Wang, Hong Su, Daqing Wang, Aolie Dou, Lu Chen, Teng Ma, Moning Liu, Jie Su, Xiaojing Xu, and et al. 2022. "Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes" Animals 12, no. 6: 792. https://doi.org/10.3390/ani12060792

APA Style

Yang, G., Wang, C., Su, H., Wang, D., Dou, A., Chen, L., Ma, T., Liu, M., Su, J., Xu, X., Yang, Y., He, T., Li, X., Song, Y., & Cao, G. (2022). Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes. Animals, 12(6), 792. https://doi.org/10.3390/ani12060792

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop