Next Article in Journal
Context-Dependent Gestural Laterality: A Multifactorial Analysis in Captive Red-Capped Mangabeys
Next Article in Special Issue
Identification of Candidate Genes for Pigmentation in Camels Using Genotyping-by-Sequencing
Previous Article in Journal
Influence of Broiler Age on the Apparent Metabolizable Energy of Cereal Grains Determined Using the Substitution Method
Previous Article in Special Issue
First Report of Antimicrobial Susceptibility and Virulence Gene Characterization Associated with Staphylococcus aureus Carriage in Healthy Camels from Tunisia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries

by
Morteza Bitaraf Sani
1,*,
Zahra Roudbari
2,
Omid Karimi
3,
Mohammad Hossein Banabazi
4,5,
Saeid Esmaeilkhanian
4,
Nader Asadzadeh
4,
Javad Zare Harofte
1,
Ali Shafei Naderi
1 and
Pamela Anna Burger
6
1
Animal Science Research Department, Yazd Agricultural and Natural Resources Research and Education Center, Agricultural Research, Education & Extension Organization (AREEO), Yazd 8915813155, Iran
2
Department of Animal Science, Faculty of Agriculture, University of Jiroft, Jiroft 7867155311, Iran
3
Department of Animal Viral Diseases Research, Razi Vaccine and Serum Research Institute, Agricultural Research, Education and Extension Organization, Karaj 3146618361, Iran
4
Animal Science Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj 3146618361, Iran
5
Department of Animal Breeding and Genetics (HGEN), Centre for Veterinary Medicine and Animal Science (VHC), Swedish University of Agricultural Sciences (SLU), 75007 Uppsala, Sweden
6
Research Institute of Wildlife Ecology, Vetmeduni Vienna, 1160 Vienna, Austria
*
Author to whom correspondence should be addressed.
Animals 2022, 12(2), 184; https://doi.org/10.3390/ani12020184
Submission received: 15 December 2021 / Revised: 9 January 2022 / Accepted: 11 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Trends in Camel Health and Production)

Abstract

:

Simple Summary

This project aimed to find biological themes affecting growth in dromedaries. Candidate SNPs associated with growth were mapped to 22 genes, and 25 significant themes were identified related to growth. The main biological functions included calcium ion binding, protein binding, DNA-binding transcription factor activity, protein kinase activity, tropomyosin binding, myosin complex, actin-binding, ATP binding, receptor signaling pathway via JAK-STAT, and cytokine activity. EFCAB5, MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes are candidates for improving growth in camel breeding programs.

Abstract

Growth is an important heritable economic trait for dromedaries and necessary for planning a successful breeding program. Until now, genome-wide association studies (GWAS) and QTL-mapping have identified significant single nucleotide polymorphisms (SNPs) associated with growth in domestic animals, but in dromedaries, the number of studies is very low. This project aimed to find biological themes affecting growth in dromedaries. In the first step, 99 candidate SNPs were chosen from a previously established set of SNPs associated with body weight, gain, and birth weight in Iranian dromedaries. Next, 0.5 kb upstream and downstream of each candidate SNP were selected from NCBI (assembly accession: GCA_000803125.3). The annotation of fragments with candidate SNPs regarding the reference genome was retrieved using the Blast2GO tool. Candidate SNPs associated with growth were mapped to 22 genes, and 25 significant biological themes were identified to be related to growth in dromedaries. The main biological functions included calcium ion binding, protein binding, DNA-binding transcription factor activity, protein kinase activity, tropomyosin binding, myosin complex, actin-binding, ATP binding, receptor signaling pathway via JAK-STAT, and cytokine activity. EFCAB5, MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes are candidates for improving growth in camel breeding programs.

1. Introduction

