BioMed Research International

BioMed Research International / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 8064613 | https://doi.org/10.1155/2017/8064613

Xueqin Liu, Pu Zhang, Gongying Zhang, Sichen Li, Long Zhang, Zhongxian Xu, Tianyuan Ma, Diyan Li, "A Comparison of Genetic Diversity of COX-III Gene in Lowland Chickens and Tibetan Chickens", BioMed Research International, vol. 2017, Article ID 8064613, 13 pages, 2017. https://doi.org/10.1155/2017/8064613

A Comparison of Genetic Diversity of COX-III Gene in Lowland Chickens and Tibetan Chickens

Academic Editor: Daniele Corsaro
Received21 Jan 2017
Revised03 May 2017
Accepted25 May 2017
Published03 Jul 2017

Abstract

To obtain a full understanding of the genetic diversity of the cytochrome oxidase III gene (COX-III) and its association with high altitude adaptation in Tibetan chickens, we sequenced COX-III in 12 chicken populations (155 Tibetan chickens and 145 other domestic chickens). We identified a total of 11 single nucleotide polymorphisms (SNPs) and 12 haplotypes (Ha1–Ha12). Low genetic diversity (haplotype diversity = 0.531 ± 0.087, nucleotide diversity = 0.00125) was detected for COX-III, and haplotype diversity of Tibetan chicken populations (0.750 ± 0.018) was markedly higher than lowland chicken populations (0.570 ± 0.028). Obvious genetic differentiation (nucleotide divergence = 0.092~0.339) and conspicuous gene communication (gene flow = 0.33~32.22) among 12 populations suggested that Tianfu black-bone fowl (white feather) was possibly introduced from Tibetan chicken. SNP m.10587 T>C affects the specific functions of the COX enzyme. Haplotype Ha3 was found in Tibetan chickens, and SNP m.10115G>A caused an amino acid substitution (Val62Ile) associated with phospholipid binding, while mutations m.10017C>A and m.10555G>A and the previously reported SNP m.10065T>C reduced the hydropathy index to some extent. Together, this indicates that the mitochondrial membrane is more hydrophobic in Tibetan chickens.

1. Introduction

Domestic chickens fulfill various roles ranging from food and entertainment to religion and ornamentation [1]. Tibetan chicken is a widely distributed aboriginal chicken breed found at altitudes ranging from 2200 to 4100 m; it has adapted well to high altitudes after over the years of living on the plateau [2]. Generally, long-term exposure to hypoxia in animals reduces metabolic activity, retards development, and increases embryo mortality [3]. However, Tibetan chicken has developed an adaptive mechanism to hypoxia, demonstrated by its increased hatchability and survival rate compared with lowland chicken breeds in high altitude areas of Tibet. With its low weight, small size, and strong chest and legs, the appearance and behaviour of the Tibetan chicken resemble those of the Cochin-Chinese red jungle fowl (Gallus gallus gallus), making it particularly good at flying and foraging on the plateau alpine region [25]. The breed is also important to the resources used to expand the industry in cold areas of high altitude in China.

Avian species living at high altitudes are characterized by the high oxygen affinity of their haemoglobin [6, 7]. Tibetan birds can improve their physiological performance by enhancing their oxygen transport capacity, which has yielded important insights into the genetic basis of adaptation involving haemoglobin as an oxygen carrier [8, 9]. Taking into account the importance of using oxygen more efficiently under hypoxic conditions, information mining for the cellular respiratory chain is a valuable way of understanding the hypoxia response and adaptation mechanisms.

Mitochondrial DNA (mtDNA) sequences are widely used in molecular evolutionary studies. These sequences are useful for estimating times of species and population divergences, comparisons of relative rates of evolution, and phylogenetic inferences within and between vertebrate species [10]. The complete sequence of the chicken mtDNA is 16,775 base pairs and contains 13 protein coding genes, two rRNA genes, and 22 tRNA genes [11]. Cytochrome c oxidase (COX), the terminal enzyme of the mitochondrial respiratory chain, contains 14 protein subunits in mammals, of which three (COX-I, COX-II, and COX-III) are synthesized in the mitochondria [12]. Mitochondrial cytochrome c oxidase subunit I (COX-I) is the central catalytic subunit of cytochrome c oxidase (complex IV), and COX-I gene is used as a standard marker for DNA bar coding to enable species identification in animals [3, 13]. The COX-III protein is an important element in regulating the efficiency of proton translocation in cytochrome oxidase over several turnovers [14]. The COX-III gene of bar-headed geese contained a nonsynonymous substitution (Trp-116→Arg) that resulted in a major functional change of amino acid class. This mutation was predicted by structural modeling to alter the interaction between COX-III and COX-I, which contributed to adaptation in mitochondrial enzyme kinetics and O2 transport capacity and may finally contribute to the exceptional ability of bar-headed geese to fly at extreme heights [15].

