Next Article in Journal
Differential Expression Profiles and Bioinformatics Analysis of tRNA-Derived Small RNAs in Muscle-Invasive Bladder Cancer in a Chinese Population
Previous Article in Journal
Pharmacogenomics: A Step forward Precision Medicine in Childhood Asthma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolutionary and Association Analysis of Buffalo FABP Family Genes Reveal Their Potential Role in Milk Performance

1
Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Huazhong Agricultural University, Ministry of Education, Wuhan 430070, China
2
National Center for International Research on Animal Genetics, Breeding and Reproduction (NCIRAGBR), Huazhong Agricultural University, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Genes 2022, 13(4), 600; https://doi.org/10.3390/genes13040600
Submission received: 15 February 2022 / Revised: 11 March 2022 / Accepted: 24 March 2022 / Published: 28 March 2022
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:
The fatty acid-binding protein (FABP) family gene encode a group of proteins that affect long-chain fatty acid (LCFAs) trafficking and play a crucial function in the regulation of milk fat synthesis. Nevertheless, little is known regarding the identification of members, theevolutionary background, and functional characteristics of FABP genes in buffalo. In this study, in silico analysis was performed to identify the members of FABPs in buffalo. The results revealed that a total of 17 FABP genes were identified. Based on their phylogenetic relationships, these sequences clustered into five groups with similar motif patterns and gene structures. According to positive selection analyses, all duplicated gene pairs containing FABPs in buffalo had Ka/Ks (nonsynonymous/synonymous) ratios that were less than 1, suggesting that they were under purifying selection. Association analysis showed that one SNP in LOC102401361 was found significantly associated with buffalo milk yield. The expression levels of several FABPs in buffalo mammary epithelial cells were regulated by palmitic and stearic acid treatment. The findings of this study provide valuable information for further research on the role of FABPs in regulating buffalo milk synthesis.

Graphical Abstract

1. Introduction

The FABP family is a group of intracellular carriers of bioactive lipids that affect the trafficking of fatty acids [1,2]. LCFAs have been found in different substructures of cells, including the mitochondria, peroxisome, endoplasmic reticulum, lipid droplets, and nucleus, assisting in fat synthesis, assembly, and storage [3]. There are several small proteins in the FABP family of approximately 15 kDa which transport fatty acids together with intracellular retinol and retinoic acid-binding protein [4]. They belong to the intracellular lipid-binding protein (iLBP) superfamily, which includes subfamilies for cellular retinol-binding proteins (CRBPs) and cellular retinoic acid-binding proteins (CRABPs) [5]. Several multiple genes have appeared in this subfamily, giving rise to many iLBPs including 12 FABPs in vertebrates [6,7,8]. However, some members of this subfamily disappeared in some species, such as teleost fishes, FABP10, and FABP11, which are only found in nonmammalian vertebrates [9]. A previous study reported that the FABP family consisted of four exons and some are distributed on a single chromosome in mammals such as mice, rats, and humans [10,11]. Structural analysis indicated that this subfamily contained a three-element fingerprint domain shared by three motifs termed FATTYACIDBP 1–3 [12]. These three motifs have a similar structure and provide an internal hollow space that assists as a binding site for hydrophobic ligands [13]. According to phylogenetic analysis, FABPs can be divided into three clusters. Cluster 1 contains FABP1 and FABP6, and cluster 2 includes FABP3, FABP4, FABP5, FABP7, FABP8, FABP9 [13], and FABP12 [1], only FABP2 belongs to cluster 3.
Since the FABP family genes facilitate fatty acids binding and transfection, they also synthesize fatty acids and milk fat. FABP3, FABP4, and FABP5 were up-regulated during lactation and involved fatty acids trafficking towards milk TAG [14]. It was found that FABP3 provided stearoyl-CoA (or other substrates such as 16:0 and 18:1) [15] to SCD, which then released oleic acid to FABP4, which then released fatty acids to other enzymes involved in TAG synthesis. Liang et al. [16] found that FABP3 can regulate milk fat synthesis by modulating the expression of SREBP1 and PPARγ in cattle mammary epithelial cells. The composition of fatty acids in bovine milk is influenced by polymorphisms in FABP3 and FABP4 [17], and genetic polymorphisms of FABP3 are associated with backfat thickness in Korean native cattle [18]. The FABP4 gene was also linked to fat depth in beef cattle [19] and fat deposition in Australian cattle longissimus [20]. Moreover, it has been reported that the inhibition of FABP4 or knockout of it reduces osteoarthritis in mice induced by high-fat diets [21], and it was also found to be a therapeutic target for treating obesity-related cancers [22], indicating that this protein may regulate fat metabolism. The expression levels of FABP1 and FABP6 were also found upregulated after parturition, promoting fatty acids uptake and intracellular transport [23]. Evidence has demonstrated that different livestock species have different fat percentages in their meat and milk due to FABPs. However, despite buffalo being the second-largest milk producer worldwide, little is known about FABPs in buffalo.
As proteins bind and traffic with fatty acids, the FABP family genes are also, at the same time, affected by fatty acids. For example, it was reported that conjugated linoleic acid persuaded K-FABP and PPAR-delta expression in the skin of mice. Retinoic acid induces the expression of the PA-FABP (psoriasis-associated fatty acid-binding protein) gene in human skin. H-FABP expression was also up-regulated by several LCFAs treatments in cultured rat muscle cells [24] and other physiological experiments [25,26]. Long-chain polyunsaturated fatty acids stimulate cellular fatty acid uptake by regulating the expression of FABP3 and other enzymes in human placental choriocarcinoma cells [27] and in bovine mammary epithelial cells [28].
Buffalo milk is renowned for its high milk fat content and quality [29]. However, as important genes related to milk fat, the buffalo FABP family genes have not been completely identified and analyzed. The completion of the buffalo genome sequence has made it possible to perform genome-wide identification and phylogenetic analysis [30]. Therefore, the current study was designed to identify the FABP family genes in the buffalo genome and analyze their classification, protein motifs, and feature structures. Furthermore, to recognize crucial markers affecting buffalo milk, we performed an association analysis of the FABP family genes with buffalo milk production traits. Next, the expression profile of the subfamily was measured after adding three kinds of LCFAs into buffalo mammary epithelial cells and hepatocytes. Our research provides future functional characterization and influence on milk production of the FABP family genes in buffalo, which is essential for buffalo dairy herd improvement.

2. Methods and Materials

2.1. Genome-Wide Identification of FABP Genes