Understanding genetic factors that influence animal growth is very important for accurate selection and genetic improvement of these quantitative traits. Growth traits, especially birth and weaning weight and gain per day, are economically important heritable traits and necessary for planning a successful camel breeding program [1,2]. Vice versa, good knowledge about genomic pathways underlying a breeding program is necessary to optimize the growth of animals [3]. Many SNPs associated with growth have been detected using GWAS and QTL mapping in livestock [4]; however, few projects have been performed yet for the genetic improvement of growth and meat production in camels.
Approximately 3,000,000 camels are slaughtered annually, and camel meat production reaches 653,000 tons (FAO 2021), although contributing approximately only 0.18% to total red meat production. There are approximately 140,000 dromedaries in Iran (FAO, 2019) that are four basic types, milk, meat, dual purpose, and riding camels [5]. The majority of dromedaries in Iran’s central desert belong to the meat type, which is characterized by a heavy and large head, short neck, wide posterior pars, and firm body [6]. Most dromedary herds are kept semi-extensive or extensive and are located around the desert or near a village. The owners usually gather calves in spring for sale and fattening. The herd sizes range from 4–5 to 100–150 heads [1].
In recent GWAS using genotyping-by-sequencing, Bitaraf Sani et al. [1] identified 99 SNPs potentially associated with important growth traits to improve camel breeding, namely birth weight, daily gain, and body weight. Taking advantage of the comprehensive dataset, here, we aimed to find biological themes and genes affecting growth in dromedaries using a previously established set of associated SNPs.

2. Materials and Methods