COX has been investigated in several biological studies [3, 6, 16, 17], but studies of COX-III have rarely been reported in chickens. However, COX-III is an important component of the respiratory chain and is very conserved among species. Adaptive changes in COX activity can alter the ATP supply derived from oxidative phosphorylation during hypoxia [9]. The quaternary structure formed by different protein subunits is stabilized mainly through hydrophobic interactions in spite of hydrogen bonding and the van der Waals force is also important. We hypothesized that the stability of the COX holoenzyme three-dimensional structure would increase in line with increases in COX-III protein hydrophobicity. Therefore, in the present study, based on the assumption that cytochrome c oxidase activity is more stable in Tibetan chickens, we analyzed COX-III SNPs in 12 chicken populations (five lowland and seven highland populations) to better understand the COX-III genetic diversity and to determine the contribution of specific SNPs to high altitude adaptations in Tibetan chickens.

2. Materials and Methods

In all experimental populations, 5 populations (Muchuan, Emei, Jiuyuan, Black-Tianfu, and White-Tianfu) belonged to lowland chickens; the other 7 populations (Haiyan, Doilungdêqên, Ganzi, Nyingchi, Diqing, Shannan, and Shigatse) which belonged to Tibetan chickens were collected (Table 1). Blood samples were collected from the wing vein. No bird was slaughtered or unexpectedly injured during sampling. The protocol was approved by the Committee on the Care and Use of Laboratory Animals of the State-Level Animal Experimental Teaching Demonstration Center of Sichuan Agricultural University (Approval ID: Decree number S20160906).


TypeBreedsPopulationAltitudeSample ()Haplotype distributions (number of birds)

Lowland chickenMuchuan black-boneMuchuan400 m26Ha1 (22) H2 (1) Ha6 (2) Ha11 (1)
Emei black fowlEmei500 m18Ha1 (4) Ha5 (10) Ha6 (4)
Jiuyuan black fowlJiuyuan800 m43Ha1 (9) Ha5 (33) Ha9 (1)
Tianfu black-bone fowl (black feather)Black-Tianfu800 m44Ha1 (38) Ha5 (2) Ha6 (3) Ha12 (1)
Tianfu black-bone fowl (white feather)White-Tianfu800 m14Ha1 (8) Ha2 (1) Ha5 (5)

Highland chickenTibetan fowlHaiyan3100 m49Ha1 (20) Ha2 (10) Ha5 (12) Ha6 (5) Ha8 (1) Ha10 (1)
Tibetan fowlDoilungdêqên3658 m33Ha1 (9) Ha2 (13) Ha4 (1) Ha5 (7) Ha6 (2) Ha7 (1)
Tibetan fowlGanzi3390 m10Ha1 (2) Ha3 (6) Ha5 (1) Ha6 (1)
Tibetan fowlNyingchi3100 m11Ha1 (1) Ha2 (6) Ha6 (4)
Tibetan fowlDiqing3280 m15Ha1 (6) Ha3 (9)
Tibetan fowlShannan3700 m29Ha1 (7) Ha2 (21) Ha5 (1)
Tibetan fowlShigatse3836 m8Ha1 (4) Ha2 (4)

2.1. DNA Extraction, Amplification, and Sequencing

We extracted mtDNA by salt extraction method [18]. PCR used the known primer pairs F9797: 5′-ACCAATAATACCATCAATCTCC-3′ and R10830: 5′ CGCTTAGTAGAAAGGATAGTGAG-3′ [19, 20]. PCR amplification was performed in a 50 µl volume with 100–150 ng of genomic DNA, 25 mM MgCl2, 2.5 mM of dNTP mixture, 2 mM each primer, 5 µl of 10x buffer, and 1.25 U LA Taq polymerase (Takara, Dalian, China) under the following conditions: denaturation at 94°C for 5 min, then 35 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 60 s, followed by a final extension at 72°C for 7 min [19]. PCR products were verified on 1.5% agarose gels, and specific bands were purified using the TIANgel Midi Purification Kit (Tiangen Biotech, Beijing, China). Purified PCR products were sequenced in both directions using the Big Dye Terminator v. 3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA) on the ABI Prism 3100 DNA sequencer (Applied Biosystems) according to the manufacturer’s instructions.

2.2. Sequence Data Analysis and Statistical Analysis

