Front. Vet. Sci.Frontiers in Veterinary ScienceFront. Vet. Sci.2297-1769Frontiers Media S.A.10.3389/fvets.2023.1134903Veterinary ScienceOriginal ResearchThe ovine HIAT1 gene: mRNA expression, InDel mutations, and growth trait associationsLuoYunyun1AkhatayevaZhanerke1MaoCui12JiangFugui2GuoZhengang3XuHongwei4*LanXianyong1*1Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China2Shandong Key Lab of Animal Disease Control and Breeding, Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, Jinan, China3Bijie Animal Husbandry and Veterinary Science Research Institute, Bijie, China4College of Life Science and Engineering, Northwest Minzu University, Lanzhou, China
Edited by: Ibrar Muhammad Khan, Fuyang Normal University, China
Reviewed by: Cemal Ün, Ege University, Türkiye; Muhammad Zahoor, University of Agriculture, Dera Ismail Khan, Pakistan
*Correspondence: Hongwei Xu xuhongwei@xbmu.edu.cnXianyong Lan lanxianyong79@nwsuaf.edu.cn
This article was submitted to Livestock Genomics, a section of the journal Frontiers in Veterinary Science
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
The hippocampal abundant transcript 1 (HIAT1) gene, also known as major facilitator superfamily domain-containing 14A (MFSD14A), encodes for a transmembrane transporter protein and has been previously shown to be associated with milk production in buffalo and sheep breeds, as well as growth traits in chicken and goats. However, tissue level distribution of the ovine HIAT1 gene, as well as its effect on body morphometric traits in sheep, has yet to be studied.
Methods
The HIAT1 mRNA expression profile of Lanzhou fat-tailed (LFT) sheep was determined by quantitative real-time PCR (qPCR). A total of 1498 sheep of three indigenous Chinese sheep breeds were PCR-genotyped for polymorphisms of HIAT1 gene. Student's t-test was used to observe the association between the genotype and sheep morphometric traits.
Results
HIAT1 was widely expressed in all examined tissues, and was particularly abundant in the testis of male LFT sheep. Additionally, a 9-bp insertion mutation (rs1089950828) located within the 5'-upstream region of HIAT1 was investigated in Luxi black-headed (LXBH) sheep and Guiqian semi-fine wool (GSFW) sheep. The wildtype allele frequency 'D' was found to be more prevalent than that of the mutant allele ‘I'. Furthermore, low genetic diversity was confirmed in all sampled sheep populations. Subsequent association analyses indicated an association between the 9-bp InDel mutation of interest and the morphometric traits of LXBH and GSFW sheep. Furthermore, yearling ewes with a heterozygous genotype (ID) demonstrated smaller body sizes, while yearling rams and adult ewes with the heterozygous genotype were found to have overall better growth performance.
Conclusion
These findings imply that functional InDel polymorphism (rs1089950828) has the potential to be utilized for marker-assisted selection (MAS) of growth traits in domestic Chinese sheep populations.
The final grown performance of sheep and other domesticated animals is affected by multiple factors including environmental, nutritional, and genetic influences (1). Using high-throughput whole genome sequencing, RNA sequencing, and molecular biology techniques, a number of candidate functional genes and genetic variants linked to important growth traits in livestock were found and used for the purposes of molecular-assisted and omics-assisted breeding (2–6). Among the numerous genome variants including SNPs (single nucleotide polymorphisms), CNVs (copy number variations), and other genetic markers, numerous studies have demonstrated that InDels (insertion/deletion) function by regulating gene translation or altering a protein and are related to differences in the characteristics of individual animals (7–9).
The hippocampal abundant transcript 1 (HIAT1) gene is also known as the major facilitator superfamily domain 14A (MFSD14A) gene (10). The HIAT1 protein contains a conserved sugar transporter motif as well as 12 transmembrane alpha helices (TMs) linked by hydrophilic loops, and belongs to the major facilitator superfamily (MFS) of secondary protein transporters (11). The MFS protein superfamily is the largest secondary transporter family present in a broad variety of species from archaea to mammals and plays an essential role in the exchange of cellular materials and energy metabolism, which facilitates the transport of many substrates across both cytoplasmic and internal cell membranes (12). The expression of both MFSD1 and MFSD3 genes has previously been shown to be significantly different in both the cerebellum and brain of mice after they were subjected to starvation and high-fat diets. This suggests that these genes are regulated at the nutrient level and may be related to nutrient intake and transport (13). Two members of the MFS gene family, MFSD6 and MFSD10 were found to be related to the energy metabolism of cells. Expression levels of these two genes were found to be down-regulated when the cellular energy consumption rate increased or when cellular energy was significantly depleted (14). These findings suggest that HIAT1 may also be related to fatty acid uptake and transport.
Chinese scholars have found that HIAT1 affects buffalo milk fat synthesis mainly through regulation of cellular membrane transport functions (15). Moreover, studies have shown that the HIAT1 protein has the capacity to transport solute from the blood stream that is required for spermiogenesis, thus disruption of this gene has been shown to cause round-headed sperm and sterility in male mice (16). HIAT1 has previously been shown to be commonly expressed in multiple tissues types and highly expressed in testis, suggesting that it has extensive biological functions (17). HIAT1 is expressed in brown adipose tissue and is negatively correlated with the process in which glucose is converted to fat (18). Additionally, polymorphisms in HIAT1 have been shown to affect milk traits in buffalo as well as growth in goats (15, 19). However, to date, there are few studies that have reported on the function of HIAT1 in indigenous Chinese sheep populations.
In order to further improve the production features of native Chinese sheep populations, the current study aimed to detect HIAT1 expression features in Lanzhou fat-tailed (LFT) sheep, explore InDel variations of the HIAT1 gene in three indigenous Chinese sheep breeds, which included Luxi black-headed (LXBH) sheep from east China, Guiqian semi-fine wool (GSFW) sheep from southwest China, and LFT sheep from northwest China, as well as conduct association analyses between sheep genotypes and growth traits to provide a basis for use of the HIAT1 gene for sheep selection and breeding.
2. Materials and methods2.1. Animals and ethics statement
Twelve tissue samples were collected from three 6-month-old male Lanzhou fat-tailed sheep and include heart, liver, spleen, lung, kidney, testis, longissimus dorsi muscle, tail adipose, perirenal fat, subcutaneous fat, small intestine, and rennet stomach (Yongjing county, Gansu Province). A total of 1498 indigenous Chinese sheep were used for DNA experiments and included Luxi black-headed sheep (n = 616) raised in Liaocheng city, Shandong province; Guiqian semi-fine wool sheep (n = 824) raised in Bijie city, Guizhou province; and Lanzhou fat-tailed sheep (n = 58) provided by Gansu Ruilin Science and Technology Breeding Company (Yongjing county, Gansu province). All tissue samples were frozen in liquid nitrogen and stored at −80°C for RNA and DNA isolation. Relevant body morphometric traits were recorded, including body weight, chest circumference, height at hip cross, body height, body length, cannon circumference, hip width, abdominal circumference, chest depth, and chest width. Animal experiments were conducted in line with the guidelines for the care and use of animals and approved by the Animal Ethics Committee of Northwest A&F University.
2.2. DNA isolation, primer design, and InDel genotyping
Genomic DNA was isolated following previously established protocols (20, 21). Total RNA and genomic DNA concentrations were determined using Nanodrop 2000. DNA was diluted to a concentration of 20 ng/μL and stored at −40°C. Thirty DNA samples were randomly selected from each breed to construct three genomic DNA pools for investigation of HIAT1 variations.
Three upstream gene variants were screened from the Ensembl database, and then several primers (Table 1) were designed using Primer Premier 5 software and synthesized by Sangon Biotech (Xi'an, China) based on reference sequences in ovine HIAT1 (NC 056054.1). Polymorphic fragments were amplified using the polymerase chain reaction (PCR) method and utilized a touch-down program and reaction volumes as described in previous studies (8, 22). PCR products were then separated by electrophoresis on a 3.5% agarose gel stained with GoldView. The genotype of each individual sheep was analyzed.
Primers used for InDel detection and qPCR analysis of ovine HIAT1.
2.3. RNA isolation, cDNA synthesis, and quantitative real-time PCR
Total RNA was extracted from tissue samples using TRIzol reagent, and cDNA was synthesized according to manufactures' instructions using an Evo M-MLV RT Kit with gDNA Clean for qPCR II (Accurate biology Co., Ltd, Changsha, China). A LightCycler 96 real-time PCR system (Roche, Switzerland) was used to run the qPCR program. qPCR reactions (10 μL) contained 5 μL 2 × ChamQ SYBR qPCR Master Mix (Vazyme, Nanjing, China), 0.2 μL of each primer, and 4.6 μL cDNA (1/100 dilution). qPCR conditions were as follows: 95°C for 180 s, followed by 40 amplification cycles lasting 40 s each (95°C for 10 s and 60°C for 30 s). Gene expression was quantitated using the 2−ΔΔCt method.
2.4. Bioinformatics and statistical analysis
Differing amino acid sequences alignment were made using DNAMAN 6.0 software. The percent similarity of amino acid sequences was analyzed by NCBI-Blast online software. Using AliBaba2.1 (http://gene-regulation.com/pub/programs/alibaba2/), predictions of transcription factor binding to the mutant ovine HIAT1 gene sequence were generated (23). Population genetic indices, as well as allelic and genotypic frequencies, were determined. χ2 test was performed to determine whether the gene variants are within Hardy-Weinberg equilibrium (HWE) and analyze the difference in genotype distributions amongst each sheep population. The statistical model used was Yij = μ + Gi + εij, where Yij = observed growth traits, μ = population mean, Gi = fixed effect of the genotype, εij = random error. If there were only two genotypes, student's t-test was performed to determine the effect of different genotypes on economic traits.
3. Results3.1. Multiple sequence alignment and mRNA expression
The results of NCBI-BLASTP sequence analysis showed that the HIAT1 in Ovis aries (Genbank: XP_004002278.1) shares 100%, 100%, 99.80%, 99.80% and 99.59% similarity with that of Capra hircus (XP_017901305.1), Bos taurus (NP_001095508.1), Mus musculus (NP_032272.2), Sus scrofa (XP_013852805.2) and Homo sapiens (NP_149044.2), respectively (Figure 1). Additionally, qPCR results revealed that HIAT1 was expressed in most sheep tissues, especially in the testis of male LFT sheep (Figure 2), which coincided with expression and distribution in human tissues.
Amino acid sequence alignment of HIAT1 in different mammalian species. The red boxes denote the parts of the sequence that are unable to be matched with the Ovis aries sequence. Note: • denotes the same amino acid.
HIAT1 mRNA expression patterns detected by qPCR in Lanzhou fat-tailed sheep.
3.2. Detection and genetic parameter analysis of InDel polymorphisms within the <italic>HIAT1</italic> gene
Electrophoresis patterning and sequence mapping of the P3-ins-9bp InDel showed polymorphisms in the 5'-upstream region of HIAT1 for both Luxi black-headed sheep and Guiqian semi-fine wool sheep. There were no detected HIAT1 polymorphisms found in Lanzhou fat-tailed sheep (Figure 3). Genotyping analyses showed that the HIAT1 P3-ins-9bp locus contained three genotypes (the homozygous insertion type, II, 148 bp; the homozygous deletion type, DD, 139 bp; and the heterozygous mutation type, ID, 148 bp and 139 bp) for both Luxi black-headed sheep and Guiqian semi-fine wool sheep; however, only one genotype (DD) was discovered in the Lanzhou fat-tailed breed.
Electrophoresis gel (A) and sequence chromatograms (B) of the P3-ins-9bp InDel within the ovine HIAT1 gene.
The genotypic frequency, allelic frequency, and population genetic indices were subsequently calculated and shown in Table 2. The DD genotype frequency was found to be much higher than both ID and II genotypes. Amongst all sheep breeds tested, the P3-ins-9bp locus was found to have low genetic polymorphism (PIC < 0.25) while maintaining Hardy-Weinberg equilibrium (P > 0.05). χ2 testing revealed that there were significant differences both in the genotypic and allelic frequency distributions of the P3-ins-9bp locus amongst all three populations (P < 0.05) (Table 3).
Genotypic frequencies, allelic frequencies, and genetic diversity of the P3-ins-9bp InDel.
Chi-square test of HIAT1 P3-ins-9bp genotype distribution in the three sheep breeds.
Breeds
LXBH
GSFW
LFT
LXBH
7e-6**
8.6e-5**
GSFW
1.4e-5**
0.013*
LFT
4.24e-4**
0.038*
Results below the diagonal line represent the Chi-square test of genotype frequency for each sheep breed; Results above the diagonal line represent the Chi-square test of allele frequency for each sheep breed.
P < 0.05,
P < 0.01.
3.3. Association between the P3-ins-9bp InDel and growth traits
The associations between the P3-ins-9bp InDel of the ovine HIAT1 gene and growth performance were analyzed (Table 4). For LXBH sheep, t-tests results confirmed that the P3-ins-9bp InDel had no remarkable influence on growth traits for lambs (P > 0.05) (not shown). Yearling rams with the homozygous DD genotype were found to have lower chest depth, body weight, cannon circumference, and body height (P < 0.05). Yearling ewes with the homozygous DD genotype had significantly higher body weight, chest depth, chest circumference, height at the hip cross and abdominal circumference than those with the heterozygous ID genotype (P < 0.05). Additionally, adult ewes exhibited significant differences in height at the hip cross, body height, and cannon circumference (P < 0.01), with the heterozygous ID genotype being the dominant individual. For GSFW sheep, adult ewes with the heterozygous ID genotype showed superior performance compared to those with the homozygous DD genotype with regards to body weight, body length, chest width, and chest circumference.
Association between the HIAT1 P3-ins-9bp locus and sheep growth traits.
Breeds
Types (n)
Traits
Observed genotypes (LSM ±SE)
P-values
LXBH
Yearling rams
ID (9)
DD (70)
(79)
Body weight (kg)
68.78 ± 7.00
54.47 ± 2.02
0.023
Body height (cm)
70.67 ± 1.76
66.53 ± 0.68
0.041
Chest depth (cm)
31.11 ± 1.15
27.63 ± 0.54
0.029
Cannon circumference (cm)
10.22 ± 0.21
9.54 ± 0.11
0.036
Yearling ewes
ID (16)
DD (82)
(98)
Body weight (kg)
39.99 ± 3.10
48.92 ± 1.52
0.018
Height at the hip cross (cm)
64.38 ± 1.19
66.94 ± 0.51
0.044
Chest depth (cm)
24.60 ± 0.83
27.47 ± 0.45
0.010
Chest circumference (cm)
80.06 ± 2.20
86.87 ± 1.05
0.009
Abdominal circumference (cm)
100.81 ± 3.07
108.52 ± 1.21
0.013
Adult ewes
ID (15)
DD (58)
(73)
Body height (cm)
71.13 ± 1.31
68.05 ± 0.47
0.009
Height at the hip cross (cm)
72.13 ± 1.18
69.16 ± 0.45
0.006
Cannon circumference (cm)
10.20 ± 0.17
9.33 ± 0.12
0.001
GSFW
Adult ewes
ID (19)
DD (155)
(174)
Body weight (kg)
55.26 ± 2.00
49.81 ± 0.50
0.001
Body length (cm)
82.68 ± 1.07
80.69 ± 0.31
0.037
Chest width (cm)
30.05 ± 3.01
28.41 ± 0.18
0.005
Chest circumference (cm)
99.63 ± 1.75
95.17 ± 0.49
0.004
4. Discussion
Recently, a study by Liu et al. (24) identified a quantitative trait locus (QTL) that affects milk traits in buffalo breed and includes SLC35A3 (Solute carrier family 35 member A3) and HIAT1 (24). The SLC35A3 gene is located adjacent to the HIAT1 gene and is involved in the transport of glucose, vitamins and other substances (25). A genome-wide association study (GWAS) on milk production traits determined that SLC35A3 may also play a significant role in sheep milk traits (26). Moreover, a missense mutation in bovine SLC35A3 has been found to lead to complex vertebral malformation (25). A recent study by Ye (15) found that HIAT1 regulates milk fat synthesis in mammary epithelial cells through the PPAR signaling pathway, and that two SNPs within HIAT1 are significantly associated with both milk fat and milk protein rates in buffalo (15). A more surprising finding is that both HIAT1 and SLC35A3 have the potential to affect both carcass and growth traits in chickens (27). In this study, we found that HIAT1 has high amino acid sequence conservation, and that its tissue expression pattern is similar to that of humans. These results suggest that HIAT1 might have a similar function amongst multiple types of mammals. Moreover, HIAT1 mRNA expression was detected in perirenal fat, subcutaneous fat, and the livers of LFT sheep, suggesting that HIAT1 may be related to fat deposition in sheep. These finding suggest HIAT1 is a key candidate gene for regulating important economic traits in livestock.
Body morphology traits reflect the growth and development of different parts of the body and are also important economic and breeding indicators for both meat and meat-wool sheep breeds (28, 29). Existing studies have showed that the HIAT1 gene influenced both growth and meat production traits in goats, and that a 15-bp insertion in HIAT1 was associated with body morphology in Shaanbei white cashmere goats (19, 30). In this study, a 9-bp insertion mutation located in the 5'-upstream region of the ovine HIAT1 gene (rs1089950828 -/GTCCAGTGG) was elucidated. It is well known that variations in the 5'-upstream region of genes has the potential to affect gene transcription and mRNA expression levels through changing transcription factor binding sites (TFBSs) and can also influence an individual's phenotype (31, 32). Thus, this research studied the association between the P3-ins-9bp InDel and sheep growth traits. Furthermore, the P3-ins-9bp InDel was strongly associated with body measurement traits in both LXBH and GSFW sheep based on association analysis. For yearling LXBH sheep, ewes with a heterozygous ID genotype demonstrated overall lower body size, while rams with a heterozygous ID genotype had better growth performance. These results indicate that the 9-bp insertion mutation may have a negative regulatory effect on both growth and development of LXBH yearling ewes as well as promote growth in yearling rams. Additionally, a heterozygous ID genotype in ewes was found to inhibit body size development in yearling sheep as well as promote growth in adult sheep. This may be due to the effect of the 9-bp InDel on the differential expression of HIAT1 at different growth stages throughout the sheep lifecycle (33). Furthermore, the PIC value is an indicator of molecular marker quality in genetic studies. Our InDel locus of interest showed low polymorphism (PIC < 0.25) amongst all tested sheep according to the PIC value, indicating that there was sufficient space in which artificial selection can take place (34, 35).
Next, we searched for the putative transcription factor binding sites that contained our InDel of interest. The P3-ins-9bp locus was found 624-bp upstream of the HIAT1 transcriptional start site (TSS), which created an additional binding site for the transcription factor Sp1 (Specificity protein 1) falling upstream of the four existing Sp1 binding sites (Figure 4). SP1 is a transcription factor of the SP/KLF family. The regulation of SP1 with regards to target genes is achieved by binding of the DNA binding domain to GC-rich motifs in the promoter region of the target gene, thereby regulating cellular processes such as autophagy, apoptosis, proliferation, differentiation and angiogenesis (36, 37). Chen et al. (38) found that Sp1 binding sites in the core promoter region are essential for positive regulation of FGF21 (Fibroblast growth factor 21) through gene transcription in both hepatic and adipose tissue (38). Furthermore, a recent study found that an SNP g.133A>C located within the SCD (Stearoyl CoA desaturase) promoter generated overall higher Sp1 binding affinity to the SCD promoter, which consequently affected gene expression (39). Moreover, a T > C mutation in the promoter region of IGF1 (Insulin-like growth factor 1) was significantly associated with the litter size of Yunshang black goats. This mutation was found to create a new binding site for SP1, thereby promoting goat granulosa cell proliferation by regulating the expression of IGF1 (40). Therefore, the interaction between Sp1 and HIAT1 leads to changes in gene expression at both the mRNA and protein levels, which in turn has the potential to affect the growth traits of sheep. In future studies this interaction is worthy of further in-depth research and exploration. Furthermore, the molecular mechanism of action and associated allelic consequences should be further validated in future studies.
Transcription factor binding sites on the 9-bp InDel locus.
5. Conclusion
This study investigated both the mRNA expression profile and genetic variations of the ovine HIAT1 gene. Our results demonstrated a high mRNA expression level of HIAT1 in the testis of Lanzhou fat-tailed sheep. Moreover, a 9-bp insertion mutation (rs1089950828) in the 5'-upstream region of ovine HIAT1 was detected. The frequency of the mutant allele “I” was found to be low. Furthermore, the 9-bp InDel locus was found to be closely associated with growth traits of LXBH and GSFW sheep, indicating that this InDel has the potential to be used as a genetic marker for assisted selection programs in domestic sheep.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
All animal experiments were conducted in accordance with the the China national standard of Guidelines on Welfare and Ethical Review for Laboratory Animals (GB/T 35892-2018). Our study was approved by Institutional Animal Care and Use Committee (IACUC) of Northwest A&F University.
Author contributions
YL, XL, and HX came up with idea and revised the manuscript. YL wrote the manuscript and performed the experiments. CM, FJ, ZG, and HX collected the sheep samples and isolated the genomic DNA. YL and ZA analyzed the data. All authors approved the final version of the manuscript for submission, contributed to the article, and approved the submitted version.
Funding
This work was supported by the National Natural and Science Foundation of China (32060741).
We thank Shaanxi Key Laboratory of Molecular Biology for Agriculture and the Life Science Research Core Services (LSRCS) of Northwest A&F University (Northern Campus) for their cooperation and support. We would like to thank the staff of Shandong Academy of Agricultural Sciences, Guizhou Xinwumeng Ecological Animal Husbandry Development Co., Ltd., and Gansu Ruilin Science and Technology Breeding Company for providing samples and phenotypic traits data for us.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesLi ZH LiHZhangHWangSZWangQGWangYX. Identification of a single nucleotide polymorphism of the insulin-like growth factor binding protein 2 gene and its association with growth and body composition traits in the chicken. J Anim Sci. (2006) 84:2902–6. 10.2527/jas.2006-14417032782BakhtiarizadehMRAlamoutiAA. RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep. Sci Rep. (2020) 10:13525. 10.1038/s41598-020-70527-832782325GuBSunRFangXZhangJZhaoZHuangD. Genome-wide association study of body conformation traits by whole genome sequencing in Dazu Black Goats. Animals. (2022) 12:548. 10.3390/ani1205054835268118Salek ArdestaniSJafarikiaMSargolzaeiMSullivanBMiarY. Genomic prediction of average daily gain, back-fat thickness, and loin muscle depth using different genomic tools in Canadian Swine populations. Front Genet. (2021) 12:665344. 10.3389/fgene.2021.66534434149806WangKLiuXQiTHuiYYanHQuL. Whole-genome sequencing to identify candidate genes for litter size and to uncover the variant function in goats (Capra hircus). Genomics. (2021) 113:142–50. 10.1016/j.ygeno.2020.11.02433276007RenFYuSChenRLvXYPanCY. Identification of a novel 12-bp insertion/deletion (indel) of iPS-related Oct4 gene and its association with reproductive traits in male piglets. Anim Reprod Sci. (2017) 178:55–60. 10.1016/j.anireprosci.2017.01.00928139300GamarraDAldaiNArakawaAde PancorboMMTaniguchiM. Effect of a genetic polymorphism in SREBP1 on fatty acid composition and related gene expression in subcutaneous fat tissue of beef cattle breeds. Anim Sci J. (2021) 92:e13521. 10.1111/asj.1352133554418WangZZhangXJiangEYanHZhuHChenH. InDels within caprine IGF2BP1 intron 2 and the 3'-untranslated regions are associated with goat growth traits. Anim Genet. (2020) 51:117–21. 10.1111/age.1287131625179WangSHLiuSRYuanTTSunXZ. Genetic effects of FTO gene insertion/deletion (InDel) on fat-tail measurements and growth traits in Tong sheep. Anim Biotechnol. (2021) 32:229–39. 10.1080/10495398.2019.168037931642366FredrikssonRNordströmKJVStephanssonOHägglundMGASchiöthHB. The solute carrier (SLC) complement of the human genome: phylogenetic classification reveals four major families. FEBS Lett. (2008) 582:3811–6. 10.1016/j.febslet.2008.10.01618948099CollinsSMartinTLSurwitRSRobidouxJ. Genetic vulnerability to diet-induced obesity in the C57BL/6J mouse: physiological and molecular characteristics. Physiol. Behav. (2004) 81:243–8. 10.1016/j.physbeh.2004.02.00615159170DittamiSMBarbeyronTBoyenCCambefortJColletGDelageL. Genome and metabolic network of “Candidatus Phaeomarinobacter ectocarpi” Ec32, a new candidate genus of Alphaproteobacteria frequently associated with brown algae. Front Genet. (2014) 5:241. 10.3389/fgene.2014.0024125120558PerlandEHellstenSVLekholmEErikssonMMArapiVFredrikssonR. The novel membrane-bound proteins MFSD1 and MFSD3 are putative SLC transporters affected by altered nutrient Intake. J Mol Neurosci. (2017) 61:199–214. 10.1007/s12031-016-0867-827981419BagchiSPerlandEHosseiniKLundgrenJAl-WalaiNKhederS. Probable role for major facilitator superfamily domain containing 6 (MFSD6) in the brain during variable energy consumption. Int J Neurosci. (2020) 130:476–89. 10.1080/00207454.2019.169402031906755YeT. Identification of Milk Fat Regulation Gene MFSD14A and Its Mechanism in Buffalo. Wuhan: Huazhong Agricultural University. (2022).DoranJWaltersCKyleVWoodingPHammett-BurkeRColledgeWH. Mfsd14a (Hiat1) gene disruption causes globozoospermia and infertility in male mice. Reproduction. (2016) 152:91–9. 10.1530/REP-15-055727107036SreedharanSStephanssonOSchiothHBFredrikssonR. Long evolutionary conservation and considerable tissue specificity of several atypical solute carrier transporters. Gene. (2011) 478:11–8. 10.1016/j.gene.2010.10.01121044875PravenecMSabaLMZídekVLandaVMlejnekPŠilhavýJ. Systems genetic analysis of brown adipose tissue function. Physiol Genomics. (2018) 50:52–66. 10.1152/physiolgenomics.00091.201729127223GaoJYSongXYWuHTangQWeiZYWangXY. Detection of rs665862918 (15-bp Indel) of the HIAT1 gene and its strong genetic effects on growth traits in goats. Animals. (2020) 10:358–66. 10.3390/ani1002035832102183AljanabiSMMartinezI. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. (1997) 25:4692–3. 10.1093/nar/25.22.46929358185LanXYPanCYChenHZhang CL LiJYZhaoMLeiCZ. An AluI PCR-RFLP detecting a silent allele at the goat POU1F1 locus and its association with production traits. Small Ruminant Res. (2007) 73:8–12. 10.1016/j.smallrumres.2006.10.009LiXJiangEZhangKZhangSJiangFSongE. Genetic variations within the bovine CRY2 gene are significantly associated with carcass traits. Animals. (2022) 12:1616. 10.3390/ani1213161635804515ZhangSKangZSunXCaoXPanCDangR. Novel lncRNA lncFAM200B: molecular characteristics and effects of genetic variants on promoter activity and cattle body measurement traits. Front Genet. (2019) 10:968. 10.3389/fgene.2019.0096834745234LiuJJLiangAXCampanileGPlastowGZhangCWangZ. Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo. J Dairy Sci. (2018) 101:433–44. 10.3168/jds.2017-1324629128211ThomsenBHornPPanitzFBendixenEPetersenAHHolmLE. missense mutation in the bovine SLC35A3 gene, encoding a UDP-N-acetylglucosamine transporter, causes complex vertebral malformation. Genome Res. (2006) 16:97–105. 10.1101/gr.369050616344554RezvannejadEAsadollahpour NanaeiHEsmailizadehA. Detection of candidate genes affecting milk production traits in sheep using whole-genome sequencing analysis. Vet Med Sci. (2022) 8:1197–204. 10.1002/vms3.73135014209ZhangHShenLYXuZCKramer LM YuJQZhangXYNaW. Haplotype-based genome-wide association studies for carcass and growth traits in chicken. Poult Sci. (2020) 99:2349–61. 10.1016/j.psj.2020.01.00932359570Al-MamunHAKwanPClarkSAFerdosiMHTellamRGondroC. Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight. Genet Selection Evol. (2015) 47:66. 10.1186/s12711-015-0142-426272623MacéTGonzález-GarcíaEPradelJParisotSCarrièreFDoulsS. Genetic analysis of robustness in meat sheep through body weight and body condition score changes over time. J Anim Sci. (2018) 96:4501–11. 10.1093/jas/sky31830085118ZhangBChangLLanXYAsifNGuanFLFuDK. Genome-wide definition of selective sweeps reveals molecular evidence of trait-driven domestication among elite goat (Capra species) breeds for the production of dairy, cashmere, and meat. GigaScience. (2018) 7:giy105. 10.1093/gigascience/giy10530165633GlanzmannBLombardDCarrJBardienS. Screening of two indel polymorphisms in the 5'UTR of the DJ-1 gene in South African Parkinson's disease patients. J Neural Transm. (2014) 121:135–8. 10.1007/s00702-013-1094-x24052110OnerYKeskinAUstunerHSoysalDKarakaşV. Genetic diversity of the 3' and 5' untranslated regions of the HSP70. 1 gene between native Turkish and Holstein Friesian cattle breeds South African. J Anim Sci. (2017) 47:424–39. 10.4314/sajas.v47i4.2CuiYYanHWangKXuHZhangXZhuH. Insertion/Deletion within the KDM6A gene is significantly associated with litter size in goat. Front Genet. (2018) 9:91. 10.3389/fgene.2018.0009129616081SerroteCMLReinigerLRSSilvaKBRabaiolliSMDSStefanelCM. Determining the polymorphism information content of a molecular marker. Gene. (2020) 726:144175–8. 10.1016/j.gene.2019.14417531726084ZhaoHDHeSWangSHZhuYJXuHWLuoRY. Two new insertion/deletion variants of the PITX2 gene and their effects on growth traits in sheep. Anim Biotechnol. (2018) 29:276–82. 10.1080/10495398.2017.137941529200321GuidoCPanzaSSantoroMAvenaPPannoMLPerrottaI. Estrogen receptor beta (ERβ) produces autophagy and necroptosis in human seminoma cell line through the binding of the Sp1 on the phosphatase and tensin homolog deleted from chromosome 10 (PTEN) promoter gene. Cell cycle. (2012) 11:2911–21. 10.4161/cc.2133622810004VellingiriBIyerMDevi SubramaniamMJayaramayyaKSiamaZGiridharanB. Understanding the role of the transcription factor Sp1 in ovarian cancer: from theory to practice. Int J Mol Sci. (2020) 21:1153. 10.3390/ijms2103115332050495ChenSLiHZhangJJiangSZhangMXuY. Identification of Sp1 as a transcription activator to regulate fibroblast growth factor 21 gene expression. BioMed Res Int. (2017) 29:8402035. 10.1155/2017/840203528466020GuMCosenzaGIannacconeMMacciottaNPPGuoYDi StasioL. The single nucleotide polymorphism g. 133A>C in the stearoyl CoA desaturase gene (SCD) promoter affects gene expression and quali-quantitative properties of river buffalo milk. J Dairy Sci. (2019) 102:442–51. 10.3168/jds.2018-1505930391181LiKLiuYHeXTaoLJiangYLanR. novel SNP in the promoter region of IGF1 associated with Yunshang Black Goat kidding number via promoting transcription activity by SP1. Front Cell Dev Biol. (2022) 10:873095. 10.3389/fcell.2022.87309535646903