Biological themes that related to growth were identified using gene-set enrichment analysis. In a previous study, we collected data on 51 herds of dromedaries in five regions of the central desert of Iran in 2018. Among the registered 964 calving between January to May, we genotyped 96 calves and collected 252 body weight records in different times. The samples were genotyped-by-sequencing using two restriction enzymes combination, EcoR1 and HinF1, and paired-end (150 bp) sequencing (10 X) on the Illumina HiSeq 2000. Among 14,522 markers, a total of 99 SNPs were associated with growth traits, and twelve haplotype blocks and 80 tag SNPs were predicted. The accuracy of GEBVs based on the 99 associated SNPs was 0.62, 0.82, and 0.57 for birth weight, daily gain, and body weight. The 99 associated markers with body weight, gain, and birth weight in dromedaries were included as candidate SNPs [1]. Next, 0.5 kb upstream and downstream of each candidate SNP were selected from NCBI. The fragments were aligned to the dromedary reference genome (assembly accession: GCA_000803125.3) using Blast2GO OmicsBox (version 1.4.11), including all (non)coding and pseudo genes [7]. Only genes with at least one significant SNP (Blast Expectation Value = 0.001) were considered significant genes. We identified the functional classes of the genes based on five databases, including GO (http://www.geneontology.org; accessed on 10 January 2022), KEGG (http://www.genome.jp/kegg; accessed on 10 January 2022), Panther (http://www.pantherdb.org; accessed on 10 January 2022), Wikipathway (http://www.Wikipathway.org; accessed on 10 January 2022), and Reactome (http://www.reactome.org; accessed on 10 January 2022). Genes located in the same functional class were considered a group of genes with some special and common features, such as participation in three ontological processes, including biological processes, molecular function, and cellular components. Finally, significant functional classes related to growth traits were tested using hyper geometric distribution and Fisher’s exact test by the software package goseq [8].

3. Results and Discussion

This study aimed to understand biological themes underlying important growth-related traits in dromedaries. We identified 22 genes involved in 25 functional themes (Table 1). Calcium ion binding, including the EFCAB5 gene, was the most significant terms (e-value = 4.74 × 10−96) related to growth. The other main significant terms were DNA-binding transcription factor activity (e-value = 1.13 × 10−40), protein kinase activity, tropomyosin binding, myosin complex, actin-binding (e-value = 5.34 × 10−24), ATP binding (e-value = 5.34 × 10−24), receptor signaling pathway via JAK-STAT, and cytokine activity (e-value = 2.93 × 10−42), which are enriched by the MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes. The seven biological processes related to growth were detected, including protein phosphorylation; regulation of transcription DNA-templated, protein-end actin filament capping; defense response to the virus; intercellular signal transduction; and response to stimulus (Figure 1). There were four cellular components, including extracellular region, cytoplasm, myosin complex, and nucleus (Figure 2). The 11 significant molecular functions were calcium ion binding, protein binding, cytokine activity, actin-binding, protein kinase activity, ATP binding, guanyl-nucleotide exchange factor activity, DNA-binding transcription factor activity, tropomyosin binding, and motor activity (Figure 3). One of the important genes in these terms was MYO3A in the myosin gene superfamily in muscular-skeletal cells [9]. This gene is a protein that plays a major role in muscle regulation and as a transcription factor that directly regulates the expression and differentiation of muscle cells. It is also a regulatory factor in the myogenesis process. MYO3A gene has a functional relationship with body mass index [10]. The EFCAB5 gene is part of the calcium ion binding pathway and is associated with body height and bone density [11]. All muscle fibers use Ca2+ as their main regulatory and signaling molecule. IFNL3 was another gene related to growth in dromedaries in this study. The interferons (IFNs) are classic antiviral cytokines that create the first defense line against infection using antiviral, antitumor, and immunomodulatory activities [12]. The PREX1 (Phosphatidylinositol-3,4,5-Trisphosphate Dependent Rac Exchange Factor 1) gene involved in the guanyl-nucleotide exchange factor activity pathway encodes a protein that is bonded with guanosine diphosphate (GDP) to release guanosine triphosphate (GTP) and regulate the cytoskeletal structure, gene expression, and reactive oxygen species (ROS) production. Moreover, PREX1 regulates neutrophil responses [13] and plays a key role in actin cytoskeletal rearrangement, adhesion, and the production of reactive oxygen species (ROS) [14]. Additionally, the activity of PREX1 protein correlates with the induction of actin-mediated membrane ruffling. Tbx15, a T-box transcription factor, is to be involved in the regulation of transcription, DNA-templated, DNA binding, DNA-binding transcription factor activity, and nucleus pathways. Tbx15 has an important role in developing the skeleton, vertebral column, and head by controlling the number of mesenchymal precursor cells and chondrocytes [15]. Moreover, it affects the development of different tissues and organs [16]. TMOD2 affects growth trait through erythrocyte development, pointed-end actin filament capping, and tropomyosin binding pathways [17].

4. Conclusions

This report has identified some terms related to growth in dromedaries using gene-set enrichment analysis. The main significant biological activities were calcium ion binding, protein binding, DNA-binding transcription factor activity, protein kinase activity, tropomyosin binding, myosin complex, actin-binding, ATP binding, receptor signaling pathway via JAK-STAT, and cytokine activity, which can be helpful in improving the understanding of mechanisms of growth. EFCAB5, MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes are candidates for improving growth in dromedaries. These results provide some insights into camel breeding programs.

Author Contributions

Conceptualization, M.B.S., J.Z.H., S.E. and M.H.B.; methodology, M.B.S., J.Z.H., M.H.B., Z.R. and A.S.N.; software, M.B.S. and Z.R.; formal analysis, M.B.S.; investigation, M.B.S., N.A. and O.K.; writing—original draft preparation, M.B.S.; writing—review and editing, S.E., M.H.B. and P.A.B.; visualization, M.B.S.; supervision, M.H.B., S.E. and P.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

P.A.B. acknowledges funding from the Austrian Science Fund (FWF): P29623-B25. This research was jointly funded by Animal Science Research Institute of Iran (ASRI), Animal Breeding Center of Iran, and Yazd Agricultural and Natural Resources Research and Education Center, grant number 34-64-1357-005-970180.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during the current study are available in the [dryad] repository, https://doi.org/10.5061/dryad.n2z34tmwf; accessed on 10 January 2022.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bitaraf Sani, M.; Zareh Harofte, J.; Banabazi, M.H.; Esmaeilkhanian, S.; Shafie Naderi, A.; Salim, N.; Teimoori, A.; Bitaraf, A.; Zadehrahmani, M.; Burger, P.A.; et al. Genomic prediction for growth using a low-density SNP panel in dromedary camels. Sci. Rep. 2021, 11, 7675. [Google Scholar] [CrossRef] [PubMed]
  2. Wei, C.H.; Liu, C.S. Encyclopedia of Modern Sheep Production Technology 70–77; China Agriculture Press: Beijing, China, 2014. [Google Scholar]
  3. Alemneh, T.; Getabalew, M. Factors influencing the growth and 440 development of meat animals. Int. J. Anim. Sci. 2019, 3, 1048. [Google Scholar]
  4. Ghiasi, H.; Abdollahi-Arpanahi, R. The candidate genes and pathways affecting litter size in sheep. Small Rumin. Res. 2021, 205, 106546. [Google Scholar] [CrossRef]
  5. Salehi, M.; Gharahdaghi, A.A. Camel Production Potential and Recent Research in Iran. 2013. Available online: http://agris.fao.org/agris-search (accessed on 10 January 2022).
  6. Khodai, S.A. The report on camel production systems and the socio-economics of camel herders in the Islamic Republic of Iran. Cardn/Acsad/Camel 2004, 107, 4–10. [Google Scholar]
  7. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. goseq: Gene Ontology testing for RNA-seq datasets. R Bioconduct. 2012, 30, 1–25. [Google Scholar]
  9. Chen, C.P.; Chen, X.; Qiao, Y.N.; Wang, P.; He, W.Q.; Zhang, C.H.; Zhao, W.; Gao, Y.Q.; Chen, C.; Tao, T.; et al. In vivo roles for myosin phosphatase targeting subunit-1 phosphorylation sites T694 and T852 in bladder smooth muscle contraction. J. Physiol. 2015, 3, 681–700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Graff, M.; Scott, R.A.; Justice, A.E.; Young, K.L.; Feitosa, M.F.; Barata, L.; Winkler, T.W.; Chu, A.Y.; Mahajan, A.; Hadley, D.; et al. Genome-wide physical activity interactions in adiposity―A meta-analysis of 200,452 adults. PLoS Genet. 2017, 13, 1006528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Cheung, C.L.; Sham, P.C.; Xiao, S.M.; Bow, C.H.; Kung, A.C. Meta-analysis of gene-based genome-wide association studies of bone mineral density in Chinese and European subjects. Osteoporos. Int. 2012, 23, 131–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Møhlenberg, M.; Terczynska-Dyla, E.; Thomsen, K.L.; George, J.; Eslam, M.; Grønbæk, H.; Hartmann, R. The role of IFN in the development of NAFLD and NASH. Cytokine 2019, 124, 154519. [Google Scholar] [CrossRef] [PubMed]
  13. Welch, H.C.; Condliffe, A.M.; Milne, L.J.; Ferguson, G.J.; Hill, K.; Webb, L.M.; Okkenhaug, K.; Coadwell, W.J.; Andrews, S.R.; Thelen, M.; et al. P-Rex1 regulates neutrophil function. Curr. Biol. 2005, 15, 1867–1873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Jaffe, A.B.; Hall, A. Rho GTPases: Biochemistry and biology. Annu. Rev. Cell Dev. Biol. 2005, 21, 247–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Singh, M.K.; Petry, M.; Haenig, B.; Lescher, B.; Leitges, M.; Kispert, A. The T-box transcription factor Tbx15 is required for skeletal development. Mech. Dev. 2005, 122, 131–144. [Google Scholar] [CrossRef] [PubMed]
  16. Papaioannou, V.E. T-box genes in development: From hydra to humans. Int. Rev. Cytol. 2001, 1, 1–70. [Google Scholar]
  17. Yamashiro, S.; Gokhin, D.S.; Kimura, S.; Nowak, R.B.; Fowler, V.M. Tropomodulins: Pointed-end capping proteins that regulate actin filament architecture in diverse cell types. Cytoskeleton 2012, 69, 337–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Gene ontology of biology processes for growth traits in dromedaries.
Figure 1. Gene ontology of biology processes for growth traits in dromedaries.
Animals 12 00184 g001
Figure 2. Gene ontology of cellular components for growth traits in dromedaries.
Figure 2. Gene ontology of cellular components for growth traits in dromedaries.
Animals 12 00184 g002
Figure 3. Gene ontology of molecular function for growth traits in dromedaries.
Figure 3. Gene ontology of molecular function for growth traits in dromedaries.
Animals 12 00184 g003
Table 1. Gene set enrichment analysis significantly (p < 0.001) associated with growth traits in dromedaries.
Table 1. Gene set enrichment analysis significantly (p < 0.001) associated with growth traits in dromedaries.
SNP ListDescription/Gene Symbol e-Valuesim meanGO IDsGO Names
Chr16:34483240EFCAB54.74 × 109692.67F: GO:0005509F: calcium ion binding
Chr14:30865111endonuclease/reverse transcriptase5.76 × 105661.08F: GO:0005515F: protein binding
Chr14:30865085endonuclease/reverse transcriptase6.36 × 105661.1F: GO:0005515F: protein binding
Chr17:75949MTIF27.48 × 105689.27
Chr14:30865065endonuclease/reverse transcriptase1.52 × 105361.53F: GO:0005515F: protein binding
Chr19:10237661IFNL32.93 × 104292.11P: GO:0007259; P: GO:0050778; P: GO:0051607; F: GO:0005125; C: GO:0005576P: receptor signaling pathway via JAK-STAT;
P: positive regulation of immune response;
P: defense response to virus;
F: cytokine activity;
C: extracellular region
Chr9:22550930TBX151.13 × 104097.45P: GO:0006355; F: GO:0003677; F: GO:0003700; C: GO:0005634P: regulation of transcription, DNA-templated;
F: DNA binding;
F: DNA-binding transcription factor activity;
C: nucleus
Chr19:10237641IFNL14.46 × 103893.51P: GO:0007259; P: GO:0050778; P: GO:0051607; F: GO:0005125; C: GO:0005576P: receptor signaling pathway via JAK-STAT;
P: positive regulation of immune response;
P: defense response to virus; F: cytokine activity;
C: extracellular region
Chr19:9631630LORF21.19 × 103161.66
Chr6:14654249TMOD32.67 × 1027100P: GO:0048821; P: GO:0051694; F: GO:0005523P: erythrocyte development;
P: pointed-end actin filament capping;
F: tropomyosin binding
Chr35:9290500myo3a5.34 × 102494.62P: GO:0006468;P: protein phosphorylation;
P: visual perception;
P: response to stimulus;
F: motor activity;
F: actin binding;
F: protein kinase activity;
F: ATP binding;
C: cytoplasm;
C: myosin complex
Chr31:16977384NUTM2E5.71 × 1022100
Chr19:11148798PREX11.29 × 102188.94P: GO:0035556; F: GO:0005085P: intracellular signal transduction;
F: guanyl-nucleotide exchange factor activity
ChrX:113249264LINE-1 retro transposable element ORF2 protein3.24 × 101684.22
Chr33:4467956hypothetical protein CB1_0005680258.33 × 101591.67
ChrX:60452363hypothetical protein HJG63_0114878.67 × 101272.06
Chr12:5711reverse transcriptase family protein3.45 × 10870.44
Chr18:29830420DUF1725 domain-containing protein9.72 × 10767.49
Chr11:72060113KCNK185.06 × 10666.67
Chr1:63040834POLR2A1.03 × 10493.09
Chr1:63040824POLR2A1.08 × 10488.89
Chr11:72356436hypothetical protein HJG59_0078182.44 × 10450
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sani, M.B.; Roudbari, Z.; Karimi, O.; Banabazi, M.H.; Esmaeilkhanian, S.; Asadzadeh, N.; Zare Harofte, J.; Shafei Naderi, A.; Burger, P.A. Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries. Animals 2022, 12, 184. https://doi.org/10.3390/ani12020184

AMA Style

Sani MB, Roudbari Z, Karimi O, Banabazi MH, Esmaeilkhanian S, Asadzadeh N, Zare Harofte J, Shafei Naderi A, Burger PA. Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries. Animals. 2022; 12(2):184. https://doi.org/10.3390/ani12020184

Chicago/Turabian Style

Sani, Morteza Bitaraf, Zahra Roudbari, Omid Karimi, Mohammad Hossein Banabazi, Saeid Esmaeilkhanian, Nader Asadzadeh, Javad Zare Harofte, Ali Shafei Naderi, and Pamela Anna Burger. 2022. "Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries" Animals 12, no. 2: 184. https://doi.org/10.3390/ani12020184

APA Style

Sani, M. B., Roudbari, Z., Karimi, O., Banabazi, M. H., Esmaeilkhanian, S., Asadzadeh, N., Zare Harofte, J., Shafei Naderi, A., & Burger, P. A. (2022). Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries. Animals, 12(2), 184. https://doi.org/10.3390/ani12020184

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