Raw sequences were aligned and edited by DNAstar software (DNAstar Inc. Madison, WI, USA). We exported all sequences as an aligned FASTA file. Sequence variations were identified using MEGA 6.0 software [21]. Standard population genetics statistics, including haplotypes and number of haplotypes, haplotype diversity within each group (Hd), nucleotide diversity (Pi), nucleotide divergence (Dxy), net genetic distance (Da), coefficient of differentiation (Gst), and gene flow (Nm), Tajima’s value neutral test, were defined using DnaSP V5 software [22], whereas median joining network analysis was performed using program network 4.611 (http://www.fluxus-engineering.com/sharenet.htm). The complete mitochondrial genome sequence of the red jungle fowl was used as the reference sequence (GenBank accession number: NC_001323).

Analysis of molecular variance (AMOVA) was estimated using Arlequin 3.0 software [23]. Bayesian inference was performed as previously described [24]. The optimal model for each data set was estimated by the program Modeltest 3.7 [25]. The program BEAUti v1.5.3 (distributed with BEAST) was used to create the input file to run in BEAST (http://beast.bio.ed.ac.uk/). Samples from the posterior were summarized on the maximum clade credibility tree using the program TreeAnnotator v1.4.8 (distributed with BEAST) and visualized using the program FigTree v1.3.1 (http://tree.bio.ed.ac.uk/software/figtree/). Statistical differences in COX-III haplotype frequencies between Tibetan chickens and lowland chickens were analyzed using Fisher’s exact test; odds ratios (OR) and 95% confidence intervals (95% CIs) were also calculated. value < 0.05 was taken into account as statistical significance. The bioinformatics platform MitoTool (http://www.mitotool.org/) was used to analyze haplotype distribution frequencies between Tibetan chickens and lowland chickens [26]. were taken to be statistically significant.

Protter [27] was used to annotate and predict protein sequence features of the COX-III protein (http://wlab.ethz.ch/protter/start/). Hydropathy plot of the COX-III protein was predicted using the Tmpred Program [28] (http://www.ch.embnet.org/software/TMPRED_form.html).

3. Results

3.1. Nucleotide Diversity of COX-III

The length of 300 COX-III sequences was truncated into 784 bp (GenBank accession numbers: NC_001323; no insertion/deletions were detected). We calculated the overall base composition of COX-III from the 12 chicken populations using MEGA 6.0. Cytosine (C) was shown to be the rarest nucleotide (16%) with guanine (G) to the most common (32%). The A + T% was around half of COX-III (52%).

A total of 11 SNPs (including six singleton sites and five parsimony-informative sites (with no insertions/deletions)) accounted for 1.403% of the total 784 bp COX-III sequence from 300 individuals. Table 2 summarized the number of haplotypes, Hd, Pi, and other information within each group. The number of variable sites in each population varies from 1 in Diqing and Shigatse to 7 in Haiyan, and the sequences from Doilungdêqên, Ganzi, Nyingchi, Shannan, were 5, 5, 3, and 2, respectively. The highest haplotype diversity was found in Haiyan and Doilungdêqên. While the average number of nucleotide differences in Emei from lowland chickens was greater than others. Taken together, a total of 8 haplotypes were identified in lowland chickens, and the overall haplotype diversity, nucleotide diversity, and average nucleotide differences were 0.570 ± 0.028, 0.00155, and 1.216, respectively. In Tibetan chickens, a total of 9 haplotypes were identified, and the overall haplotype diversity, nucleotide diversity, and average nucleotide differences were 0.750 ± 0.018, 0.0016, and 1.250, respectively. The result showed that the genetic diversity of the Tibetan chicken was noticeably higher than that of lowland chickens. After Tajima’s value neutral test, all values were greater than 0.1; therefore, the 12 populations belong to neutral mutations.


PopulationsNumber of variable sitesNumber of parsimony-informative sitesNumber of haplotypesHaplotype diversity (Hd)Nucleotide diversity (Pi)Average number of nucleotide differences ()Tajima’s

Lowland chickenMuchuan5340.000920.720−1.27664
Emei4430.002091.6211.18115
Jiuyuan3230.000990.7770.25955
Black-Tianfu5440.000790.559−1.27905
White-Tianfu2230.001301.0221.69598
All lowland chickens7480.001551.216−0.07994

Highland chickenHaiyan7460.001991.2930.37142
Doilungdêqên5460.001471.152−0.17397
Ganzi5240.001901.489−0.63193
Nyingchi3230.001531.2000.58729
Diqing1120.000660.5141.37595
Shannan2130.000570.448−0.24788
Shigatse1120.000730.5711.44416
All highland chickens8590.00161.250−0.27902

3.2. Nucleotide Divergence and Net Genetic Distance among Populations

Nucleotide polymorphism among populations can be represented by nucleotide divergence (Dxy) and net genetic distance in nucleotides (Da). Within the 12 chicken populations, the average Dxy was 0.062% (range 0.092%–0.339%) and the average Da was 0.196% (range 0%–0.216%) (Table 3). The largest Da (0.216%) was found between Diqing and Jiuyuan black chickens (between Tibetan chickens and lowland chickens), and the smallest Da (0%) was observed between Haiyan and Doilungdêqên chickens; the smallest Dxy (0.092%) was observed between Black-Tianfu and Muchuan black chickens (within lowland chickens) and the largest Dxy (0.339%) was found between Diqing and Emei chickens (between Tibetan chickens and lowland chickens).


MuchuanEmeiJiuyuanBlack-TianfuWhite-TianfuHaiyanDoilungdêqênGanziNyingchiDiqingShannanShigatse

Muchuan0.2580.2350.0920.1440.1740.1810.1980.2320.140.1340.104
Emei0.0920.1840.2580.2140.2260.2110.3380.2390.3390.1930.214
Jiuyuan0.1280.0130.2330.1620.1950.1810.3140.2570.3080.1620.184
Black-Tianfu−0.0020.0990.1350.1390.1710.180.1910.2350.1310.1620.101
White-Tianfu0.0200.0250.0340.0230.1690.1660.2340.2340.1980.1290.122
Haiyan0.0280.0160.0440.0330.0010.1770.2630.2150.2370.1380.141
Doilungdêqên0.0470.0140.0430.0550.01100.2680.1900.250.1180.132
Ganzi0.040.1160.1520.0420.0550.0620.0800.3180.140.2260.2
Nyingchi0.0950.0380.1160.1070.0760.0350.0220.1260.3190.150.175
Diqing0.0520.1870.2160.0510.0880.1050.1310.0030.1970.1990.157
Shannan0.0510.0450.0740.0740.0240.0110.0040.0880.0320.130.076
Shigatse0.0120.0580.0880.0170.0090.0050.0090.0530.0490.0790.004

Note. Lower diagonal was nucleotide divergence Da; upper diagonal was net genetic distance Dxy; all the values were enlarged 100 times.

To pinpoint the most and least closely related populations, we used net genetic distance in nucleotides (Da) to demonstrate that Haiyan and Doilungdêqên are most similar in highland chickens (Da  =  0%, Dxy  =  0.177%) and that Muchuan black-bone fowl and Tianfu black-bone fowl (black feather) are most similar in lowland chickens (Da  =  0.002%, Dxy  =  0.092%). Conversely, Tibetan chickens (Diqing) appear more distantly related to Jiuyuan black fowl (Da  =  0.216%, Dxy  =  0.308%) and Emei black fowl (Da  =  0.187%, Dxy  =  0.339%).

These findings are conclusive of definite genetic differentiation between different groups of chickens, with the strongest differentiation seen between Tibetan chickens and lowland chickens. Net genetic distances in nucleotides (Da) are substantially consistent with the outcome of nucleotide differences (Dxy).

3.3. Coefficient of Differentiation (Gst) and Gene Flow (Nm) between Populations

The coefficient of differentiation (Gst) can reveal the extent of gene flow and genetic drift to some extent, while gene flow can uncover possible gene infiltration among populations. In lowland and Tibetan chicken populations, distinct gene exchange (Gst  =  0.01, Nm  =  25.21; Gst  =  0.05, Nm  =  4.9) was detected between White-Tianfu (lowland chickens) and Haiyan and Doilungdêqên (Tibetan chickens). However, at the same altitude, obvious genetic differentiation was found in lowland geographical populations () while conspicuous gene communication (Gst  =  0.01  Nm  =  29.34) was only detected in Muchuan and Black-Tianfu chickens. Obvious genetic differentiation appears to have occurred in Tibetan chickens (Table 4). This indicates that differences in varieties caused genetic differentiation and that the introduction of varieties of different regions led to gene exchange.


MuchuanEmeiJiuyuanBlack-TianfuWhite-TianfuHaiyanDoilungdêqênGanziNyingchiDiqingShannanShigatse

Muchuan0.650.37−29.342.632.461.170.660.540.760.441.75
Emei0.289.230.664.395.883.051.41.190.850.671.2
Jiuyuan0.400.030.332.141.661.230.90.660.560.40.86
Black-Tianfu−0.010.270.432.891.880.930.710.540.760.381.64
White-Tianfu0.090.050.110.0825.214.91.310.981.230.932.94
Haiyan0.090.040.130.120.0132.222.813.632.132.047.7
Doilungdêqên0.180.080.170.210.050.012.076.551.465.4513.33
Ganzi0.280.150.220.260.160.080.111.0116.950.761.06
Nyingchi0.320.170.280.320.200.060.040.200.724.573.46
Diqing0.250.230.310.250.170.110.150.020.260.591.12
Shannan0.360.270.390.400.210.110.040.250.520.3010.16
Shigatse0.130.170.230.130.080.030.020.190.070.180.02

Note. Lower diagonal was nucleotide divergence Gst, upper diagonal was net genetic distance Nm; meant significant difference between lowland chicken and Tibetan chicken breed.
3.4. Analysis of Molecular Variance

AMOVA showed that the percentage of variation within populations (70.43%) was greater than that between populations (29.57%). Fst value was 0.2957 () which implied that the genetic divergence within populations was significant. The results indicate that the twelve geographic populations do not produce largely genetic differentiation, while the genetic diversity in COX-III gene mainly comes from within the populations (Table 5).


Source of variationSum of squaresdfVariance componentsPercentageFixation index

Among populations56.807110.19 (Va)29.57
Within populations132.7312880.46 (Vb)70.43

Note. Fixation index: 0.2957; .
3.5. Sequence Variations in COX-III Gene

We detected seven synonymous mutations (m.10081A>G, m.10162G>A, m.10270G>A, m.10336A>G, m.10369G>A, m.10587 T>C, and m.10809C>T) and four nonsynonymous substitutions (m.10017C>A, m.10112 G>A, m.10115G>A, and m.10555G>A) in COX-III. Table 6 shows the observed allele frequencies in each polymorphic site between Tibetan chicken and lowland chicken breeds. After using Pearson chi-square test, we found that three SNPs compared with lowland chickens in COX-III gene (m.10115G>A, m.10270G>A, and 10587 T>C) were significantly different with the allele frequency of 9.7%, 91.0%, and 86.5% () in Tibetan chicken, respectively. Similarly, in lowland chicken, the SNP m.10081A>G was significantly different from Tibetan chicken with the allele frequency of 57.2%. The four SNPs (m.10017C>A and m.10369G>A, m.10555G>A, and m.10809C>T) were only distributed in Tibetan chicken breeds with the allele frequency of only 0.6% and 7.7%, while the other three SNPs (m.10112 G>A, m.10162G>A, and m.10336A>G) were only distributed in lowland chicken with the allele frequency of 0.7%. We consider that the seven SNPs showed nonsignificant difference between Tibetan chicken and lowland chicken breeds ().


SNP sitesAllele distribution value in Pearson chi-square test
AlleleTCLC

10017C154 (99.4%)145 (100%)0.33263
A1 (0.6%)0 (0%)
10081A90 (58.1%)62 (42.8%)0.00805
G65 (41.9%)83 (57.2%)
10112G155 (100%)144 (99.3%)0.30037
A0 (0%)1 (0.7%)
10115G140 (90.3%)145 (100%)0.00012
A15 (9.7%)0 (0%)
10162G155 (100%)144 (99.3%)0.30037
A0 (0%)1 (0.7%)
10270G14 (9%)136 (93.8%)
A141 (91.0%)9 (6.2%)
10336A155 (100%)144 (99.3%)0.30037
G0 (0%)1 (0.7%)
10369G154 (99.4%)145 (100%)0.33263
A1 (0.6%)0 (0%)
10555G154 (99.4%)145 (100%)0.33263
A1 (0.6%)0 (0%)
10587T21 (13.5%)94 (64.8%)
C134 (86.5%)51 (35.2%)
10809C143 (92.3%)136 (93.8%)0.60255
T12 (7.7%)9 (6.2%)

Note: relative content in parentheses meant the number of birds in corresponding allele in the specific polymorphic site; meant significant difference between lowland chicken and Tibetan chicken breeds; TC was the abbreviations for Tibetan chicken; LC was the abbreviations for lowland chicken.
3.6. Median Joining Network of Haplotypes and Phylogenetic Analysis

We identified 12 haplotypes (Ha1–Ha12) in 300 chickens from the 12 different populations (Table 7). The median joining network was constructed using the 12 haplotypes. Three clusters (A, B, and C) were clearly defined from the network with substantial mutation distances visible between the clusters (Figure 1). Ha1, Ha3, Ha8, Ha9, and Ha12 were restricted to cluster A: Ha1 was the dominant haplotype, present in 87.8% of all individuals in cluster A (130/148 = 87.8%); Ha2, Ha5, Ha10, and Ha11 were restricted to cluster B: Ha5 was the dominant haplotype, present in 55.03% of all individuals in cluster B (71/129 = 55.03%); Ha4, Ha6, and Ha7 were restricted to cluster C: Ha6 was the dominant haplotype, present in 91.3% of all individuals in cluster B (21/23 = 99.3%). The ancestral haplotype was mainly distributed in the center of the median joining network, with derivative haplotypes spreading outwards from it. This indicates that Ha1, Ha2, and Ha4 are the earliest haplotypes.


Haplotype reference sequencePosition of variable sites in COX-III geneNumber of birds in highland chickens (%)Number of birds in lowland chickens (%)Number of birds
1001710081101121011510162102701033610369105551058710809
CAGGGGAGGTC

Ha1G49 (31.6%)81 (55.9%)130
Ha254 (34.8%)2 (1.4%)56
Ha3GA15 (9.7%)15
Ha41 (0.6%)1
Ha521 (13.5%)50 (34.5%)71
Ha612 (7.7%)9 (6.2%)21
Ha71 (0.6%)1
Ha8GAC1 (0.6%)1
Ha9GG1 (0.7%)1
Ha10A1 (0.6%)1
Ha11A1 (0.7%)1
Ha12G1 (0.7%)1
AA SubstSer29TyrVal61MetVal62IleAla162Thr

Note. Black dot () indicates the same base as the reference, Ha1–Ha12 are haplotypes, and dashes (—) represented the absence of certain haplotype in the population amino acid substitutions (AA Subt.) that were listed below the nucleotide information and marked in italic; position of variable sites in the COX-III gene of chickens was obtained using the complete mitochondrial genome sequence of RJF as the reference sequence.

We also used the Meleagris gallopavo as an outside group to construct a phylogenetic tree using the Bayesian method. The correlation of the 12 haplotypes is shown in Figure 2, and three clusters can also be seen in the Bayesian tree.

3.7. Association between Haplotype Distribution and Altitude Adaptation

Haplotypes of sample sizes under five were not taken into account here. After reviewing the level of significance by Bonferroni correction, haplotype Ha2 was found to be significantly associated with high altitude adaptation at the 0.05 level ( value, ; OR, 30.375, 95% CI, 7.332–125.837). Haplotypes Ha1 and Ha5 also appear to be significantly associated with lowland adaptation (H1: value, 0.001; OR, 1.918; 95% CI, 1.286–2.860; H5: value, 0.00022, OR, 2.685, CI, 1.567–4.598), whereas haplotype Ha6 does not seem significantly associated with altitude adaptation at the 0.05 level ( value, 0.658; OR, 0.74736, 95% CI, 0.310–1.801) (Table 8).


HaplotypeNumber of Tibetan chickensNumber of lowland chickens valueOR95% CIs

Ha149810.000350.4840.32515–0.72169
Ha254230.3757.332–125.837
Ha31500.00006
Ha4101.000
Ha521500.000080.3490.20373–0.59712
Ha61290.661781.2570.52175–3.030
Ha7101.000
Ha8101.000
Ha9010.483
Ha10101.000
Ha11010.483
Ha12010.483

Total155145

Note: OR: odds ratio. OR < 1 mean haplotype may be negatively associated with high-altitude adaptation; OR = 1 mean haplotype is not associated with high-altitude adaptation; OR > 1, haplotype may be surely associated with high-altitude adaptation. CIs: confidence intervals; meant significant difference between lowland chicken and Tibetan chicken breeds.
3.8. Prediction and Analysis of Secondary Structure Changes in the COX-III Protein

Protein sequencing showed that the Tibetan chicken-specific nonsynonymous COX-III variants m.10017C>A and m.10555G>A are located in the transmembrane helical structure, while SNP m.10115G>A is located in the outer surfaces of the inner mitochondrial membrane (Figure 3). The hydrophobicity of the COX-III protein was not changed by the Val62Ile amino acid change caused by variant m.10115G>A (Figure 4), but this change at a key site of phospholipid binding may affect the combination of phospholipids related to adaptations to a hypoxic environment. Similarly, nonsynonymous substitutions Ser29Tyr (caused by SNP m. 10017C>A) and Ala162Thr (caused by m. 10555G>A) reduced the hydropathy index to some extent.

4. Discussion

Oxygen is one of the critical determinants for normal embryonic and foetal development. In avian embryos, a lack of oxygen causes high foetal mortality, heteroplasia, and cardiovascular dysfunction. The Tibetan chicken breed is native to Tibet and can survive with high hatchability regardless of the negative effects of hypoxia. Animals adapted to high altitudes are characterized by high haemoglobin concentrations and oxygen affinity [6], while lowland chicken breeds suffer polycythemia and ventricular hypertrophy at high altitudes [2, 4, 29].

By analyzing population haplotype diversity and the structure of chicken breeds in southwestern China based on COX-III sequences, we found that the Tibetan chicken has a higher level of haplotype diversity (0.750 ± 0.018) than lowland chickens (0.570 ± 0.028). For the results observed, we think that artificial selection leads to reduced nucleotide diversity in lowland chickens. Although 11 SNP sites were identified in COX-III from 12 geographic populations, the average genetic distance was only 0.196%, revealing a low level of genomic polymorphisms in all populations. This low level of genetic diversity indicates that this gene is functionally important and hence has an evolutionary constraint. The 12 populations belong to neutral mutations (). This suggests that these populations are stable and have not undergone a population expansion in the past few years.

The analysis of the coefficient of differentiation and gene flow between populations suggests that some chicken populations have undergone genetic differentiation at places of equal altitude. This may reflect the introduction of varieties of different regions. For instance, the Tianfu white-bone fowl from a lower altitude has a close population genetic relationship with the Tibetan chicken (Gst  =  0.01, Nm  =  25.21; Gst  =  0.05, Nm  =  4.9), indicating that the Tianfu white-bone fowl was introduced from high altitude area. Moreover, Tibetan chickens from Haiyan and Doilungdêqên regions were most similar (Da  =  0%, Dxy  =  0.177%) suggesting that the genetic distance of the two regions is very close, and gene exchange is very rich (Gst  =  0.01, Nm  =  32.22).

The regulating mechanism in COX-III gene difference is still ambiguous between lowland and highland populations. A previous study detected SNPs in three mitochondrially encoded subunit genes of chicken COX, including only one in COX-III (m.10081A>G) between an expanded sample of 56 Tibetan chickens and 152 lowland birds [20]. Another study identified [19] eight SNPs, of which five (m.10081A>G, m.10115G>A, m.10270G>A, m.10336A>G, and m.10447C>T) showed significant differences between Tibetan chickens and lowland chickens. Only the synonymous mutation m.10081A>G was found to differ between haplotypes H4 and H5, and chickens with the A allele at m.10081A>G had a probability of being over 2.6 times better adapted to hypoxia than those with the G allele indicating that m.10081A>G may be a prerequisite for shaping high altitude adaptation-specific haplotypes. Our focus on COX-III SNPs to explore the different haplotypes detected a novel mutation associated with high altitude adaptations. AMOVA showed that COX-III variation mainly existed within a population. Of the 12 defined haplotypes, the existence of Ha4, Ha7, Ha8, and Ha10 only in highland chickens and Ha9, Ha11, and Ha12 only in lowland chickens indicates different degrees of genetic divergence between Tibetan chickens and lowland chickens. Ha1, Ha2, and Ha4 were found to be the earliest ancestors (Figure 1), while haplotype Ha1 was common to all populations suggesting that it is more stable and capable of adapting to new environmental selection. Moreover, Ha2 had significant relationship with high altitude adaptation ( value, ; OR, 30.375, 95% CI, 7.332–125.837), with the C allele at m.10587 T>C was found to have a probability of being over 30.375 times better adapted to hypoxia than the T allele. We propose that m.10587 T>C affects the functions of the COX enzyme in a similar way to the effect of m.10081A>G [19] on high altitude adaptation of the Tibetan chicken. However, haplotype Ha1 simultaneously contained mutations m.10081A>G and m.10587 T>C and was negatively associated with high altitude adaptation. The function of these mutations warrants need to further research.

Three of the four nonsynonymous mutations identified in the present study (m.10017C>A, m.10115G>A, and m.10555G>A) were peculiar to highland chickens. Mutation m.10115G>A (with the allele frequency of 9.7%, ), shared by populations Ganzhi and Diqing, caused the amino acid mutation Val62Ile (Figure 1). A previous study reported SNP m.10115G>A as an uncommon missense mutation in the Tibetan chicken mtDNA genome [30], while another study found that more than one-third of Tibetan chickens (44/125, 35.2%) harboured this mutation, suggesting that it might be associated with high altitude adaptation [19]. In the current study, we found that mutation m.10115G>A is a key site for phospholipid binding, suggesting that it impacts on the function of the mitochondrial membrane in Tibetan chicken. Whereas the other two nonsynonymous mutations (m.10017C>A m.10555G>A) and the reported SNP m.10065T>C [19] both reduced hydropathy index to a certain extent.

5. Conclusions

In conclusion, we found largest genetic differentiation between Tibetan and lowland breeds and identified haplotype Ha2 is associated with Tibetan chicken populations. The possible association between increased hydrophobicity/reduced hydrophilic characteristics of the mitochondrial membrane and high altitude adaptation could provide a theoretical reference for poultry genetics. Therefore, our results provide a theoretical basis for future research into fowl breeding.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by China Agricultural Research System (CARS-41), the 12th Five-Year Plan for Breeding Program in Sichuan-Selective Breeding of New Breeds and the Synthetic Strains in Laying Hens (2011NZ0099-7), National Natural Science Foundation of China (31402063), and Sichuan Provincial Department of Science and Technology Program (2015JQO023).

References

  1. Y.-W. Miao, M.-S. Peng, G.-S. Wu et al., “Chicken domestication: An updated perspective based on mitochondrial genomes,” Heredity, vol. 110, no. 3, pp. 277–282, 2013. View at: Publisher Site | Google Scholar
  2. H. G. Bao, C. J. Zhao, J. Y. Li et al., “Sequencing and alignment of mitochondrial genomes of Tibetan chicken and two lowland chicken breeds,” Science in China, Series C: Life Sciences, vol. 51, no. 1, pp. 47–51, 2008. View at: Publisher Site | Google Scholar
  3. X. Zhao, N. Wu, Q. Zhu, U. Gaur, T. Gu, and D. Li, “High-altitude adaptation of Tibetan chicken from MT-COI and ATP-6 perspective,” Mitochondrial DNA, vol. 27, no. 5, pp. 3280–3288, 2016. View at: Publisher Site | Google Scholar
  4. H. Zhang, C. X. Wu, and Y. Z Chamba, “Adaptability to high altitude and NOS activity of lung in Tibetan chicken,” Journal of China Agricultural University, vol. 11, no. 1, pp. 35–38, 2006. View at: Google Scholar
  5. M. Li and C. Zhao, “Study on Tibetan Chicken embryonic adaptability to chronic hypoxia by revealing differential gene expression in heart tissue,” Science in China, Series C: Life Sciences, vol. 52, no. 3, pp. 284–295, 2009. View at: Publisher Site | Google Scholar
  6. C. Liu, L. F. Zhang, M. L. Song, H. G. Bao, C. J. Zhao, and N. Li, “Highly efficient dissociation of oxygen from hemoglobin in Tibetan chicken embryos compared with lowland chicken embryos incubated in hypoxia,” Poultry Science, vol. 88, no. 12, pp. 2689–2694, 2009. View at: Publisher Site | Google Scholar
  7. Y. Liang, Z. Hua, X. Liang, Q. Xu, and G. Lu, “The crystal structure of bar-headed goose hemoglobin in deoxy form: The allosteric mechanism of a hemoglobin species with high oxygen affinity,” Journal of Molecular Biology, vol. 313, no. 1, pp. 123–137, 2001. View at: Publisher Site | Google Scholar
  8. J. Zhang, Z. Hua, J. R. H. Tame, G. Lu, R. Zhang, and X. Gu, “The crystal structure of a high oxygen affinity species of haemoglobin (bar-headed goose haemoglobin in the oxy form),” Journal of Molecular Biology, vol. 255, no. 3, pp. 484–493, 1996. View at: Publisher Site | Google Scholar
  9. Z.-Y. Zhang, B. Chen, D.-J. Zhao, and L. Kang, “Functional modulation of mitochondrial cytochrome c oxidase underlies adaptation to high-altitude hypoxia in a Tibetan migratory locust,” Proceedings of the Royal Society B: Biological Sciences, vol. 280, no. 1756, Article ID 20122646, 2013. View at: Publisher Site | Google Scholar
  10. S. Kumar, “Patterns of nucleotide substitution in mitochondrial protein coding genes of vertebrates,” Genetics, vol. 143, no. 1, pp. 537–548, 1996. View at: Google Scholar
  11. P. Desjardins and R. Morais, “Sequence and gene organization of the chicken mitochondrial genome. A novel gene order in higher vertebrates,” Journal of Molecular Biology, vol. 212, no. 4, pp. 599–634, 1990. View at: Publisher Site | Google Scholar
  12. E. Balsa, R. Marco, E. Perales-Clemente et al., “NDUFA4 is a subunit of complex IV of the mammalian electron transport chain,” Cell Metabolism, vol. 16, no. 3, pp. 378–386, 2012. View at: Publisher Site | Google Scholar
  13. R. P. Kapur, “Emery and Rimoin's Principles and Practice of Medical Genetics, Fourth edition,” Pediatric and Developmental Pathology, vol. 6, no. 2, pp. 204-205, 2003. View at: Publisher Site | Google Scholar
  14. M. BRUNORI, G. ANTONINI, F. MALATESTA, P. SARTI, and M. T. WILSON, “Cytochrome‐c oxidase: Subunit structure and proton pumping,” European Journal of Biochemistry, vol. 169, no. 1, pp. 1–8, 1987. View at: Publisher Site | Google Scholar
  15. G. R. Scott, P. M. Schulte, and S. Egginton, “Molecular evolution of cytochrome C oxidase underlies high-altitude adaptation in the bar-headed goose,” Molecular Biology Evolution, vol. 28, no. 1, pp. 351–363, 2011. View at: Publisher Site | Google Scholar
  16. U.-R. Heinrich, J. Brieger, O. Selivanova et al., “COX-2 expression in the guinea pig cochlea is partly altered by moderate sound exposure,” Neuroscience Letters, vol. 394, no. 2, pp. 121–126, 2006. View at: Publisher Site | Google Scholar
  17. J. L. Park, J. A. Shayman, and L. Shu, “Differential involvement of COX1 and COX2 in the vasculopathy associated with the a-galactosidase A-knockout mouse,” American Journal of Physiology - Heart and Circulatory Physiology, vol. 296, no. 4, pp. H1133–H1140, 2009. View at: Publisher Site | Google Scholar
  18. S. A. Miller, D. D. Dykes, and H. F. Polesky, “A simple salting out procedure for extracting DNA from human nucleated cells,” Nucleic Acids Research, vol. 16, no. 3, article 1215, 1988. View at: Publisher Site | Google Scholar
  19. J. Sun, H. Zhong, S.-Y. Chen, Y.-G. Yao, and Y.-P. Liu, “Association between MT-CO3 haplotypes and high-altitude adaptation in Tibetan chicken,” Gene, vol. 529, no. 1, pp. 131–137, 2013. View at: Publisher Site | Google Scholar
  20. H. Bao, C. Zhao, L. Zhang, J. Li, and C. Wu, “Single-nucleotide polymorphisms of mitochondrially coded subunit genes of cytochrome c oxidase in five chicken breeds,” Mitochondrial DNA, vol. 19, no. 5, pp. 461–464, 2008. View at: Publisher Site | Google Scholar
  21. K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011. View at: Publisher Site | Google Scholar
  22. P. Librado and J. Rozas, “DnaSP v5: a software for comprehensive analysis of DNA polymorphism data,” Bioinformatics, vol. 25, no. 11, pp. 1451-1452, 2009. View at: Publisher Site | Google Scholar
  23. L. Excoffier, G. Laval, and S. Schneider, “Arlequin (version 3.0): an integrated software package for population genetics data analysis,” Evolutionary Bioinformatics Online, vol. 1, no. 4A, pp. L418–L421, 2005. View at: Google Scholar
  24. J. W. Streicher, A. J. Crawford, and C. W. Edwards, “Multilocus molecular phylogenetic analysis of the montane Craugastor podiciferus species complex (Anura: Craugastoridae) in Isthmian Central America,” Molecular Phylogenetics and Evolution, vol. 53, no. 3, pp. 620–630, 2009. View at: Publisher Site | Google Scholar
  25. D. Posada and K. A. Crandall, “MODELTEST: testing the model of DNA substitution,” Bioinformatics, vol. 14, no. 9, pp. 817-818, 1998. View at: Publisher Site | Google Scholar
  26. L. Fan and Y.-G. Yao, “MitoTool: A web server for the analysis and retrieval of human mitochondrial DNA sequence variations,” Mitochondrion, vol. 11, no. 2, pp. 351–356, 2011. View at: Publisher Site | Google Scholar
  27. U. Omasits, C. H. Ahrens, S. Müller, and B. Wollscheid, “Protter: Interactive protein feature visualization and integration with experimental proteomic data,” Bioinformatics, vol. 30, no. 6, pp. 884–886, 2014. View at: Publisher Site | Google Scholar
  28. K. Hofmann, “TMBASE-A database of membrane spanning protein segments,” Biol. Chem. Hoppe-Seyler, vol. 374, no. 1, p. 1, 1993. View at: Publisher Site | Google Scholar
  29. C. Monge, F. Le, and F. León-Velarde, “Le≤n-Velarde, Physiological adaptation to high altitude: Oxygen transport in mamals and birds,” Physiological Reviews, vol. 71, no. 4, pp. 1135–1172, 1991. View at: Google Scholar
  30. X. M. Tong, Y. Liang, and W. Wang, “[Complete sequence and gene organization of the Tibetan chicken mitochondrial genome],” Hereditas, vol. 28, no. 7, pp. 769–777, 2006. View at: Google Scholar

Copyright © 2017 Xueqin Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views971
Downloads408
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.