Seasonal and sexual variation in mRNA expression of selected adipokine genes affecting fat deposition and metabolism of the emu (Dromaius novaehollandiae)

Emus are farmed for fat production. Oil rendered from their back and abdominal fat pads has good anti-oxidant and anti-inflammatory properties and has ingredients that promote cell growth. Our objective is to examine the mRNA expression of 7 emu adipokine genes (eFABP4, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2, eLEP and eLepR) to identify gene markers that may help improve emu fat production. Back and abdominal fat tissues from 11 adult emus were biopsied at four time points (April, June, August and November). Total RNA was isolated and cDNA was synthesized. Gene specific primers were designed for partial cloning fragments to amplify the open reading frame of the 7 genes. eLEP was not expressed in emu fat tissue. Nucleotides and amino acids sequences of the 6 expressed gene were compared with homologs from other species and phylogenetic relationships established. Seasonal mRNA expression of each gene was assessed by quantitative RT-PCR and differential expression analysed by the 2−ΔΔCT method. The 6 expressed genes showed seasonal variation in expression and showed association of expression level with back fat adiposity. More whole-genome scanning studies are needed to develop novel molecular markers that can be applied to improve fat production in emus.


Scientific Reports
| (2022) 12:6325 | https://doi.org/10.1038/s41598-022-10232-w www.nature.com/scientificreports/ that regulates the synthesis of unsaturated fatty acids 15 ; (3) adiponectin (AdipoQ), an adipokine hormone that is mainly secreted from mammalian adipose tissue, is involved with lipogenesis and insulin resistance 16 . In birds, adiponectin receptors are expressed in a diversity of tissues and its function may be altered from that of mammals [17][18][19] ; (4) leptin, an adipocyte-derived hormone that regulates feeding behavior in mammals where energy expenditure via its interaction with the leptin receptor (LepR) belongs to the class I cytokine receptor superfamily 20 . The FABPs are abundant intracellular proteins that play important roles in the transportation and metabolism of long-chain fatty acids 14,21,22 . FABP family proteins could be used as tissue specific injury markers because they have high tissue specificity, abundance in the tissue, and low molecular weight (approx. 15 kDa) 21,23 . The development and growth of adipose tissue are due to increases of both adipocyte cell size and cell number. The FABP4 has been extensively used as a marker for differentiated adipocytes 14,24 .
AdipoQ is an adipokine hormone that influences several metabolic functions including glucose utilization, lipogenesis, energy homeostasis and immunity by signaling through two distinct receptors, AdipoR1 and AdipoR2. AdipoR1 is abundantly expressed in skeletal muscle, whereas AdipoR2 is predominantly expressed in the liver. AdipoR1 and AdipoR2 genes are ubiquitously expressed in chicken tissues and their expression is altered by feed deprivation in the anterior pituitary gland and adipose tissue 25 .
Leptin (LEP) is an adipokine hormone that is the central mediator in a negative feedback loop that regulates energy homeostasis through the hypothalamus. In mammals, LEP administration leads to reduced food intake, increased energy expenditure and weight loss 26 . Unlike in mammals, LEP shows no expression in adipose tissue of the few avian species examined so far 27,28 , while its receptor (LepR) is still weakly expressed with no correlation to adiposity in chickens 29 . Previous studies revealed that avian LepR shares signal transduction pathway via administration of mammalian LEP 30 .
In this study, we first cloned these adipokine genes, FABP4, SCD1, AdipoQ, AdipoR1, AdipoR2, and LepR, from emu adipose tissue and investigated the seasonal gene expression profile associated fat deposition with the intention to identify genetic markers for improving fat production in the emu.
We tested the following null hypotheses:

Methods
The study was conducted in accordance with the relevant guidelines and regulations. Methods are reported in accordance with ARRIVE guidelines (https:// arriv eguid elines. org).
Animal tissue. We tracked 11 adult emus (7 males, 4 females) over one breeding season (TryHarder farm, Saskatchewan, Canada). Birds were weighed and back fat and abdominal fat tissues were biopsied with a tissue punch with plunger (diameter core size: 6.0 mm, TED PELLA Ltd.) at four time points (April, June, August and November 2011). These time points were chosen as April was the time when birds started to put on fat, June was the time when birds were gaining fat, August was the time when maximum fat was put on and between August and November there were minimum fat gain, November was the start of the breeding season and birds started to use back fat as an energy source. In addition, for June, August and November, 6 birds from the same flock were also sampled. Different 6 birds were sampled for each of these 3 time points (See Supplemental  _F1: 5′-GCC TGA CAA AAT GTG CGA C-3′  JN663389 (399 bp)   FABP_R1: 5′-AAG AGT TTA CGA AAG AGC ATG AGG AA-3′   eSCD1   SCD_F1: 5′-CAC ATG CCT GCG CAC TTG CTACA-3′  JN663390 (1083 bp)   SCD_R1: 5′-GAC TAC TCC ACC AGT GAG TTT GGC TGGC-3′   SCD_F2: 5′-GGA ATA TCA TCC TCA TGA GCC TGC TGCA-3′   SCD_R2: 5′-TGG GAG TCA CAA GAG CGG CTG AGT TC-3′   eAdipoQ   AdipoQ_F1: 5′-ACG TTT ACC GCT CCG CCT TCA GCG T - Table 2. Primers used for amplifying eLept from emu adipose tissue*. *Because in the few avian species that have been examined, LEP expression was not found in adipose tissue, we have designed primers based on the conserved region of LEP from 6 different species (wild mallard, Japanese eel, catfish, mouse, chicken and turkey) to determine whether LEP expression can be found in emu adipose tissue.

Fatty acids analysis.
Fat samples collected in November were thawed, weighed and ground (cold pressed).
The samples were then placed in a double boiler pan and heated at 70 °C for 30 min to yield the rendered oil. The resulting oil samples were then filtered, weighed and stored under nitrogen in white plastic bottles at 4 °C. Total lipids in the oil samples were extracted with chloroform:methanol (2:1 v/v) containing 0.01% butylated hydroxytoluene (BHT) as the antioxidant. Fatty acid methyl esters (FAME) were prepared by transesterification with Boron trifluoride (BF3) in methanol following the method described by Kitts et al. 34 , based on the procedure described by Ackman et al. 35 . After the extraction with hexane, FAMEs were analyzed by gas liquid chromatograph (GC-17A, Shimadzu Scientific Instruments Inc., Columbia, Maryland), equipped with flame ionization detector and an auto injector (AOC1400, Shimadzu Scientific Instruments Inc., Columbia, Maryland). Samples were injected onto a capillary column (30 m × 0.25 mm; 0.25 μm film thickness; liquid phase: J and W DB 23), with helium as the carrier gas. Temperature programming was used according to the method described by Budge et al. 36 with minor modification. The column temperature was initially set at 153 °C for 2 min, then increased to 174 °C at 2.3 °C/min, and then to a final temperature of 220 °C at 2 °C/min with a final hold time of 2 min. Detector and injector temperatures were both set at 250 °C. Chromatographic peaks were integrated and identified using the Shimadzu software package (version 7.2.1 SP1, https:// www. shima dzu. com/ an/ produ cts/ softw are-infor matics/ softw are-option/ index. html), which were compared to known standards supplied by Nu-Chek Prep (Elysian, MN). Individual fatty acids are reported as weight percent of total fatty acids using mass response factors relative to C18:0.

Statistical analysis. Repeated Measures Mixed Model Least Squares ANOVA tests and multiple regression
analyses were performed using JMP 8.0 (SAS Institute, Cary, North Carolina, USA). Tukey's HSD was used for mean separation and the level of significance was defined at P < 0.05.
Statistical models:

Seasonal variation in gene expression
where Y ijkl = the mRNA expression level of an adipokine gene, μ = the grand mean, Sex i = the effect of sex, i = male or female, e ij = the effect of replication and the sub-plot error term for testing the sex effect, j = 1, 2, 3…17, Season k = the effect of season, k = April, June, August, or November time point, (Season x Sex) ik = the interaction between sex and season, e ijk = model error term. Both Season and Sex and their interaction were treated as fixed terms. Replication is considered as a random effect nested in Sex.

Association between fat gain and level of gene expression
where Ŷ is the is the predicted value of the dependent variable (i.e., amount (kg) of fat gain during the period specified (April-June, June to August, or August to November)); β 0 is the value of Y when all of the independent variables equal to zero (i.e. the Y-intercept); β 1 -β6 are the coefficients of the independent variables listed in the equation (eFABP4, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2, eLepR which are the mRNA expression level of the adipokine genes at a month specified (April, June, or August).

Association between fatty acid profile and gene expression
where Ŷ is the predicted value of proportion of a Fatty Acid in the oil extracted from emu fat tissue collected in November; β 0 is the value of Y when all of the independent variables are equal to zero (i.e. the Y-intercept); β 1 -β 7 are the coefficients of the independent variables listed in the equation; Sex (female) is a one for female birds and zero for male birds; eFABP4, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2, and eLepR are the mRNA expression level of the adipokine genes in November.

Ethics approval.
All experiments were performed in accordance with protocols reviewed and approved by the UBC Animal Care Committee (Certificate # A10-0106).

Results
Isolation, amplification and cloning. The full-length cDNA of eFABP4, eSCD1, eAdipoQ, eAdpoiR1, eAdipoR2, eLepR, and eβ-actin were isolated and cloned from emu adipose tissue with gene specific primers conserved by other avian species EST database ( Table 3). Nomenclature of each gene was based on identities of the primary gene structure to other homologs (Table 1), and was assigned the following GenBank accession numbers: , and housekeeping gene eβ-actin (JN663391). Nucleic acid sequence of eFABP4 (92-94% similarity) and eAdipoR2 (92-94% similarity), showed the highest similarity with orthologues of other avian species (wild mallard, Greylag goose, Swan goose, zebra finch, chicken, turkey, and pheasant) ( Table 4). eAdipoR1 also showed high similarity (92-94%) with other avian species with the exception of Zebra finch (only 90% similar). eSCD1 and eAdipoQ were the next high with 88-91% and 80-85%, respectively. We were not able to amplify any eLEP mRNA in emu adipose tissue. eLepR showed the most divergent among the 6 genes examined with similarity ranging from 76 (zebra finch) to 88% (Greylag and Swan geese).
Amino acid sequence similarity. The primary protein structure assigned the nomenclature and Genbank accession numbers are: eFABP4 (AET74082), eSCD1 (AET74083), eAdipoQ (AFF19461), eAdopoR1 (AFF19462), eAdipoR2 (AFF19463), eLepR (AFF19464) and eβ-actin (AET74084). The eFABP4 encodes a protein of 132 amino acids and shares high similarities with other avian species (average 96.1%) ( Table 5). The eSCD1 encodes a protein of 360 amino acids and also shares high similarities with other avian species (average 94%). eAdipoQ encodes a protein of 245 amino acids and shows average 88.5% similarities with other avian species. eAdipoR1 and eAdipoR2 are the most conservative proteins we have examined in this study. eAdipoR1 encodes a protein of 352 amino acids, shares 100% similarity with Swan goose and an average 89.8% similarities with all other species examined. eAdipoR2 encodes a protein of 385 amino acids and shows an average of 92.3% similarities with all other species. The eLepR encodes 1,151 amino acids and shares 90% similarities with the waterfowl group, 86% similarities with chicken, turkey and pheasant, and 79% with zebra finch. (Table 5).
Phylogeny tree. Phylogeny trees of the eFABP1, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2 and eLepR in this study in association with other selected species were constructed to illustrate their genetic relatedness (Fig. 1). Generally speaking, clustering of amino acid sequences resulted in trees showing that emu was closer to zebra finch and migrating waterfowl (swan goose, graylag goose, and mallard) than domestic and gallinaceous birds (turkey, chicken and pheasant). The eFABP4 and eSCD1 are less divergent than eAdipoQ. eAdipoR1 and eAdi-poR2 are both less divergent than eAdipoQ, while eAdipoR1 is the least divergent of the emu proteins examined. On the other hand, eLepR is the most divergent from the non-avian proteins.
Seasonal variation in emu back and abdominal fat weight gain. We recorded the body weights of the birds in April, June, August and November. After the birds were slaughtered in November, the back and abdominal fat pads were dissected out and weighed. We also weighed the bird carcasses without the fat pads. The carcass weight (34.8 ± 1.2 kg) was not significantly (P < 0.91) different from the April body weight (34.6 ± 1.2 kg). We therefore concluded that the difference in body weight between two time points would be a good estimate of fat gained between the two time points. In this study (Fig. 2), emus gained fat from April to August. From August to November (beginning of breeding season), there was very little fat gain. The mean fat gain from April to November was 11.3 ± 3.28 kg. There was no significant difference between males and females in fat gain.   (Table 6 and Fig. 3A). Sex effect and the Season X Sex interaction were not significant.

Regression of fat gained (kg) on eFABP4 expression levels. The amount of fat gained between April
and June regressed significantly and positively on eFABP4 expression level in April (Table 8 and Supplemental  Fig. 1A). The amount of fat gained between June and August regressed significantly and positively on eFABP4 expression level in June. However, the amount of fat gained between August and November regressed significantly but negatively on eFABP4 expression level in August.

Seasonal variation in eSCD1 expression.
There was a significant difference between males and females in their seasonal variation (Season X Sex interaction) in eSCD1 expression level (Table 7 and Fig. 3A). Male eSCD1 expression level was slightly higher (not statistically significant) in April and remained low from June to November. Female eSCD1 level was significantly higher in August but not significantly different than male levels in April, June and November.
Regression of fat gained (kg) on eSCD1 expression levels. Fat gain was not associated with eSCD1 expression level except from June to August. The amount of fat gained during that period significantly regressed on June eSCD1 expression level (Table 8 and Supplemental Fig. 1B).
Seasonal variation in eAdipoQ expression. eAdipoQ expression level was not significantly (P = 0.23) affected by season and there was no significant (P = 0.10) difference between males and females ( Table 6 and Fig. 3A).

Regression of fat gained (kg) on eAdipoQ expression level.
Fat gain from April to June was not associated with eAdipoQ expression level. Fat gain from June to August significantly regressed on June eAdipoQ level (Table 8 and Supplemental Fig. 1C). Fat gain from August to November significantly but negatively on August eAdipoQ level.
Seasonal variation in eAdipoR1 expression. eAdipoR1 expression was significantly higher in April than the rest of the year (Table 6 and Fig. 3B). Sex effect and the Season X Sex interaction were not significant.

Regression of fat gained (kg) on eAdipoR1 expression level. The amount of fat gained from April
to June significantly but negatively regressed on April eAdipoR1 level (Table 7 and Supplemental Fig. 1D). From June to August, fat gain regressed significantly on June eAdipoR1 level.
Seasonal variation in eAdipoR2 expression. eAdipoR2 expression level was highest in April, intermediate in June and August, and lowest in November (Table 6 and Fig. 3B). Sex effect and the Season X Sex interaction were not significant.     (Table 7 and Supplemental Fig. 1E).

Seasonal variation in eLepR expression. eLepR expression level was highest in April, intermediate in
June and August, and lowest in November (Table 6 and Fig. 3B). Sex effect and the Season X Sex interaction were not significant.  (Table 7 and Supplemental Fig. 1F).

Discussion
In order to better understand the genetics of adiposity in emu, a bird that has huge seasonal variation in fat deposition and metabolism, we have selected 7 adipokine genes for examination that have been established to be involved with fat deposition and metabolism. We were able to reject all the null hypotheses that we set out to test. Of the 6 genes (eFABP4, eSCD1, eAdipoQ, eAdipoR1, eAdipoR2, and eLepR) expressed, they shared high nucleotides sequences and amino acids sequences similarity with other avian species (H 0 1 rejected). Except for eAdipoQ, there were seasonal variation in the mRNA expression between April and November (H 0 2 rejected), and the expression level of all 6 genes were associated with back fat weight gain during some time in the season (H 0 3 rejected). eAdipoR1 expression level was associated with the % of Palmitic acid in emu oil while eAdipoR2 level was associated with % Lenoleic acid and % α-Linolenic Acid in emu oil (H 0 5 rejected). There were significant differences between males and females in % Palmitic acid and % Oleic Acid in emu oil (H 0 6 rejected).  www.nature.com/scientificreports/ We were not able to detect any eLEP mRNA expression in the emu adipose tissue. Based on the lack of LEP expression in the adipose tissue of zebra finch 37 , jungle fowl 38 , several lines of commercial chickens 39 , rock dove 40 , and quail 29 , Friedman and Seroussi 28 suggested that LEP is not expressed in avian adipose tissue. More avian species have to be examined to confirm this observation. Our result from emu, a ratite that is phylogenetically distant from the birds examined so far, provided support to their observations (H 0 4 rejected). In birds, LEP is expressed in brain tissue, adrenal glands and gonads, but is not expressed in the liver and is generally not detectable in the blood. LepR receptors are predominantly expressed in the pituitary. Seroussi et al. 29 reported that in chicken, ducks, and quail adipose tissue, LEP and LepR were scarcely transcribed, and the expression level was not correlated to adiposity. They proposed that LEP in birds may act as an autocrine or paracrine instead of being a circulating hormone as in mammals. These observations, mostly from chicken studies, allowed Friedman and Seroussi 28 to speculate that avian adipose tissue does not control appetite, insulin resistance, or inflammation.       www.nature.com/scientificreports/ We have detected low expression levels of eLepR in emu fat tissue. Expression level was highest in April and stepwise decreased to the lowest level in November, which was opposite to the fat weight gain trend. The amount of back and retroperitoneal fat gain between June and August regressed significantly but negatively on eLepR June expression level. Since there was no eLEP expression in emu fat tissue, it seems likely that LEP in emu is still a circulating hormone that affects fat deposition and metabolism. Our phylogenetic analysis found that emu LepR is closer to those of migrating waterfowl than other bird species examined. In mammals, LEP specially repressed the expression of SCD1 and reduced the accumulation of hepatic triglycerides, cholesterol esters and VLDL synthesis 26,41,42 . It is suspected that the role of LEP in governing adipose tissue regulation of appetite and energy expenditure has been altered in birds 28 . Never the less, the relationship between LEP and the loss of appetite over the winter breeding period in emus remains to be studied. eSCD1 is expressed in emu adipose tissue. There was little seasonal variation in expression except that fat samples collected from females in August had a significant 35-fold increase in expression. Fat weight gained between June and August regressed significantly on June eSCD1 expression level. Individuals that had high fat weight gain would have high eSCD1 expression but low eLepR expression and vice versa for individuals that had little fat weight gain. This would indicate that emu leptin can suppress the expression of eSCD1 as seen in mammals. SCD1 is also transcriptionally regulated by a number of factors in mammals, including sterol regulatory element-binding protein-1 (SREBP-1) and polyunsaturated fatty acids 43,44 .
SCD1 is predominately located in the endoplasmic reticulum and catalyzes the rate-limiting step in the cellular synthesis of mono-unsaturated fatty acids from saturated fatty acids 15,45 . SCD-1 converts the saturated fatty acids, palmitic acid (16:0) and stearic acid (18:0), to generate the mono-unsaturated fatty acids, palmitoleic (16:1 n7) and oleic acid (18:1 n9), which are accumulated as triglycerides in adipose tissues [45][46][47][48] . Oleic acid is the predominant fatty acid in emu adipose tissue. A proper ratio of saturated fatty acids to mono-unsaturated fatty acids contributes to membrane fluidity. In mice, SCD-1, known as a lipid synthesis enzyme, also plays a role in upregulating lipid mobilization through its desaturation product, oleic acid 49 . Specific unsaturated fatty acids are preferentially used during metabolism over saturated fatty acids [50][51][52] .
Fat storage and usage in birds are mainly for survival, migration and reproductive performance [53][54][55][56] . Catbirds increased adipose storage during spring and autumn migration, showing increased rates of basal lipolysis during migration and tropical overwintering 57 . In our study, emus started gaining fat in April and the rate of gain was maximized between June and August. Fat weight gain between June and August significantly regressed on June eSCD1 expression. From August to November fat gain was minimal. It was during this period when female emus were getting ready to lay eggs. There was a 35-fold increase in eSCD1 expression in females in August. They may be optimizing the fatty acid composition of the adipose tissue to get ready for the mobilization of lipids into the ovary for formation of the egg yolk. Unfortunately, we only had the fatty acid profile of emu oil collected in November and by that time eSCD1 expression was extremely low and we found no association of eSCD1 expression and fatty acid profile of emu oil.
FABPs are a family of proteins known as intracellular lipid chaperones that regulate lipid trafficking and responses in cells 58 . FABP genes have been shown to be associated with lipid metabolism (lipogenesis and lipolysis), homeostasis in adipocytes, marbling and back fat deposition [59][60][61] . FABP4 is highly expressed in adipocytes and its expression can be highly induced during adipocyte differentiation which is transcriptionally controlled by peroxisome proliferator-activated receptor (PPAR) γ agonists, fatty acids, dexamethasone and insulin 57,62 . It has also been postulated that FABP4 can activate hormone sensitive lipase (HSL) in adipocytes to regulate lipolysis 63,64 . In chickens, earlier studies that examined the relationship of FABP4 with growth and fat accumulation reported results ranging from no association with fat accumulation in hybrid chickens 65 , significant positive correlation with abdominal fat in Luyuan chickens 66 , to correlation with growth depression in Arbor Acre genotype but strong positive association with growth performance of Cobb genotype 61 . In our study, eFABP4 expression in emu adipose tissue was high both in April and November and relatively low in June and August. Fat gain from April to August regressed positively on April and June eFABP-4 expression, respectively. However, fat gain from August to November regressed negatively on August eFABP4 expression. From August to November, fat gain was minimal and a couple of birds even had negative fat gain. By this time of the year, emus started to draw on the energy from the accumulated fat and the role of FABP4 switched from lipogenesis to lipolysis 57 . eFABP4 expression was highest in November and this may be an indication that the birds were more and more dependent on fat for energy because they have very little feed intake during breeding. In geese, FABP4 was found to be involved in lipid transportation and metabolic process, follicle development and final egg production. FABP4 was upregulated in the laying group compared with the pre-laying group 67 .
AdipoQ has been originally identified as a protein secreted and expressed exclusively in adipose tissue [68][69][70] . AdipoQ showed many functions like expanding fatty acids oxidation, controlling glucose level and managing receptor activity. In humans, Adiponectin is known to stimulate the expression of FABP 16 . In chickens, adiponectin plays important roles in energy homeostasis, body weight, lipid metabolism, and insulin sensitivity [71][72][73] . In broiler chickens, Tahmoorespur et al. 74 showed that AdipoQ expression in adipose tissue was inversely related to chicken abdominal fat deposition levels. Adiponectin has an effect on the impairment of adipocyte differentiation, which contributes to the negative regulation of fat deposition in chicken 71 . In adipose tissue of adult chickens, AdipoQ expression is higher in females than males, but AdipoR1 expression was higher in males than females 17 . In female birds, Adiponectin is secreted into the blood from adipocytes with a higher serum level 75 . In emus, eAdipoQ expression did not show any seasonal (April to November) nor sexual variation. Emu fat gain from June to August regressed positively on June eAdipoQ expression but fat gain from August to November regressed negatively on August eAdipoQ expression. Similarly, fat gain from June to August also regressed positively on June eAdipoR1 expression while fat gain from August to November regressed negatively on August eAdipoR2 expression. White-throated sparrows increase fat deposits during pre-migratory periods and rely on these fat stores to fuel migration. In the adipose tissue, there was a significant change in the biological control of adipokine www.nature.com/scientificreports/ expression from pre-migratory conditions to migratory conditions. It was proposed that Adiponectin may play a role in the switch from fat deposition to lipid metabolism as the main source of energy to fuel migratory flight in birds 76 . In emus, eAdipoR1/R2 expression was highest in April, before the birds started gaining fat. eAdipoR1/ R2 expression declined until the lowest level in November. Interestingly, in the oil extracted from emu fat in November, the % of Palmitic Acid (C16:0) regressed negatively on November eAdipoR1 expression. Both Lenoleic Acid (C18:2n6) and α-linolenic acid (C18:3n3) regressed negatively on November eAdipoR2 expression. These observations support the role of Adiponectin in lipid metabolism for converting the stored fat to energy during the period of low feed intake 50,51 . In emus, back fat showed a higher level of protein, cholesterol, C16:1 and the elements K, P, Si, Na, Ca, Mg, Fe, Zn, Se and Cu than abdominal fat. Abdominal fat was characterized by higher content of fat and ash, as well as Mn and Ba. Regardless of back or abdominal fat, there was generally high content of MUFA and PUFA. Males have higher content of Si, Ca, Cu, Sr in the adipose tissue than female 77 . In chickens, the most promising candidate genes affecting polyunsaturated fatty acids percentage were FADS2, DCN, FRZB, OGN, PRKAG3, LHFP, CHCHD10, CYTL1, FBLN5, and ADGRD1 78 . There are two major methods of quantitative trait loci (QTL) determination, the candidate gene approach and the whole-genome scanning. The candidate gene approach is used to detect QTL (Quantitative Trait Loci) responsible for genetic variation in the traits of interest. In chickens, FABP4 gene polymorphism has been associated with abdominal fat weight and percentage of abdominal fat, and FABP4 gene could be a candidate locus or linked to a major gene(s) that affects abdominal fat content 79 . In pigs, SCD-1 expression plays a critical role in adipocyte differentiation and has been identified as the promising candidate gene for less back fat deposition 80,81 . Fat deposition in emu is seasonal 77 . Under natural conditions, increased fat deposition in birds is for energy storage to cope with migration or periods when food is scarce, and is not associated with obesity. In addition to being an energy storage organ, the adipose is also an endocrine organ influencing reproduction, feeding behaviour, insulin sensitivity and disease resistance 75 . Whether genetic selection for increased fat deposition would lead to complications with obesity must be considered. In broiler chickens and turkeys, the selection for fast growth rate inevitably led to increase abdominal fat deposition and drastic reduction in breeder fertility 82,83 . Additionally, one has to consider whether selection for increased fat deposition would alter the fatty acids composition and other bioactive ingredients in the adipose tissue and thus affecting the efficacy of the emu oil 5 .

Conclusion
There has not been any reported genetic selection studies for increasing subcutaneous fat deposition in farm animals and no genetic improvement of emu productive traits has been carried out due to their short history of domestication compared to other livestock species 84 . Our study has laid down the groundwork for identifying promising candidate genes for such purpose. A follow-up whole-genome scanning study 85 explored the gene networks in emu adipose tissue affecting fat deposition and utilization and identified marker genes that deserve further analyses to develop novel molecular markers that can be applied to improve fat production in emus 86,87 .