Data resources of the genome, proteome, and annotation of river buffalo (Bubalus bubalis, assembly UOA_WB_1) [31] were downloaded from the Genome database of NCBI (https://www.ncbi.nlm.nih.gov/, accessed on 1 September 2021). Both the hidden Markov model (HMM) and Basic Local Alignment Search Tool (BLAST) were used to identify all possible FABPs in the present research [32]. One hundred and nine protein sequences of buffalo (Bubalus ubalis, 10), cattle (Bos Taurus, 40), goat (Ovis aries, 20), horse (Capra hircus, 18), and sheep (Equus caballus, 21) were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov/, accessed on 1 September 2021) (Supplemental File S1). Using these sequences as queries, potential candidate FABPs were searched via BLASTP with a threshold of e-value = 10−6. Moreover, HMM profiles were constructed for detecting FABP homologous sequences in HMMER 3.2 [33] (http://hmmer.org/, ccessed on 1 September 2021) with the default setting. The sequences obtained from both approaches were considered candidate FABP family genes. Then, the sequences were submitted to ExPASy (https://web.expasy.org/compute_pi/, accessed on 1 September 2021) to calculate their molecular weight and isoelectric point.

2.2. Phylogenetic Analysis of FABPs in Different Organism

Candidate FABP family genes in the other four main milk-supplying livestock species including cattle (21), sheep (18), goat (18), and horse (23) were identified using the same method as in buffalo. The candidate sequence together with buffalo FABPs, a total of 106 amino acid sequences, were aligned by MUSCLE [34] and a neighbor-joining tree and maximum likelihood tree were assembled in MEGA 7.0 with a bootstrap test implemented with 1000 replications (random seed) [35].

2.3. Structural Features Analysis

MEGA 7.0 was used to analyze phylogenetic relationships of all FABP amino acid sequences identified to confirm that they belong to the buffalo FABP family genes in order to cluster the FABP amino acid sequences into different subfamilies [36]. To further access the structural diversity of buffalo FABP family members, the sequences were submitted to the MEME 5.0 website (https://meme-suite.org/meme/, accessed on 1 September 2021) to detect conserved motifs [37] and visualized by TBtools (https://github.com/CJ-Chen/TBtools, accessed on 1 September 2021) [38]. To further analyze the detail and function of motifs from the MEME analysis, motif sequences were submitted to the Pfam database (http://pfam.xfam.org/, accessed on 1 September 2021) to search their Pfam matches [39]. The amino acid sequences of the buffalo FABP family were also submitted to Conversed Domains Database (https://www.ncbi.nlm.nih.gov/cdd, accessed on 1 September 2021) to search conversed domains [39] and visualized by TBtools.

2.4. Chromosomal Distribution and Gene Duplication Analysis

Gene duplication events with the default parameters of all the identified FABP family members were detected by the Multiple Collinearity Scan Toolkit (MCScanX) software [40]. Then, a chromosome position map with a duplication relationship was constructed using TBtools with the buffalo annotation GFF3 file [38]. The divergence times of duplicated FABP gene pairs were calculated by divergence time = Ks/2λ [41], where Ks is the synonymous substitution rate and λ = 1.26 × 10−8 is the clock-like rate in buffalo [42]. To exhibit the positive selection relationship of each tandem duplicated FABP gene pair in buffalo, their synonymous (Ks) and nonsynonymous (Ka) numbers of substitution analysis were performed using TBtools software [38]. Finally, a collinearity map was constructed to detect orthologous gene pairs in buffalo and cattle.

2.5. Association Analysis of SNP and Buffalo Milk Traits

Phenotype resources, including 1424 lactation records of 489 Mediterranean Italian buffalos born between 2000 to 2011 and reared in 4 herds in southern Italy, were reported in our previous study [43]. Single nucleotide polymorphisms (SNPs) within and 2000 base pairs upstream of the buffalo FABP family genes were obtained from the genotyping data conducted at Delta Genomics (Edmonton, AB, Canada) using the 90K Axiom Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA, USA) [43]. Haploview 4.2 was used to compute predicted heterozygosity (He) and Hardy–Weinberg equilibrium (HWE) of the identified SNPs [44]. By using the PIC CALC software, we calculated the polymorphism information content (PIC) of each locus [45]. Six lactation traits comprising peak milk yield, total milk yield, total protein yield, protein percentage, total fat yield, and fat percentage were adjusted to 270-day records to eliminate the effects of environmental factors. The association between each SNP with six 270-day adjusted buffalo milk production traits was performed by the lme4 R-package with the lmer procedure using the following model: Yijkl = μ + Bi + Pj + HYSk + Gl + eijkl, where Yijkl = trait observation, μ = overall mean, Bi = random effect of buffalo individual, Pj = fixed effect of the jth parity (nine classes, 1–9), HYSk = fixed effect of contemporary group constructed with the effects of herd-season and year, Gl = fixed effect of the lth genotype, and eijkl = random residual [46] For the pairwise comparisons among different levels of fixed effects included in the model, the Bonferroni correction for multiple F-testing was applied to the least square means ± SE of the findings for FABPs genotypes.

2.6. Cell Culture and Fatty Acid Treatment

The cell culture and medium preparation to culture buffalo mammary epithelial cells (BuMECs) and bovine mammary epithelial cells (BoMECs) were followed with brief modifications as reference [47]. Cell culture medium was prepared using DMEM/F12 (Hyclone, Logan, UT, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco, Billings, MT, USA), 5 ug/mL bovine insulin (Sigma, Ronkonkoma, NY, USA), 1 µg/mL hydrocortisone (Sigma, Ronkonkoma, NY, USA), 1 µg/mL transferrin (Sigma, Ronkonkoma, NY, USA), 1 µg/mL progesterone (Solarbio, Beijing, China), 10 ng/mL EGF (Proteintech, Rosemont, IL, USA), 100 U/mL penicillin, and 100 mg/mL streptomycin. BuMECs were cultured in an incubator at 37 °C under 5% CO2. For fatty acid treatment, BuMECs and BoMECs were seeded at 1 × 105/mL into a 6-well plate 48 h before treatment. Three kinds of fatty acids including palmitic acid, stearic acid, and oleic acid (Sigma, Ronkonkoma, NY, USA) were first dissolved with 20 mg/mL in ethanol and then mixed with complete medium to a final concentration of 25 µM, 50 µM, 75 µM, 100 µM, and 125 µM, respectively. Then, the cells were harvested after 6 h in an incubator.

2.7. Isolation and Culture of BuMECs

A similar protocol to that used by Vijay Anand et al. was followed to isolate buffalo mammary gland epithelial cells (BuMECs) [47]. In summary, after slaughter, macerated parenchymal tissues from disease-free buffaloes were harvested and transported to the laboratory in DMEM containing 5 µg/mL of streptomycin, 100 U/m of penicillin, and 50 ng/mL of amphotericin. Uteri were washed three times using PBS and trimmed connective tissue and fat. The tissue was cut into small pieces and digested by 0.05% collagenase (Sigma, Ronkonkoma, NY, USA) and 0.05% Hyaluronidase (Sigma, Ronkonkoma, NY, USA) for 5 h at 37 °C. For further digestion, trypsin-EDTA (0.25%) and Dispase (1%; Proteintech, Rosemont, IL, USA) were added to the digested tissue, which was then incubated for 30 min at 37 °C, and the extract was filtered through a 40-cell strainer (Whb-bio, Shanghai, China). The filtrate was centrifuged at 80× g for 5 min and washed thrice by PBS. BuMECs were cultured in a culture plate (Corning, Corning, NY, USA) in a growth medium, as described in [47]. After incubation for 5 days, selective trypsinization steps were used to remove the fibroblast cells and get purified BuMECs. The cells were treated 0.25% trypsin-EDTA (Gibco, Billings, MT, USA) and trypsinization was stopped immediately after incubation for 3 min. The detached fibroblast cells in the supernatant were removed. The purified BuMECs were suspended in a freezing medium constituting 90% FBS (Gibco, Billings, MT, USA) and 10% DMSO (Sigma, Ronkonkoma, NY, USA) for the next experiments.
For LCFAs treatment, a density of 1 × 105 cells/cm2 was cultured in 6-well plates (Thermo Fisher, Waltham, MA, USA). The medium was changed after one day containing 0, 25 µM, 50 µM, 75 µM, and 100 µM LACFs including palmitic acid, stearic acid, and oleic acid. Then, the cells were allowed to grow for one day and harvested for RNA extraction.

2.8. qRT-PCR Analysis

According to the manufacturer’s instructions, cells were extracted with a total RNA kit. Extracted RNA was evaluated using agarose gel electrophoresis and a Nanodrop 2000 spectrophotometer. A total of 2 μg DNase-treated RNA was used for the cDNA synthesis. The qRT-PCR was conducted on the Bio-Rad CFX Maestro (BIO-RAD, Philadelphia, PA, USA) instrument using the QuantiNova SYBR Green PCR Kit (QIAGEN, Hilden, Germany). The GAPDH gene was used as an internal control. The 2ΔΔCT method was used to calculate the expression level [48]. The primers used are listed in Table S1.

3. Results

3.1. Genome Identification of FABP Family Members

A total of 26 non-redundant buffalo proteins encoded by 17 genes were identified using the BLAST program and HMMER software (Table 1). These genes include 2 genes encoding the cellular retinoic acid-binding protein (CRABP), 9 genes encoding the FABPs, 4 genes encoding the retinol-binding protein, a gene encoding the myelin P2 protein, and an uncharacterized protein. The amino acid length of FABP protein isoforms ranged from 116 to 348, with the predicted molecular weight from 13.16 kDa to 39.13 kDa. The isoelectric point of the transforms ranged from 4.88 to 5.76.

3.2. Structural Features of Buffalo FABP Family Members

A phylogenetic analysis was performed on buffalo FABPs to analyze their motif patterns, gene structures, and conserved domains. The phylogenetic analysis showed that buffalo FABPs can be divided into 5 groups (Figure 1A). Ten motifs were identified in buffalo FABPs by MEME research (Figure 1B). Motifs 1 and 2, which were annotated as the lipocalin/cytosolic fatty-acid binding protein family by Pfam search, were the most common motifs around the family members (Table S2). The gene structural analysis showed that almost all genes except the uncharacterized gene had a complete genome no longer than 30,000 bp. Meanwhile, the FABPs in the same group had a similar structure and exon number (Figure 1C). Furthermore, the conserved domain analysis showed that most of the family members contain domains related to the FABP superfamily (Figure 1D).

3.3. Phylogenetic Relationship Analysis of FABP Protein in Five Mammals

To assess the evolutionary relationship of the FABP protein in buffalo and other mammals, a neighbor-joining tree map was built by MUSCLE alignment [34] using all 106 identified protein sequences (Figure 2). Equivalent to buffalo, the FABP protein can be divided into 5 groups. The first group included the most FABP protein sequences (n = 45), while the fourth group contained only one FABP gene, including five protein sequences. Besides, by comparing the evolutionary relationship between different organisms, we found that buffalo had a particularly closer relationship with cattle than with other species.

3.4. Chromosomal Distribution and Collinearity Analysis of FABP Genes

Based on the mapping data of buffalo with the duplication events analysis, all identified FABP genes were randomly distributed across 11 chromosomes (Figure 3A). Most of FABP family members had 5 genes distributed very closely on chromosome 15. Meanwhile, except for different isoforms of the same genes, many duplication events occurred in different genes. Interestingly, 3 pairs of tandem duplication genes, including FABP12-FABP4, FABP4-FABP9, and FABP9-MP2P, were distributed on chromosome 15. In addition, the uncharacterized genes LOC102401361 and RBP2 were also tandem duplication pairs. For the tandem duplication events of FABP12-FABP4, FABP4-FABP9, and FABP9-MP2P, the value of the nonsynonymous substitutions per nonsynonymous site (Ka)/the number of synonymous substitutions per synonymous site (Ks) ratios ranged from 0.104 to 0.209 (Table 2), indicating that the FABP family genes might have experienced intense purifying selective pressure during evolution. For the tandem duplication events of LOC102401361-RBP2 and LOC102401361-MP2P, the Ka/Ks ratios were more than 1, showing that LOC102401361 might be under positive selection. The divergence time of duplicated FABP pairs ranged from 51.328 to 98.92 Mya.
Because the buffalo FABP genes had a particularly closer relationship with cattle than with other species, we further performed a collinearity analysis between cattle and buffalo. The results indicated that 34,710 collinear gene pairs were identified, which account for 81.95% of the total genes (Figure 3B). Although buffalo have different chromosome numbers than cattle, the syntenic gene blocks covered almost all of their chromosomes. Especially, we found 13 pairs of orthologous genes between buffalo and cattle (Table 3). Our findings suggested that these orthologous genes might have a conserved function between the two species during evolution.

3.5. Analyses of Association between Traits Related to Buffalo Milk Production

Based on our previous study on 489 buffalo with 1424 lactation records and 60,387 SNPs after quality control for individuals with a call rate ≥97%., we used genotypic and phenotypic data to identify potential markers or genes influencing milk traits in buffalo FABP genes [43]. According to the genotyping dataset, a total of 7 SNPs within 3 FABP genes were filtered after quality control and used in the present study. The allelic and genotypic frequencies, observed and expected heterozygosity, p-value of HWE, and PIC are shown in Table 4.
For each SNP, we executed an association analysis with the six milk production traits, then the AX-85106417 located on LOC102401361 was found significantly associated with MY270. Moreover, the least square mean of individuals with GG (2700 ± 126 kg) was significantly lower than those with AA (3047 ± 70 kg) and AG (3074 ± 72 kg) (Table 5).

3.6. Effect of LCFAs on Expression of FABPs

Here, we found that CRABP1, CRABP2, FABP3, FABP4, FABP5, FABP7, and PMP2 were highly expressed in buffalo mammary epithelial cells. After palmitic acid treatment, CRABP1, CRABP2, FABP3, FABP4, and FABP5 mRNA levels were upregulated dramatically with the increase in concentrations, then the FABP7 mRNA level increased after adding 25, 50, and 75 μM palmitic acid, and the MP2P was almost unaffected. As for stearic acid treatment, CRABP2, FABP3, FABP4, and FABP5 mRNA levels increased dramatically with the increase in concentrations, but the MP2P declined dramatically with the increase in concentrations; the mRNA level of CRABP1 reached the lowest when the concentration was 50 uM. It was also found that stearic acid had little effect on FABP7. However, oleic acid treatment had no effect or minimal effect on the expression of FABPs (Figure 4).

4. Discussion

Buffalo are important livestock in the agricultural economy because they supply milk, meat, and draught for plowing [49]. However, low milk yield has seriously restricted the development of the buffalo industry [50]. It is known that FABPs bind free fatty acids for transport to different organelles for lipid metabolism [51]. Besides, it has also been found that the FABP protein plays an essential role in bovine milk fat synthesis [14]. FABP family proteins have been well characterized in some nonmammalian vertebrates, such as teleost fishes [9], and chicken [1]. However, as important genes related to the synthesis of milk fat, available information on the FABP family genes in milk-producing livestock species, especially in water buffalo, is still limited. The present study identified 26 FABP family protein sequences encoded by 17 genes based on the complete buffalo genome sequence. The identified FABPs were classified into five groups according to their structural features and evolutionary relationships. Consistent with a previous study of cellular FABPs [5,12], this superfamily also included cellular retinol-binding proteins (CRBPs) and cellular retinoic acid-binding proteins (CRABPs) subfamilies. In addition, our motif and conserved domain analysis also proved the classification. We searched 10 motifs from all identified FABP protein sequences, three of which matched to the lipocalin/cytosolic fatty-acid binding protein family, while most of the protein sequences shared at least one of them. All identified protein sequences contained at least one conserved domain related to fatty acid binding. The results were also supported by a previous study [12]. Moreover, we also identified a gene, myelin P2 protein, in the buffalo annotation file; its protein structure is similar to the conserved structure of FABPs [11], and was therefore considered as FABP8 in several organisms such as humans [12,52] and bovine [53]. Thus, it was reasonable to determine the myelin P2 protein as a member of the FABP family genes.
The phylogenetic analysis of FABP family genes displayed a deep understanding of the evolution among them. From the neighbor-joining tree, it was more clear that all FABP family genes could be classified into five groups which is consistent with the result observed in buffalo. Segmental and tandem duplication were common phenomena providing a possibility for novel gene function acquisition in genome evolution [54,55]. A total of 17 identified buffalo FABP family genes were distributed across 13 chromosomes. Here, we found four tandem duplication pairs, including LOC102401361/RBP2, FABP9/MP2P, FABP4/FABP9, and FABP12/FABP4. For them, the ratio of three pairs of tandem duplication genes (FABP9/MP2P, FABP4/FABP9, and FABP12/FABP4) was less than 0.3, which indicated that these genes might experience strong purifying selection pressure.
Moreover, we observed that the value of the Ka/Ks ratio for LOC102401361/RBP2 was more than 1, indicating that the LOC102401361 gene might be under positive selection during evolution. It is well known that the ratio of Ka and Ks is usually used to determine whether there is selection pressure acting on the protein-coding gene [56,57,58]. We also found that the divergence time of this duplicated pair occurred at 58.515 Mya. These results suggest that the LOC102401361 gene has the potential to form a novel function. However, this needs to be confirmed further.
To identify more FABPs that affect buffalo milk production traits, we constructed a mixed linear model to perform SNP-traits association analysis. Many SNPs detected in different FABP genes have been significantly associated with milk traits in dairy cows, such as FABP4 and FABP3 [17]. However, the present study did not detect relevant genetic variations within these genes due to the minor SNP density of the genotyping arrays. Consequently, only seven SNPs were found in three FABP members, and only one of them was suggested as a potential marker for screening buffalo milk performance. The identification of a significant association for the AX-85106417 marker in this work could indicate that this gene has a potential influence on buffalo milk yield. However, the limited number of animals analysed in this study will make necessary that further research confirm the potential effect of polymorphism on the LOC102401361 gene on buffalo milk yield.
It is well known that mammary gland epithelial cells play a vital role in the biosynthesis of milk fat in buffalo. Several FABPs have reported to show a high expression in bovine mammary gland tissue and are remarkably up-regulated during lactation [59]. As reported in previous studies, several FABP members served as functional genes for milk fat percentage by regulating the transportation of fatty acid [17,60]. In the present study, we found high expression of CRABP1, CRABP3, FABP3, FABP4, FABP5, FABP7, and MP2P in buffalo mammary epithelial cells, and most of them were regulated by palmitic acid and stearic acid treatment but not oleic acid treatment, demonstrating that these genes may be involved in the uptake and metabolism of these two kinds of fatty acids. It has been also reported that sheep FABPs were involved in the cellular uptake and metabolism of LCFAs [61]. For example, other study have shown that the expression of FABP3 in cow mammary epithelial cells was increased after adding 50 and 75 μM oleic acid, 100 μM stearic acid, and 125 μM palmitic acid [16] and it was also affected by short-chain fatty acids (SCFAs) in goat mammary epithelial cells [62]. Those results are in agreement with the results here reported, which suggest that palmitic acid and stearic acid could up-regulate the expression level of FABP3 in buffalo mammary gland epithelial cells. Furthermore, the expression of several FABPs was found to be affected by palmitic acid and stearic acid, suggesting that FABPs may regulate the metabolism of milk fat by influencing the uptake of LCFAs.

5. Conclusions

A total of 17 FABP genes were identified in buffalo genome. Phylogenetic analysis performed here classified these sequences into five groups, showing similar motifs and gene structures within each group. Our analysis suggested that FABPs underwent purifying selection in buffalo and cattle during the evolutionary process. Only one SNP in LOC102401361 was found significantly associated with buffalo milk yield due to limited number of animals analysed in this study. LOC102401361 was also found under positive selection during evolution, indicating that the uncharacterized protein might be involved in milk synthesis. The results described in this study showed that the expression of most FABPs in buffalo mammary epithelial cells was regulated by palmitic and stearic acid treatment. The current study provides preliminary results for further investigation into the potential roles that FABP family genes play in the regulation of buffalo milk synthesis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13040600/s1, Table S1: Primers used for the qRT-PCR analysis; Table S2: Ten different motifs commonly observed in the buffalo FABP family.

Author Contributions

T.Y. and L.Y. (Liguo Yang) designed the study; T.Y., L.Y. (Lv Yang) and Y.Z. analyzed the data; T.Y. and L.Y. (Lv Yang) interpreted the data; T.Y. and C.C. wrote the drafted manuscript; A.S., L.Y. (Liguo Yang) reviewed the paper. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by China Agriculture Research System of MOF and MARA (CARS-36) and Fundamental Research Funds for the Central Universities (2662019PY014).

Institutional Review Board Statement

Experiments regarding animals were carried out after being approved by the institutional research ethics and animal welfare committee, Huazhong Agricultural University (HZAUMO_2015-12) (Wuhan, China).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Z.; Yue, Y.X.; Liu, Z.M.; Yang, L.Y.; Li, H.; Li, Z.J.; Li, G.X.; Wang, Y.B.; Tian, Y.D.; Kang, X.T.; et al. Genome-Wide Analysis of the FABP Gene Family in Liver of Chicken (Gallus gallus): Identification, Dynamic Expression Profile, and RegulatoryMechanism. Int. J. Mol. Sci. 2019, 20, 5948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Marion, M.; Hamilton, J.; Richardson, B.; Roeder, N.; Figueiredo, A.; Nubelo, A.; Hetelekides, E.; Penman, S.; Owada, Y.; Kagawa, Y.; et al. Environmental enrichment sex-dependently rescues memory impairment in FABP5 KO mice not mediated by brain-derived neurotrophic factor. Behav. Brain Res. 2022, 425, 113814. [Google Scholar] [CrossRef] [PubMed]
  3. Donnelly, K.L.; Smith, C.I.; Schwarzenberg, S.J.; Jessurun, J.; Boldt, M.D.; Parks, E.J. Sources of fatty acids stored in liver and secreted via lipoproteins in patients with nonalcoholic fatty liver disease. J. Clin. Investig. 2005, 115, 1343–1351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zheng, Y.; Blair, D.; Bradley, J.E. Phyletic distribution of fatty acid-binding protein genes. PLoS ONE 2013, 8, e77636. [Google Scholar] [CrossRef] [Green Version]
  5. Glatz, J.F.; van der Vusse, G.J. Cellular fatty acid-binding proteins: Their function and physiological significance. Prog. Lipid Res. 1996, 35, 243–282. [Google Scholar] [CrossRef]
  6. Venkatachalam, A.B.; Thisse, C.; Thisse, B.; Wright, J.M. Differential tissue-specific distribution of transcripts for the duplicated fatty acid-binding protein 10 (fabp10) genes in embryos, larvae and adult zebrafish (Danio rerio). FEBS J. 2009, 276, 6787–6797. [Google Scholar] [CrossRef]
  7. Agulleiro, M.J.; André, M.; Morais, S.; Cerdà, J.; Babin, P.J. High transcript level of fatty acid-binding protein 11 but not of very low-density lipoprotein receptor is correlated to ovarian follicle atresia in a teleost fish (Solea senegalensis). Biol. Reprod. 2007, 77, 504–516. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, R.Z.; Li, X.; Godbout, R. A novel fatty acid-binding protein (FABP) gene resulting from tandem gene duplication in mammals: Transcription in rat retina and testis. Genomics 2008, 92, 436–445. [Google Scholar] [CrossRef] [Green Version]
  9. Parmar, M.B.; Venkatachalam, A.B.; Wright, J.M. The evolutionary relationship of the transcriptionally active fabp11a (intronless) and fabp11b genes of medaka with fabp11 genes of other teleost fishes. FEBS J. 2012, 279, 2310–2321. [Google Scholar] [CrossRef]
  10. Schaap, F.G.; van der Vusse, G.J.; Glatz, J.F. Evolution of the family of intracellular lipid binding proteins in vertebrates. Mol. Cell. Biochem. 2002, 239, 69–77. [Google Scholar] [CrossRef]
  11. Furuhashi, M.; Hotamisligil, G.S. Fatty acid-binding proteins: Role in metabolic diseases and potential as drug targets. Nat. Rev. Drug Discov. 2008, 7, 489–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ono, T.; Odani, S. Initial studies of the cytoplasmic FABP superfamily. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2010, 86, 220–228. [Google Scholar] [CrossRef] [Green Version]
  13. Chmurzyńska, A. The multigene family of fatty acid-binding proteins (FABPs): Function, structure and polymorphism. J. Appl. Genet. 2006, 47, 39–48. [Google Scholar] [CrossRef] [PubMed]
  14. Bionaz, M.; Loor, J.J. Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics 2008, 9, 366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Whetstone, H.D.; Hurley, W.L.; Davis, C.L. Identification and characterization of a fatty acid binding protein in bovine mammary gland. Comp. Biochem. Physiol. B 1986, 85, 687–692. [Google Scholar] [CrossRef]
  16. Liang, M.Y.; Hou, X.M.; Qu, B.; Zhang, N.; Li, N.; Cui, Y.J.; Li, Q.Z.; Gao, X.J. Functional analysis of FABP3 in the milk fat synthesis signaling pathway of dairy cow mammary epithelial cells. Vitr. Cell. Dev. Biol. Anim. 2014, 50, 865–873. [Google Scholar] [CrossRef]
  17. Nafikov, R.A.; Schoonmaker, J.P.; Korn, K.T.; Noack, K.; Garrick, D.J.; Koehler, K.J.; Minick-Bormann, J.; Reecy, J.M.; Spurlock, D.E.; Beitz, D.C. Association of polymorphisms in solute carrier family 27, isoform A6 (SLC27A6) and fatty acid-binding protein-3 and fatty acid-binding protein-4 (FABP3 and FABP4) with fatty acid composition of bovine milk. J. Dairy Sci. 2013, 96, 6007–6021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Cho, S.; Park, T.S.; Yoon, D.H.; Cheong, H.S.; Namgoong, S.; Park, B.L.; Lee, H.W.; Han, C.S.; Kim, E.M.; Cheong, I.C.; et al. Identification of genetic polymorphisms in FABP3 and FABN and putative association with back fat thickness in Korean native cattle. BMB Rep. 2008, 41, 29–34. [Google Scholar] [CrossRef] [Green Version]
  19. Zhihua, J.; Michal, J.J. Polymorphisms in Fatty Acid Binding Protein 4 (“fabp4”) Gene and their Associations with Measures of Marbling and Subcutaneous Fat Depth in Beef Cattle. US 2007/0020658 A1. 2007. Available online: https://patentscope2.wipo.int/search/en/detail.jsf;jsessionid=C98272504AA78E3E22B622E9D1098F7B.wapp1nC?docId=AU181382563_cid=P12-K6CZY1-99124-1 (accessed on 1 September 2021).
  20. Barendse, W.; Bunch, R.J.; Thomas, M.B.; Harrison, B.E. A splice site single nucleotide polymorphism of the fatty acid binding protein 4 gene appears to be associated with intramuscular fat deposition in longissimus muscle in Australian cattle. Anim. Genet. 2009, 40, 770–773. [Google Scholar] [CrossRef]
  21. Zhang, C.; Yan, C.; Xu, A.; Chan, B.; Chiu, K. Knocking out or pharmaceutical inhibition of fatty acid binding protein 4 (fabp-4) alleviates osteoarthritis induced by a very high fat diet in mouse. OsCar 2017, 25, S14. [Google Scholar] [CrossRef]
  22. Majchrzak, K.; Piotrowska, M.; Krajewska, J.; Fichna, J. Adipocyte fatty acid binding protein (A-FABP) as a potential new therapeutic target for treatment of obesity—associated cancers. Curr. Drug Targets 2021, 22, 1. [Google Scholar] [CrossRef] [PubMed]
  23. Ghaffari, M.H.; Alaedin, M.T.; Sadri, H.; Hofs, I.; Koch, C.; Sauerwein, H. Longitudinal changes in fatty acid metabolism and in the mitochondrial protein import system in overconditioned and normal conditioned cows: A transcriptional study using microfluidic quantitative PCR. J. Dairy Sci. 2021, 104, 10338–10354. [Google Scholar] [CrossRef] [PubMed]
  24. Chang, W.; Rickers-Haunerland, J.; Haunerland, N.H. Induction of cardiac FABP gene expression by long chain fatty acids in cultured rat muscle cells. Mol. Cell. Biochem. 2001, 221, 127–132. [Google Scholar] [CrossRef] [PubMed]
  25. van Breda, E.; Keizer, H.A.; Vork, M.M.; Surtel, D.A.; de Jong, Y.F.; van der Vusse, G.J.; Glatz, J.F. Modulation of fatty-acid-binding protein content of rat heart and skeletal muscle by endurance training and testosterone treatment. Pflugers Arch. 1992, 421, 274–279. [Google Scholar] [CrossRef]
  26. van der Lee, K.A.; Vork, M.M.; De Vries, J.E.; Willemsen, P.H.; Glatz, J.F.; Reneman, R.S.; Van der Vusse, G.J.; Van Bilsen, M. Long-chain fatty acid-induced changes in gene expression in neonatal cardiac myocytes. J. Lipid Res. 2000, 41, 41–47. [Google Scholar] [CrossRef]
  27. Johnsen, G.M.; Weedon-Fekjaer, M.S.; Tobin, K.A.; Staff, A.C.; Duttaroy, A.K. Long-chain polyunsaturated fatty acids stimulate cellular fatty acid uptake in human placental choriocarcinoma (BeWo) cells. Placenta 2009, 30, 1037–1044. [Google Scholar] [CrossRef]
  28. Kadegowda, A.K.; Bionaz, M.; Piperova, L.S.; Erdman, R.A.; Loor, J.J. Peroxisome proliferator-activated receptor-gamma activation and long-chain fatty acids alter lipogenic gene networks in bovine mammary epithelial cells to various extents. J. Dairy Sci. 2009, 92, 4276–4289. [Google Scholar] [CrossRef] [Green Version]
  29. Khedkar, C.; Kalyankar, S.; Deosarkar, S. Buffalo milk. In Encyclopedia of Food and Health; Academic Press: Waltham, MA, USA, 2016; pp. 522–528. [Google Scholar] [CrossRef]
  30. Liu, J.; Wang, Z.; Li, J.; Li, H.; Yang, L. Genome-wide identification of Diacylglycerol Acyltransferases (DGAT) family genes influencing Milk production in Buffalo. BMC Genet. 2020, 21, 26. [Google Scholar] [CrossRef] [Green Version]
  31. Mintoo, A.A.; Zhang, H.; Chen, C.; Moniruzzaman, M.; Deng, T.; Anam, M.; Emdadul Huque, Q.M.; Guang, X.; Wang, P.; Zhong, Z. Draft genome of the river water buffalo. Ecol. Evol. 2019, 9, 3378–3388. [Google Scholar] [CrossRef]
  32. Li, Y.; He, L.; Li, J.; Chen, J.; Liu, C. Genome-Wide Identification, Characterization, and Expression Profiling of the Legume BZR Transcription Factor Gene Family. Front. Plant Sci. 2018, 9, 1332. [Google Scholar] [CrossRef]
  33. Potter, S.C.; Luciani, A.; Eddy, S.R.; Park, Y.; Lopez, R.; Finn, R.D. HMMER web server: 2018 update. Nucleic Acids Res. 2018, 46, W200–W204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Strimmer, K.; von Haeseler, A. Quartet Puzzling: A Quartet Maximum-Likelihood Method for Reconstructing Tree Topologies. Mol. Biol. Evol. 1996, 13, 964. [Google Scholar] [CrossRef] [Green Version]
  37. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  38. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  39. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  41. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [Green Version]
  42. Ma, X.; Liang, S.; Liang, A.; Rushdi, H.E.; Deng, T. Evolutionary Analysis of OAT Gene Family in River and Swamp Buffalo: Potential Role of SLCO3A1 Gene in Milk Performance. Genes 2021, 12, 1394. [Google Scholar] [CrossRef]
  43. Liu, J.; Liang, A.; Campanile, G.; Plastow, G.; Zhang, C.; Wang, Z.; Salzano, A.; Gasparrini, B.; Cassandro, M.; Yang, L. Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo. J. Dairy Sci. 2018, 101, 433–444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhang, X.C.; Zhao, J.; Li, W.; Liu, Y.H.; Wei, C.Q.; Wang, Y.K.; Ye, Z.Y.; Zhu, X.P. Development and characterization of 27 SNP markers in the Mauremys mutica transcriptome. Conserv. Genet. Resour. 2018, 10, 667–670. [Google Scholar] [CrossRef]
  46. Ye, T.; Deng, T.; Hosseini, S.M.; Raza, S.H.A.; Du, C.; Chen, C.; Zhang, X.; Hu, X.; Yang, L. Association analysis between FASN genotype and milk traits in Mediterranean buffalo and its expression among different buffalo tissues. Trop. Anim. Health Prod. 2021, 53, 366. [Google Scholar] [CrossRef]
  47. Anand, V.; Dogra, N.; Singh, S.; Kumar, S.N.; Jena, M.K.; Malakar, D.; Dang, A.K.; Mishra, B.P.; Mukhopadhyay, T.K.; Kaushik, J.K. Establishment and characterization of a buffalo (Bubalus bubalis) mammary epithelial cell line. PLoS ONE 2012, 7, e40469. [Google Scholar] [CrossRef] [Green Version]
  48. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  49. Warriach, H.; McGill, D.; Bush, R.; Wynn, P.; Chohan, K. A review of recent developments in buffalo reproduction—A review. Asian-Australas. J. Anim. Sci. 2015, 28, 451. [Google Scholar] [CrossRef]
  50. Du, C.; Deng, T.; Zhou, Y.; Ye, T.; Zhou, Z.; Zhang, S.; Shao, B.; Wei, P.; Sun, H.; Khan, F. Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis). Anim. Genet. 2019, 50, 207–216. [Google Scholar] [CrossRef]
  51. Sánchez-Gurmaches, J.; Cruz-Garcia, L.; Gutiérrez, J.; Navarro, I. mRNA expression of fatty acid transporters in rainbow trout: In vivo and in vitro regulation by insulin, fasting and inflammation and infection mediators. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2012, 163, 177–188. [Google Scholar] [CrossRef]
  52. Zenker, J.; Stettner, M.; Ruskamo, S.; Domènech-Estévez, E.; Baloui, H.; Médard, J.J.; Verheijen, M.H.; Brouwers, J.F.; Kursula, P.; Kieseier, B.C.; et al. A role of peripheral myelin protein 2 in lipid homeostasis of myelinating Schwann cells. Glia 2014, 62, 1502–1512. [Google Scholar] [CrossRef]
  53. Kitamura, K.; Suzuki, M.; Suzuki, A.; Uyemura, K. The complete amino acid sequence of the P2 protein in bovine peripheral nerve myelin. FEBS Lett. 1980, 115, 27–30. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, G.E.; Ventura, M.; Cellamare, A.; Chen, L.; Cheng, Z.; Zhu, B.; Li, C.; Song, J.; Eichler, E.E. Analysis of recent segmental duplications in the bovine genome. BMC Genom. 2009, 10, 571. [Google Scholar] [CrossRef] [Green Version]
  55. Liao, X.; Bao, H.; Meng, Y.; Plastow, G.; Moore, S.; Stothard, P. Sequence, structural and expression divergence of duplicate genes in the bovine genome. PLoS ONE 2014, 9, e102868. [Google Scholar] [CrossRef] [PubMed]
  56. Zhang, Z.; Li, J.; Zhao, X.Q.; Wang, J.; Wong, G.K.; Yu, J. KaKs_Calculator: Calculating Ka and Ks through model selection and model averaging. Genom. Proteom. Bioinform. 2006, 4, 259–263. [Google Scholar] [CrossRef] [Green Version]
  57. Liu, J.; Wang, Y.; Zhang, R.; Shi, C.; Lu, W.; Li, J.; Bai, M. Three Complete Mitochondrial Genomes of Erotylidae (Coleoptera: Cucujoidea) with Higher Phylogenetic Analysis. Insects 2021, 12, 524. [Google Scholar] [CrossRef] [PubMed]
  58. Mori, S.; Matsunami, M. Signature of positive selection in mitochondrial DNA in Cetartiodactyla. Genes Genet. Syst. 2018, 93, 65–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Bionaz, M.; Loor, J.J. ACSL1, AGPAT6, FABP3, LPIN1, and SLC27A6 are the most abundant isoforms in bovine mammary tissue and their expression is affected by stage of lactation. J. Nutr. 2008, 138, 1019–1024. [Google Scholar] [CrossRef]
  60. Cohen-Zinder, M.; Lipkin, E.; Shor-Shimoni, E.; Ben-Meir, Y.; Agmon, R.; Asher, A.; Miron, J.; Shabtay, A. FABP4 gene has a very large effect on feed efficiency in lactating Israeli Holstein cows. Physiol. Genom. 2019, 51, 481–487. [Google Scholar] [CrossRef]
  61. Calvo, J.H.; Marcos, S.; Jurado, J.J.; Serrano, M. Association of the heart fatty acid-binding protein (FABP3) gene with milk traits in Manchega breed sheep. Anim. Genet. 2004, 35, 347–349. [Google Scholar] [CrossRef]
  62. Sun, Y.; Luo, J.; Zhu, J.; Shi, H.; Li, J.; Qiu, S.; Wang, P.; Loor, J.J. Effect of short-chain fatty acids on triacylglycerol accumulation, lipid droplet formation and lipogenic gene expression in goat mammary epithelial cells. Anim. Sci. J. 2016, 87, 242–249. [Google Scholar] [CrossRef]
Figure 1. Representations of the recognized FABP proteins and genes isoforms in buffalo. (A) The phylogenetic tree was constructed by the N-J method. (B) The structure of amino acid sequence boxes represents ten conserved motifs. In the gene structure plot (C), the green box, black line, and orange box represent untranslated region (UTR), intron, and coding sequencing (CDS), respectively. Conserved domains (D) were searched by NIBI-CCD search.
Figure 1. Representations of the recognized FABP proteins and genes isoforms in buffalo. (A) The phylogenetic tree was constructed by the N-J method. (B) The structure of amino acid sequence boxes represents ten conserved motifs. In the gene structure plot (C), the green box, black line, and orange box represent untranslated region (UTR), intron, and coding sequencing (CDS), respectively. Conserved domains (D) were searched by NIBI-CCD search.
Genes 13 00600 g001
Figure 2. Phylogenetic relationship of FABP proteins in five mammals. Different colored lines indicate different groups. Circles with different colors indicate different groups.
Figure 2. Phylogenetic relationship of FABP proteins in five mammals. Different colored lines indicate different groups. Circles with different colors indicate different groups.
Genes 13 00600 g002
Figure 3. Gene duplication of buffalo FABP genes (A) and collinear analysis between buffalo and cattle (B). A red line marks the orthologous FABP genes pairs between buffalo and cattle.
Figure 3. Gene duplication of buffalo FABP genes (A) and collinear analysis between buffalo and cattle (B). A red line marks the orthologous FABP genes pairs between buffalo and cattle.
Genes 13 00600 g003
Figure 4. Relative mRNA levels of FABPs in buffalo mammary gland epithelial cells at 24 h after adding palmitic acid, stearic acid, or oleic acid (0, 50, 75, 100, 125, 150 μM), respectively.
Figure 4. Relative mRNA levels of FABPs in buffalo mammary gland epithelial cells at 24 h after adding palmitic acid, stearic acid, or oleic acid (0, 50, 75, 100, 125, 150 μM), respectively.
Genes 13 00600 g004
Table 1. Features of the identified FABP protein sequences in buffalo.
Table 1. Features of the identified FABP protein sequences in buffalo.
ListProtein IsoformGene IDProtein IDAmino AcidsIsoelectric PointMw/kDaProduct
1CRABP1102392457XP_006044444.11375.2615.59Cellular retinoic acid-binding protein 1
2CRABP2102406669XP_006052022.11385.3715.73Cellular retinoic acid-binding protein 2
3FABP1102407733XP_006074804.11277.7814.19Fatty acid-binding protein%2C liver
4FABP12102410779XP_006068112.21215.5113.8Fatty acid-binding protein 12
5FABP2102414836XP_006067721.11325.9415.09Fatty acid-binding protein%2C intestinal
6FABP3.1102394447NP_001277811.11336.7314.77Fatty acid-binding protein%2C heart
7FABP3.2102394447XP_025150979.11718.6318.91Fatty acid-binding protein%2C heart isoform X1
8FABP4102410448NP_001277890.11325.0414.76Fatty acid-binding protein%2C adipocyte
9FABP5102409117XP_006068107.11357.5815.07Fatty acid-binding protein%2C epidermal
10FABP6102412445XP_006073509.11286.9114.37gastrotropin
11FABP7-X1102409019XP_006047648.11325.3814.95Fatty acid-binding protein%2C brain isoform X1
12FABP7-X2102409019XP_025150573.11185.1713.42Fatty acid-binding protein%2C brain isoform X2
13FABP7-X3102409019XP_006047649.11165.1713.16Fatty acid-binding protein%2C brain isoform X3
14FABP9102410109XP_006068110.11329.0714.87Fatty acid-binding protein 9
15LOC102401361102401361XP_025146712.13466.939.13LOW QUALITY PROTEIN: uncharacterized protein LOC102401361
16MP2P102409775XP_006068109.11329.6714.95Myelin P2 protein
17RBP1102389538XP_006055984.11354.8815.69Retinol-binding protein 1
18RBP2.1102401674XP_006070028.11345.7615.7Retinol-binding protein 2
19RBP2.2102401674XP_006070029.11345.7615.7Retinol-binding protein 2
20RBP5-X1.1102393311XP_025138922.11865.7921.5Retinol-binding protein 5 isoform X1
21RBP5-X1.2102393311XP_025138923.11865.7921.5Retinol-binding protein 5 isoform X1
22RBP5-X2102393311XP_025138924.11585.6818.72Retinol-binding protein 5 isoform X2
23RBP5-X3102393311XP_025138925.11475.2917.44Retinol-binding protein 5 isoform X3
24RBP5-X4102393311XP_025138926.11355.9315.96Retinol-binding protein 5 isoform X4
25RBP5-X5102393311XP_025138927.11335.4615.82Retinol-binding protein 5 isoform X5
26RBP7102408740XP_006076768.113c46.8215.54LOW QUALITY PROTEIN: retinoid-binding protein 7
Table 2. Ka/Ks ratios estimated in the present work for duplicated FABP genes in buffalo.
Table 2. Ka/Ks ratios estimated in the present work for duplicated FABP genes in buffalo.
Seq_1Seq_2KaKsKa_KsDivergence Time (Mya)
FABP12FABP40.2712328851.2934543860.20969652151.328
FABP4FABP90.2613527772.492795340.10484325598.92
LOC102401361MP2P2.9252823751.61261781.81399608563.993
LOC102401361RBP21.993548971.4745846421.35193932858.515
MP2PFABP90.2343702412.2487702970.10422151289.237
Note. Ka: nonsynonymous substitution rate; Ks: synonymous substitution rate; Ka_Ks: nonsynonymous/synonymous; Mya: million years ago.
Table 3. Ka/Ks ratios estimated in the present work for buffalo and cattle orthologous FABP gene pairs.
Table 3. Ka/Ks ratios estimated in the present work for buffalo and cattle orthologous FABP gene pairs.
BuffaloCattleKaKsKa_Ks
LOC102401361RBP1NaNNaNNaN
FABP3FABP30.1076716460.2279996760.472244734
FABP3FABP70.310656251.6792322790.184998975
RBP5RBP50.085004550.1796021490.473293616
RBP7RBP70.0062565530.03798280.164720704
CRABP2CRABP20.0031274480.0440089570.071063894
RBP1RBP10.0413349630.0909220010.454620033
FABP2FABP2NaNNaNNaN
FABP6FABP60.0235643510.037390920.630215853
FABP7FABP700.0492990570
FABP1FABP10.0099558980.0533558250.186594405
FABP5FABP500.0578321060
FABP12PMP20.2577839071.8491164230.139409236
CRABP1CRABP1NaNNaNNaN
Note. Ka: nonsynonymous substitution rate; Ks: synonymous substitution rate; Ka_Ks: nonsynonymous/synonymous
Table 4. Characterization of the SNPs within FABP genes that were genotyped in the Buffalo population of 489 individuals studied in the present work and that passed the quality control filtering.
Table 4. Characterization of the SNPs within FABP genes that were genotyped in the Buffalo population of 489 individuals studied in the present work and that passed the quality control filtering.
GeneProbe Set IDLocationGenotypeNumberFrequencyAllelesRateObserved HePredicted HeHWE (p-Value)PIC
LOC102401361AX-85097756IntroAA7215.7%A0.389 0.4640.4750.6770.362
AG21446.5%G0.611
GG17437.8%
AX-85049047IntroTT4610.0%T0.307 0.4130.4250.6160.335
TC19141.4%C0.693
CC22448.6%
AX-85116471IntroAA12727.5%A0.508 0.4650.5000.1570.375
AG21546.5%G0.492
GG12026.0%
AX-85106417IntroAA27559.5%A0.779 0.3680.3440.1770.285
AG17036.8%G0.221
GG173.7%
AX-85072673IntroAA194.1%A0.183 0.2820.2980.3060.254
AG13028.3%G0.817
GG31167.6%
RBP2AX-85109932IntroTT25855.8%T0.753 0.3900.3720.3770.303
TC18039.0%C0.247
CC245.2%
RBP5AX-85111933IntroTT245.2%T0.231 0.3570.3551.0000.292
TC16535.7%C0.769
CC27359.1%
Frequency (%), the frequency of individuals with each haplotype among the population.
Table 5. SNP association analysis for six milk production traits in buffalo.
Table 5. SNP association analysis for six milk production traits in buffalo.
GeneProbe Set ID Traits (LSM ± SE)
GenotypePM270/kgMY270/kgPY270/kgPP270/%FY270/kgFP270/%
LOC102401361AX-85097756AA15.8 ± 0.43041 ± 82249 ± 78.21 ± 0.13139 ± 44.59 ± 0.04
AG15.7 ± 0.33040 ± 71247 ± 68.17 ± 0.12138 ± 34.57 ± 0.04
GG15.6 ± 0.43062 ± 74248 ± 68.16 ± 0.12139 ± 34.58 ± 0.04
P-value0.6660.8620.8960.8670.9010.766
AX-85049047CC15.6 ± 0.33048 ± 71247 ± 68.16 ± 0.12139 ± 34.58 ± 0.04
TC15.8 ± 0.33047 ± 73249 ± 68.2 ± 0.12138 ± 34.56 ± 0.04
TT15.8 ± 0.43032 ± 88249 ± 78.15 ± 0.14140 ± 44.60 ± 0.04
P-value0.6460.970.7710.8520.8520.343
AX-85116471AA15.4 ± 0.43035 ± 75245 ± 68.16 ± 0.12138 ± 34.59 ± 0.04
AG15.8 ± 0.33053 ± 72249 ± 68.17 ± 0.12139 ± 34.56 ± 0.04
GG15.8 ± 0.43045 ± 77249 ± 68.19 ± 0.13139 ± 34.58 ± 0.04
P-value0.1730.9190.570.940.9550.387
AX-85106417AA15.6 ± 0.33047 ± 70 a247 ± 68.16 ± 0.12138 ± 34.57 ± 0.04
AG15.9 ± 0.33074 ± 72 a250 ± 68.18 ± 0.12140 ± 34.57 ± 0.04
GG15.2 ± 0.62700 ± 126 b230 ± 108.23 ± 0.20130 ± 54.61 ± 0.06
P-value0.1970.0040.0630.9130.1250.804
RBP2AX-85072673AA15.6 ± 0.52997 ± 116245 ± 98.20 ± 0.19135 ± 54.54 ± 0.06
AG15.9 ± 0.43067 ± 76249 ± 68.17 ± 0.13139 ± 34.57 ± 0.04
GG15.6 ± 0.33042 ± 70247 ± 68.17 ± 0.12139 ± 34.58 ± 0.03
P-value0.5620.7430.8040.9880.7250.606
RBP5AX-85111933CC15.7 ± 0.33034 ± 70247 ± 68.17 ± 0.12138 ± 34.57 ± 0.03
TC15.7 ± 0.33070 ± 74250 ± 68.18 ± 0.12140 ± 34.58 ± 0.04
TT15.1 ± 0.53015 ± 111244 ± 98.15 ± 0.18137 ± 54.55 ± 0.06
P-value0.2970.6270.440.9820.5080.79
LSM ± SE represents the least square means ± standard error; PM270, 270-day peak milk yield; MY270, 270-day total milk yield; FY270, 270-day fat yield; FP270, 270-day fat percentage; PY270, 270-day protein yield, and PP270, 270-day protein percentage. Values with different superscripts differ significantly for each SNP each trait at Bonferroni corrected p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ye, T.; Shaukat, A.; Yang, L.; Chen, C.; Zhou, Y.; Yang, L. Evolutionary and Association Analysis of Buffalo FABP Family Genes Reveal Their Potential Role in Milk Performance. Genes 2022, 13, 600. https://doi.org/10.3390/genes13040600

AMA Style

Ye T, Shaukat A, Yang L, Chen C, Zhou Y, Yang L. Evolutionary and Association Analysis of Buffalo FABP Family Genes Reveal Their Potential Role in Milk Performance. Genes. 2022; 13(4):600. https://doi.org/10.3390/genes13040600

Chicago/Turabian Style

Ye, Tingzhu, Aftab Shaukat, Lv Yang, Chao Chen, Yang Zhou, and Liguo Yang. 2022. "Evolutionary and Association Analysis of Buffalo FABP Family Genes Reveal Their Potential Role in Milk Performance" Genes 13, no. 4: 600. https://doi.org/10.3390/genes13040600

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