Sequence Characterization of DSG3 Gene to Know Its Role in High-Altitude Hypoxia Adaptation in the Chinese Cashmere Goat

The Tibetan cashmere goat is one of the main goat breeds used by people living in the plateau. It exhibits the distinct phenotypic characteristics observed in lowland goats, allowing them to adapt to the challenging conditions at high altitudes. It provides an ideal model for understanding the genetic mechanisms underlying high-altitude adaptation and hypoxia-related diseases. Our previous exome sequencing of five Chinese cashmere breeds revealed a candidate gene, DSG3 (Desmoglein 3), responsible for the high-altitude adaptation of the Tibetan goat. However, the whole DSG3 gene (44 kbp) consisting of 16 exons in the goat genome was not entirely covered by the exome sequencing. In this study, we resequenced all the 16 exons of the DSG3 gene in ten Chinese native goat populations. Twenty-seven SNP variants were found between the lowland and highland goat populations. The genetic distance (F_ST) of significant SNPs between the lowland and highland populations ranged from 0.42 to 0.58. By using correlation coefficient analysis, linkage disequilibrium, and haplotype network construction, we found three non-synonymous SNPs (R597E, T595I, and G572S) in exon 5 and two synonymous SNPs in exons 8 and 16 in DSG3. These mutations significantly segregated high- and low-altitude goats in two clusters, indicating the contribution of DSG3 to the high-altitude hypoxia adaptation in the Tibetan goat.

Introduction

Low oxygen or hypoxia is the hardest environmental challenge for humans and animals existing at high altitude. Whole-genome resequencing analyses have been carried out to investigate the genetic basis of the hypoxia adaptation in many domestic animals, such as yak (Qiu et al., 2012), Tibetan pig (Ai et al., 2014), dog (Li et al., 2014), and birds (Cai et al., 2013). Various candidate genes have been identified for high-altitude adaptation, including EPAS1 (endothelial PAS domain protein 1) and HBB (hemoglobin beta) (Gou et al., 2014; Fan et al., 2015; Song et al., 2016). In prior studies, both EPAS1 and HBB revealed six non-synonymous mutations potentially affecting the gene function and influencing high-altitude hypoxic adaptation in dogs (Fan et al., 2015). These studies proved that the domestic animals are a useful animal model to explore high-altitude adaptation.

The domestic goat is found at sea level (30 m) to the high plateau (4700 m). In contrast to its low-altitude counterpart, the Tibetan goat has unique anatomical and physiological characteristics such as higher hemoglobin concentration, a larger heart, and bigger lungs that equip it to live at high altitudes (Li et al., 2004b). Our previous exome-sequencing analysis also reveals that the cardio-vascular system-related genes may play a crucial role in high-altitude adaptation (Song et al., 2016). Among the candidate genes under strong selection, the DSG3 (Desomoglein 3) loci was the most significant. DSG3 is a desmoglein gene located in a cluster on goat chromosome 24 containing 44 kbp that encodes 16 exons. It is expressed in stratifying epithelia, but its precise role is not fully understood in the cardiovascular system (Delva et al., 2009). Here, we extended our previous analyses on large goat populations from different altitudes to identify if the DSG3 candidate gene has a role in high-altitude hypoxia adaptation. By using correlation coefficient analysis, linkage disequilibrium, and haplotype network construction, we explored both the non-synonymous and synonymous mutation sites of the DSG3 gene. These mutations significantly segregated goat populations between high- and low-altitude groups.

Materials and Methods

Populations Sampling

A total of 125 Chinese cashmere goats were used for Sanger sequencing and covered the entire genomic region of DSG3 (44 kbp) consisting of 16 exons. Isolated genomic DNA used from Ritu (RT, 4700 m), Bange (BG, 400 m), Nanjiang (NJ, 1700 m), Inner Magnolia (IM, 1500 m), and Liaoning (LN, 30 m) cashmere goats selected from 330 individuals was analyzed previously for Exome Sequencing by the Illumina Hiseq2000 Platform (Song et al., 2016), which were randomly selected from 4 different locations (Table 1 and Supplementary Table S1 and Supplementary Figure S1). Furthermore, five goat populations including Dulan (DL, 3000 m), Hanshan (HS, 1000 m), Guangfeng (GF, 500 m), Hainan (HN, 120 m), and Changjiangsanjiaozhou (CZ, 50 m) were selected from five different geographical locations in China (Table 1 and Supplementary Table S1 and Supplementary Figure S1). The tissue samples of these goat populations were taken by cutting ear tissue, using an ear-cutting instrument without applying an anesthetic or analgesic agent. After the collection, the tissue samples were preserved in a tube containing 75% alcohol, which was kept in the ice box during transportation to the laboratory. The tissue samples were collected from all the animals undergoing the experiment according to the recommendations and guideline of the Institute of Animal Sciences (IAS, CAAS), with full approval from the Animal Care and Use Committee of Chinese Academy of Agricultural Sciences and the Ministry of Agriculture of the People’s Republic of China.

TABLE 1

TABLE 1. Sample locations and major allele frequency of three non-synonymous variants and two synonymous mutations in DSG3 gene.

Genomic DNA was isolated from ear tissue using a modified commercial Kit protocol (Promega, the Wizard^®). The quality and integrity of the purified DNA from each goat per sample location were determined using a Nanodrop 2000 (Thermo Fisher Scientific, DE) and pooled according to the concentration bases. The validations in the larger goat population were carried out by following sample data and locations. Highland group consisted of Ritu (RT, n = 17) and Bange (BG, n = 17) goats from the Ritu and Bange region of Tibet. The lowland group consisted of Dulan (DL, n = 8) and Nanjiang (NJ, n = 14) cashmere goats from Qinhai, Dulan, and Aksu areas of Xinjiang region, Inner Magnolia (IM, n = 13) from the Earlangshan region of Inner Mongolia, Hanshan (HS, n = 8) from Inner Magnolia city, Guangfeng (GF, n = 15) from Shangrao, Jiangxi., Hainan (HN, n = 8) from Hainan city, Changjiangsanjiaozhou (CZ, n = 13) from Changjiangsanjiaozhou city and Liaoning (LN, n = 15) cashmere goats from Gai Zhou, Liaoning, as shown in Table 1, Supplementary Table S1 and Supplementary Figure S1.

Polymerase Chain Reaction (PCR) and Single-Nucleotide Polymorphisms (SNPs) detection

Sixteen pairs of primers were used for DSG3 gene amplification and their sequences are shown in Table 2. Polymerase chain reaction (PCR) was optimized with a 50 ng DNA template using 2xTaq master mix (BioMed: Beijing Co., Ltd.), supplemented with ddH20 solutions and 10 pmol/l of forward and reverse primers and in a 30 ul mixture. The PCR was performed by initial denaturation for 5 min at 95°C, followed by 40 cycles at 95°C for 30 s, annealing at 60–63°C for 30 s and an extension at 72°C for 45 s, and a final extension at 72°C for 5 min and stored at 4°C. The PCR products were analyzed by running on 1.5% of gel electrophoresis that was made by a mixture of 0.75 g Regular Agarose Biowest (Gene Company) and 50 ml of 1x TAE buffer. The mixture was also added with 5 μl gold-view chemical for viewing the gel on UV transilluminator. The electrophoresis tank (Bio-rad Sub-Cell^® GT) and the PowerPac Basic (Bio-rad^®) were used for running the electrophoresis. Five microliters of the PCR products were mixed with 10X loading dye, which is composed of 0.25% Bromophenol blue, 0.25% Xylene Cyanol FF, 15% Ficoll, and H2O. The mixture was loaded into each running well of 1.5% gel. The electrophoresis was run under 1X TAE buffer, carried out at a 75W constant, 147 V for 40 min. When electrophoresis was finalized, the gels were set on the UV transilluminator. All the purified PCR products were directly sequenced by the Sanger Sequence methods, using the ABI 3730 sequence analyzer (Applied Bio-System, United States) in both the directions, and using the forward primer. The sequences were assembled and analyzed for single-nucleotide polymorphism (SNPs), using the Seq Man program of Laser gene software. The SNPs were detected from Primer 4, Primer 5, Primer 6, Primer 8, Primer 9, Primer 11, Primer 14, and Primer 16 in the DSG3 gene through DNA pooled strategies. Furthermore, these primers were used in large populations on an individual basis to identify polymorphism and the validation of the SNPs.

TABLE 2

TABLE 2. Sequences of the primers.

Statistical Analyses

Pairwise genetic distance (F_ST) was calculated between High-Lowland and High-Midland and Mid-lowland goat populations to find out the genetic divergence as described by Weir and Cockerham (1984). Fisher’s exact test was estimated respectively, to observe a statistical significance. LD-block structure and parameter (D’ and r²) were figured using Haploview 4.2 software (Barrett et al., 2005). To explore the signature of selection and evolutionary change, the six SNPs (four Non-Synonymous and two Synonymous variants) were analyzed for the diversity of haplotype and their frequency was estimated using phase 2.1 (Stephens et al., 2001; Crawford et al., 2004). To obtain a reliable result, we used option – X 10 to increase final runs time and – C option for low- and high-land goat groups to ensure that any similarity to the haplotype is not taken into account. The haplotype was embedded in DNASP Version5 (Librado and Rozas, 2009) for network construction and phylogenetic diversity among the identified haplotypes was inferred through a median-joining network analysis using the method recommended by Bandelt et al. (1999). To detect the functions of three major non-synonymous SNPs in the DSG3 gene, the SIFT software¹ has been used. The protein sequence of Tibetan goats was compared with the reference protein sequence of goats desmoglein-3. The cutoff value in the SIFT program is a tolerance score of ≥ 0.05. The higher the tolerance score, the lower the expected functional impact of a particular amino acid substitution.

Results

Allele Frequency and Genetic Differentiation

In this study, we sampled 128 Chinese indigenous goats belonging to ten populations living across a wide distribution of altitudes. The ten populations ranged from sea level (30 m) to high altitude (4700 m) (Table 1 and Supplementary Table S1). To find out the polymorphism and to measure allele frequencies of the DSG3 gene accurately, we resequenced the entire genomic region of the DSG3 gene. Among the 27 SNPs that were identified in DSG3, five SNPs showed the most significant differentiation in the analysis of allele frequency and global F_ST, including three non-synonymous SNPs, SNP1 (R597E, Chr24: 25794694, exon 5), SNP2 (T595I, Chr24: 25794695, exon 5), and SNP3 (G572S, Chr24: 25794771, exon 5), and two synonymous SNPs, SNP4 (Chr24: 25799255, exon 8), and SNP5 (Chr24: 25817330, exon 16) (Table 3). A striking change in allele frequency with elevation in altitude was observed in goat population. For example, low R597E allele frequencies in DSG3 were found in lowland populations including CZ (0.00, 50m), HN (0.00, 120 m), and HS (0.00, 1000 m), but were elevated in the IM (0.65, 1500 m), DL (0.56, 3000 m), BG (0.67, 4000 m), and reached the highest frequency in the RT (0.97, 4700 m) population (Table 1, Supplementary Table S2, and Figure 1A). The other four SNPs in DSG3 exhibited a similar pattern of change (Table 1, Supplementary Table S2, and Figures 1, 2). All these results indicated that the frequencies of the five variants in the DSG3 locus remain rare in the lowland goat population but increased with the elevation of altitude.

TABLE 3

TABLE 3. Correlation and global F_ST value of five mutation loci in DSG3 gene.

FIGURE 1

FIGURE 1. The pattern of allele frequency of SNP1 (A), SNP2 (C) and SNP3 (E) of DSG3 has continuously changed from sea level at 30 m (LN goat population) to the elevation at 4700 m (RT goat population). Blue, a frequency of the mutant allele, Green, a frequency of a reference allele in each goat population, respectively. The plot of the correlation analysis between major allele frequency of SNP1 (B), SNP2 (D) and SNP3 (F) and sampling locations of 10 cashmere goat populations, r represents the correlation coefficient (r = 0.8529, P < 0.05). Asterisks (^∗∗) represents very significantly.

FIGURE 2

FIGURE 2. The pattern of allele frequency of SNP4 (A) and SNP5 (C) of DSG3 has continuously changed from sea level at 30 m (LN goat population) to the elevation at 4700 m (RT goat population). Blue, a frequency of the mutant allele, Green, a frequency of a reference allele in each goat population, respectively. The plot of the correlation analysis between major allele frequency of SNP4 (B), SNP5 (D) and sampling locations of 10 cashmere goat populations, r represents the correlation coefficient (r = 0.9061, P < 0.05). Asterisks (^∗∗∗) represents extremely significantly.

Using allele frequencies of the three non-synonymous and two synonymous SNPs measured above, we measured the correlation co-efficiency of variant allele frequency per population versus the altitude of sampling location using Graph Pad Prism 5 software. The result of this analysis showed a significant positive linear correlation between the elevated altitude and the frequency of a variant allele in DSG3 (Table 3 and Supplementary Table S3). The linear correlation coefficient between the elevated altitude and the frequencies of variant allele was r = 0.85, P < 0.05 for SNP1 (Figure 1B), r = 0.83, P < 0.05 for SNP2 (Figure 1D), r = 0.84, P < 0.05 for SNP3 (Figure 1F). The linear correlation coefficient was r = 0.91, P < 0.05 for SNP 4 (Figure 2B), and r = 0.90, P < 0.05 for SNP5 (Table 3, Supplementary Table S3, and Figure 2D). The significant positive correlation between the elevated altitudes and the variant allele frequency at the DSG3 locus in ten Chinese native goat populations strongly suggested that the candidate gene DSG3 potentially contributed to the high-altitude adaptation in goats.

To further understand the genetic differentiation that occurred between the lowland and highland goat populations, we first estimated the global F_ST value for all the identified loci. The mean F_ST value was 0.24, indicating the high genetic divergence in these ten Chinese native goat populations, which is mainly attributed to the divergence between the lowland and highland goat populations. The global F_ST of all the 27 loci ranged from 0.0 to 0.61 (Supplementary Table S2). The three most divergent non-synonymous SNPs had the value that ranged from 0.51 (SNP1) and 0.52 (SNP2) to 0.55 (SNP3), whereas the top two synonymous SNPs had a global F_ST of 0.47 (SNP4) and 0.61 (SNP5) (Table 3 and Supplementary Table S2). Subsequently, the pairwise genetic distance (F_ST) was measured between High-Low, High-Mid, and Mid-Low altitude goat populations (Table 4). The F_ST values of non-synonymous variants between high-altitude and low-altitude populations were 0.42 (SNP2, P < 0.001) and 0.45 (SNP1/SNP3, P < 0.001), which reflected a remarkable genetic differentiation (Table 4). Similarly, the F_ST values of the two synonymous substitutes were 0.47 (SNP4, P < 0.001) and 0.58 (SNP5, P < 0.001). These results showed high genetic differentiation between High-Low altitude goat populations compared with the High-Mid and Mid-Low altitude goat populations (Table 4). Taken together, both of global F_ST values and pairwise genetic distance (F_ST) values demonstrate a remarkable population differentiation between the high- and low-altitude goat populations.

TABLE 4

TABLE 4. Pairwise genetic distances (F_ST) & Fisher exact test (p-value) calculated between highland and reference goat populations.

Linkage Disequilibrium (LD) and Haplotype Network Analysis

The LD haploblock pattern was used to analyze the DSG3 loci (Supplementary Figure S2). We constructed LD blocks using the five significant SNPs and calculated D°, r², and an algorithm of the odds (LOD). The variant allele of these five SNPs was tightly linked as a Haplo-Block on chromosome 24 in highland goat populations (RT-BG), with the three non-synonymous SNPs in the main position and showed strong correlation by D°= 1.0 -1.0, r² = 0.63–1.00, and LOD = 5.19–9.82 (Table 5 and Figure 3A). In contrast, the linkage of three non-synonymous SNPs and two synonymous SNPs was loosely found in lowland goat populations (DL-NJ-AB-LN-GF-CZ-HN-HS). However, the LD level of SNPs varied from D = 0.58–1, r² = 0.10–0.94, and LOD = 2.85–31.71 (Table 5 and Figure 3B). The haplotype pattern was analyzed to detect the effect of natural selection on selected genes, including DSG3. Twelve haplotypes were obtained from four non-synonymous (SNP1/R392Q, SNP2/T595I, SNP3/G572S, and T480N) and two from synonymous SNPs (SNP4 and SNP5) from 128 goats from lowland and highland goat populations. The top three haplotypes were detected: HL1, HL2, and L7 with a frequency of 42%, 40%, and 30%, respectively (Figure 3C). Among these haplotypes, HL1 showed a remarkably higher haplotype frequency of a mutant allele (TCCCAA) in highland population than in lowland populations (18.0%). However, HL2 indicated a much lower haplotype frequency of the reference allele (CGTAGG) in the highland (7.0%) than in the lowland population (42%). L7 showed a haplotype frequency of (CGTCGG) in lowland (31%) and indicated a separate haplotype cluster in lowland populations. These results suggested that the haplotype TCCCAA of the DSG3 gene has been selected during the high-altitude adaptation of the Tibetan goat.

TABLE 5

TABLE 5. Pairwise linkage disequilibrium of the non-synonymous and synonymous SNPs of DSG3 in the highland and lowland goat populations.

FIGURE 3

FIGURE 3. (A) Linkage disequilibrium of non-synonymous variants (SNP1, SNP2, and SNP3) and synonymous substitutions (SNP4 and SNP5) in highland goat populations; RT and BG. (B) Linkage disequilibrium of non-synonymous variants (SNP1, SNP2, and SNP3) and synonymous substitutions (SNP4 and SNP5) in lowland goat populations; LN, CZ, HN, GF, HS, IM, NJ, and DL. Strengths of the LD between two SNPs were indicated by the color scheme (a bright red diamond indicates a strong correlation) by the pairwise D° values given as the percentage in each diamond (a bright red diamond indicates D° = 1.00). (C) Haplotype networks of non-synonymous variants (SNP1, SNP2, and SNP3, T480N) and synonymous substitutions (SNP4 and SNP5) among lowland and highland goat populations. Each node represents a haplotype, with the size of the circle proportional to frequency. The length of a branch is proportional to the number of nucleotide differences. Circles are color code according to the population. Blue: highland (RT and BG), Green: lowland (LN, CZ, HN, GF, HS, IM, NJ, and DL), abbreviation: HL, high and low altitude; H, high altitude; L, low altitude.

Prediction of Three Non-synonymous SNPs

An online tool SIFT was used to investigate the potential impact of the non-synonymous substitutions on protein structure and function. SIFT aligns a query sequence and a subject amino acid sequence to determine the effect of an amino acid substitution. SIFT software takes protein FASTA sequences as input and is aligned with PSI-BLAST. In our result analysis, we found three nsSNPs to be predicated, as shown in Table 6, by using SIFT software.

TABLE 6

TABLE 6. Predications of non-synonymous SNPs function by SIFT software.

Discussion

There are many animal species that live on the Tibetan plateau. The plateau consists of 25% of Mainland China, with an average altitude exceeding 4,500 m (Thompson et al., 2000; Wei et al., 2016). Up to now, genomic studies indicated several significant pathways and functional categories including energy metabolism and oxygen transmission, DNA repair, and ATPase production, in response to hypoxia at the Tibetan plateau (Qiu et al., 2012; Ge et al., 2013; Li et al., 2013; Qu et al., 2013; Xiang et al., 2013). Moreover, SNP data have been obtained from population surveys to find the mechanism underlying plateau adaptation (Wei et al., 2016) and have successfully identified candidate genes EGLN1, PPARA, and EPAS1 as transcription factors with a pivotal role in adaptation (Buroker et al., 2012; Xiang et al., 2013; Gou et al., 2014; Wang et al., 2014).

In our previous exome sequencing of 330 cashmere goats, we identified a genomic region containing a desmosome gene DSG3 (Desomoglein3) under strong selection. DSG3 gene is a desomoglein protein-coding gene located in a cluster on goat chromosome 24, consisting of 16 exons.

In this study, the entire region of (44 kbp) DSG3 was sequenced. Subsequently, the single-nucleotide polymorphism was identified in a large population of indigenous Chinese goats to determine the genetic difference between the lowland and highland goat populations. Our results showed polymorphism in DSG3 (a) segregated indigenous Chinese goat population in two clusters; (b) a significant positive linear correlation between the elevated altitude and the frequency of a mutant allele; (c) non-synonymous candidate mutations indicating the contribution of gene in high-altitude adaptation of Tibetan goats.

Moreover, our previous studies suggested that the positive directional selections of genes such as EPAS1 (Endothelial Per-ARNT-Sim (PAS) domain protein 1), EDNRA (Endothelin-1 receptor Precursor), SIRT1 (Sirtuin type 1), and ryanodine receptor 1 (RYR1) are linked to hypoxia. EPAS1, RYR1, DSG2 (Desomgelin2) are related to cardiomyopathy and PTPRJ (Receptor-type tyrosine-protein phosphatase eta), FUT1 (fucosyltransferase1), HEG1 (heart development protein with EGF-like domains 1), PTPRZ1 (tyrosine phosphatase receptor-type Z polypeptide 1), SIGLEC1 (Sialic acid-binding Ig-like lectin 1 sialoadhesin), NPC1L1 (Niemann-Pick disease-type C1 gene-like 1), and NES (Nestin) are connected with the cardiovascular system, whereas DSG3 is important in maintaining the normal structure and function of hairs (Song et al., 2016). Also, some studies suggested that the DSG3 gene is associated with skin development (Koch et al., 1998; Song et al., 2016).

A high phenotypic variation exists in the Cashmere goat introduced by natural and artificial selection for over 10,000 years, in line with human demand (Naderi et al., 2008; Wang et al., 2016). Owing to evolutionary and artificial selection, the goat has been adapting to different harsh environmental conditions (Song et al., 2016). Multiple genes are described in the significant role of hair and fleece growth in cashmere goats, such as the POU1F1 and PRL genes associated with cashmere yield, diameter, and length, whereas the LHX2 and LIM genes regulate the generation of hair. The FGF5 gene is a key regulator that controls hair length and the FGF9 gene to promote the regeneration of hair follicle after wounding. The WNT2 gene is involved in the initiation of hair follicle (Wang et al., 2016). Furthermore, the DSG1 gene is involved in the communication of hair follicle cell and morphogenesis of hair follicle (Hanakawa et al., 2004; Dong et al., 2013). The DSG4 gene affects wool traits and is responsible for coat color in Cashmere goats (E et al., 2016). Furthermore, the DSG3 gene may have linked with additional fiber growth in the cashmere goats. Other studies have reported that the DSG1, DSG3, and DGS4 genes were expressed in keratinocytes of the basal and immediate suprabasal cell layers (Koch et al., 1998; Bazzi et al., 2006; Dusek et al., 2006). According to Nazari-Ghadikolaei et al. (2018), keratinocytes produce keratin, the main protein for hair, nail, and skin synthesis (Nazari-Ghadikolaei et al., 2018).

The biological function of DSG3 is keratinocyte cell-to-cell adhesion in the basal and suprabasal layers of stratified squamous epithelia. DSG3 has been shown to be the self-antigen for autoantibodies in the pemphigus vulgaris (PV) disease and is known as PV antigen (PVA) (Koch et al., 1997). In past studies of mice and humans, DSG3 was predominantly expressed in the outer layer of skin tissues such as basal and suprabasal cell layers of epidermis (Kljuic et al., 2003), whereas other members of the desmoglein gene family DSG1 are only expressed in suprabasal cell layers and those of DSG2 in the basal cell layer (Koch et al., 1997). DSG3 provide desmogleing compensations with DSG1 in pemphigus foliaceus (PF), which is a potentially fatal autoimmune blistering skin disease in which autoantibodies against DSG3 and DSG1 cause loss of keratinocyte cell adhesion (Payne et al., 2005).

Similarly, desmoglein 3 (DSG3), a SLC24A5 gene, is involved in lighter skin pigmentation in Europeans, showing strong signals of positive selection in the high-altitude populations (Huerta-Sanchez et al., 2013). Furthermore, desmoglein genes (DSG1-4) are expressed in the epidermis and myocardium tissue and their disruption causes some autoimmune diseases that affect the skin, heart, and mucous membranes in humans (Garrod and Chidgey, 2008; Thomason et al., 2010). Desmoglein3 (DSG3) has been identified as one of the autoantigens in an autoimmune blistering skin disease called pemphigus vulgaris (PV) (Stanley and Amagai, 2006; Hartlieb et al., 2013). In this disease, circulating autoantibodies targeting DSG3 induce loss of cell cohesion within the epidermis and mucous membranes. DSG3 was considered to be the new gene marker for terminal differentiation and was found to be a highly upregulated gene during early phase acute hypobaric hypoxia in the rats (Sharma et al., 2015). There were a few studies related to the functional role of DSG3; however, its clear function has not been established yet. It has been suggested that DSG3 may involve more than cell-cell adhesion (Tsang et al., 2010). The expression of the DSG3 gene mostly occurred in the stratified epithelial tissues (Stanley and Amagai, 2006). An epithelial tissue plays an essential role in the diseases of altitude because these tissue cells are the lining of the alveoli of lungs (air sacs) and the blood vessels (endothelium). These tissues have no blood vessels and thus must receive oxygen from the blood vessels in the adjacent connective tissue. For this reason, the tissues are particularly susceptible to damage from low oxygen at high altitude. High-altitude hypoxia impacts structures of vital cells such as sodium and potassium pumps and transcription of the genes slows the sodium and potassium pumps in the alveoli of the lungs. Decreasing the number of these pumps contributes to the accumulation of water in the alveoli and causes high-altitude pulmonary edema (Litch, 2006). Several candidate genes have a key role in controlling the complex cellular processes, energy metabolism, oxygen respiration chain, and reform adenosine triphosphate (ATP) in the cell tissues (Niu et al., 2017). During the evolutionary process, the genes functioning within the oxygen respirations chain have been considered under positive selections. These genes are important for adaptation to high-altitude environments (Xu et al., 2007). Our current study has also shown the importance of Desmosomes genes in high-altitude hypoxia adaptation which has been supported by the SNP data validation analysis.

To our knowledge, there are no research publications related to polymorphism of the DSG3 gene and its role for adaptation to hypoxia at high altitude. In this study, twenty-seven SNPs sites were found from the whole sequence of the DSG3 gene. The allele frequency of variants was observed to be very low in lowland than in highland goat populations. These results have indicated that allele frequencies of variants might be associated with evolution, selection, genetic variation, and ecological conditions (Li et al., 2004a; Wang et al., 2011). Pearson correlation analysis of SNP1, SNP2, SNP3, SNP4, and SNP5 showed a significant relationship between allele frequency and elevation of altitude. These results are in line with the EPAS1 gene in the Cashmere goat and Tibetan dog (Fan et al., 2015; Song et al., 2016).

Similarly, the global F_ST value varied among the Chinese indigenous goat population at all the SNPs sites. The pairwise genetic distances were an ideal index to measure the polymorphism between the groups. Some exons and introns were excluded from further analysis due to a lower pairwise F_ST value. Three non-synonymous and two synonymous mutations have shown large polymorphisms and pairwise genetic distance in the DGS3 gene. Furthermore, this segregated the goat population in distinct branches, as discussed in previously reported studies (Xiang-Long and Valentini, 2004; Di et al., 2011; Ling et al., 2012). Linkage disequilibrium and haplotype analysis showed that the non-synonymous SNPs are highly linked with synonymous SNPs and had higher haplotype frequency variations in the genome region of DSG3 in highland goat populations (Nordborg and Tavaré, 2002; Conrad et al., 2006; Li et al., 2006).

In this study, three non-synonymous mutations are found from the DSG3 gene, out of which one mutation G597S is similar to the mutation G305S at EPAS1 and G14S at HBB (Gou et al., 2014; Fan et al., 2015) in the dog. The G305S is an important mutation with a protein functional change in the PAS domain and has a major effect on the blood flow resistance (Adzhubei et al., 2010; Gou et al., 2014). Moreover, it has been demonstrated as a strongly conserved amino acid site at EPAS1 in the dog (Fan et al., 2015). Therefore, the non-synonymous mutation G597S is likely to affect the function of the DSG3 genes (Figures 4A–C) and suggests a role in high-altitude hypoxic adaptation.

FIGURE 4

FIGURE 4. (A) Protein sequence alignment among different mammalian aligned by cluster omega (Sievers et al., 2011) predicated non-synonymous mutation in Tibetan goat, a conserver site or polymorphism. The dashes indicate deletions. (B) Sanger sequence of blue and black region showed nucleotide position changed in lowland goat, blue indicated nucleotide C, while black indicated nucleotide G. (C) Sanger sequence of red and blue region indicated nucleotide position changed in Tibetan goat, red indicated nucleotide T, while Blue indicated nucleotide C.

Gene mutations and polymorphism are important to understand selection, adaptation, and biological evolutionary processes and to identify the genes related to genetic disease and particular traits (Akey et al., 2002; Nielsen, 2005). The single-nucleotide polymorphism (SNPs) variations in the exonic region of DNA sequences have shown many evolutionary changes. Although the substitution of non-synonymous mutations changed the coding sequence and altered the function of the protein, whereas the synonymous SNP did not change the coding sequence, it may affect the timing of cotranslational protein folding (Sun et al., 2013). In our studies, we have explored five significant SNPs in the most important skin gene DSG3, including three non-synonymous and two synonymous substitutes. They have shown a significant positive linear correlation between mutant allele frequency and elevated altitude in the Chinese indigenous goat population.

Future investigations of the DSG3 gene in goats should focus on the expression level of the gene associated with mutations and altitude. This study has provided the information on the potential function of an allele of DSG3 in high adaptation and strengthens the understanding of adaptation of the Tibetan goat to the extremely harsh environment of high altitude.

Conclusion

Our research has, for the first time, shown the role of the Desmosomes gene for the high-altitude hypoxia adaptation and abundant genetic diversity between goat populations. We found that three non-synonymous candidate mutations (R597E, T595I, and G572S) and two synonymous substitutions significantly segregated these goat populations. Moreover, these results indicated the contribution of DSG3 in high-altitude adaptation and provided new insights for Tibetan cashmere goats.

Author Contributions

CK and LJ wrote the manuscript. XH, QZ, and YP contributed to sample collections and designing the manuscripts. CK and SS has performed the data analysis. YM and LJ conceived the study design. YM, KKM, and AAK interpreted the results and revised the manuscript. All the authors read and approved the final manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Nos. 31472064 and 31601910), the Special Fund for Agro-Scientific Research in the Public Interest (201303059), and the earmarked fund for Modern Agro-industry Technology Research System (CARS-40-01). LJ was supported by the Elite Youth Program in the Chinese Academy of Agricultural Sciences.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00553/full#supplementary-material

Footnotes

1. ^http://sift.bii.a-star.edu.sg/

References

Adzhubei, I. A., Schmidt, S., Peshkin, L., Ramensky, V. E., Gerasimova, A., Bork, P., et al. (2010). A method and server for predicting damaging missense mutations. Nat. Methods 7, 248–249. doi: 10.1038/nmeth0410-248

PubMed Abstract | CrossRef Full Text | Google Scholar

Ai, H., Yang, B., Li, J., Xie, X., Chen, H., and Ren, J. (2014). Population history and genomic signatures for high-altitude adaptation in Tibetan pigs. BMC Genomics 15:834. doi: 10.1186/1471-2164-15-834

PubMed Abstract | CrossRef Full Text | Google Scholar

Akey, J. M., Zhang, G., Zhang, K., Jin, L., and Shriver, M. D. (2002). Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 12, 1805–1814. doi: 10.1101/gr.631202

PubMed Abstract | CrossRef Full Text | Google Scholar

Bandelt, H. J., Forster, P., and Rohl, A. (1999). Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48. doi: 10.1093/oxfordjournals.molbev.a026036

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2005). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi: 10.1093/bioinformatics/bth457

PubMed Abstract | CrossRef Full Text | Google Scholar

Bazzi, H., Getz, A., Mahoney, M. G., Ishida-Yamamoto, A., Langbein, L., Wahl, J. K. III., et al. (2006). Desmoglein 4 is expressed in highly differentiated keratinocytes and trichocytes in human epidermis and hair follicle. Differentiation 74, 129–140. doi: 10.1111/j.1432-0436.2006.00061.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Buroker, N. E., Ning, X. H., Zhou, Z. N., Li, K., Cen, W. J., Wu, X. F., et al. (2012). EPAS1 and EGLN1 associations with high altitude sickness in han and tibetan chinese at the qinghai-tibetan plateau. Blood Cells Mol. Dis. 49, 67–73. doi: 10.1016/j.bcmd.2012.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Q., Qian, X., Lang, Y., Luo, Y., Xu, J., Pan, S., et al. (2013). Genome sequence of ground tit Pseudopodoces humilis and its adaptation to high altitude. Genome Biol. 14:R29. doi: 10.1186/gb-2013-14-3-r29

PubMed Abstract | CrossRef Full Text | Google Scholar

Conrad, D. F., Jakobsson, M., Coop, G., Wen, X., Wall, J. D., Rosenberg, N. A., et al. (2006). A worldwide survey of haplotype variation and linkage disequilibrium in the human genome. Nat. Genet. 38, 1251–1260. doi: 10.1038/ng1911

PubMed Abstract | CrossRef Full Text | Google Scholar

Crawford, D. C., Bhangale, T., Li, N., Hellenthal, G., Rieder, M. J., Nickerson, D. A., et al. (2004). Evidence for substantial fine-scale variation in recombination rates across the human genome. Nat. Genet. 36, 700–706. doi: 10.1038/ng1376

PubMed Abstract | CrossRef Full Text | Google Scholar

Delva, E., Tucker, D. K., and Kowalczyk, A. P. (2009). The desmosome. Cold Spring Harb. Perspect. Biol. 1:a002543. doi: 10.1101/cshperspect.a002543

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, R., Vahidi, S. M., Ma, Y. H., He, X. H., Zhao, Q. J., Han, J. L., et al. (2011). Microsatellite analysis revealed genetic diversity and population structure among chinese cashmere goats. Anim. Genet. 42, 428–431. doi: 10.1111/j.1365-2052.2010.02072.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, Y., Xie, M., Jiang, Y., Xiao, N., Du, X., Zhang, W., et al. (2013). Sequencing and automated whole-genome optical mapping of the genome of a domestic goat (Capra hircus). Nat. Biotechnol. 31, 135–141. doi: 10.1038/nbt.2478

PubMed Abstract | CrossRef Full Text | Google Scholar

Dusek, R. L., Getsios, S., Chen, F., Park, J. K., Amargo, E. V., Cryns, V. L., et al. (2006). The differentiation-dependent desmosomal cadherin desmoglein 1 is a novel caspase-3 target that regulates apoptosis in keratinocytes. J. Biol. Chem. 281, 3614–3624. doi: 10.1074/jbc.M508258200

PubMed Abstract | CrossRef Full Text | Google Scholar

E, G. X., Zhao, Y. J., Ma, Y. H., Cao, G. L., He, J. N., Na, R. S., et al. (2016). Desmoglein 4 diversity and correlation analysis with coat color in goat. Genet. Mol. Res. 15:15017814. doi: 10.4238/gmr.15017814

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, R., Liu, F., Wu, H., Wu, S., Zhu, C., Li, Y., et al. (2015). A positive correlation between elevated altitude and frequency of mutant alleles at the EPAS1 and HBB loci in chinese indigenous dogs. J. Genet. Genomics 42, 173–177. doi: 10.1016/j.jgg.2015.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrod, D., and Chidgey, M. (2008). Desmosome structure, composition and function. Biochim. Biophys. Acta 1778, 572–587. doi: 10.1016/j.bbamem.2007.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, R. L., Cai, Q., Shen, Y. Y., San, A., Ma, L., Zhang, Y., et al. (2013). Draft genome sequence of the Tibetan antelope. Nat. Commun. 4:1858. doi: 10.1038/ncomms2860

PubMed Abstract | CrossRef Full Text | Google Scholar

Gou, X., Wang, Z., Li, N., Qiu, F., Xu, Z., Yan, D., et al. (2014). Whole-genome sequencing of six dog breeds from continuous altitudes reveals adaptation to high-altitude hypoxia. Genome Res. 24, 1308–1315. doi: 10.1101/gr.171876.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanakawa, Y., Li, H., Lin, C., Stanley, J. R., and Cotsarelis, G. (2004). Desmogleins 1 and 3 in the companion layer anchor mouse anagen hair to the follicle. J. Invest. Dermatol. 123, 817–822. doi: 10.1111/j.0022-202X.2004.23479.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartlieb, E., Kempf, B., Partilla, M., Vigh, B., Spindler, V., and Waschke, J. (2013). Desmoglein 2 is less important than desmoglein 3 for keratinocyte cohesion. PLoS One 8:e53739. doi: 10.1371/journal.pone.0053739

PubMed Abstract | CrossRef Full Text | Google Scholar

Huerta-Sanchez, E., Degiorgio, M., Pagani, L., Tarekegn, A., Ekong, R., Antao, T., et al. (2013). Genetic signatures reveal high-altitude adaptation in a set of ethiopian populations. Mol. Biol. Evol. 30, 1877–1888. doi: 10.1093/molbev/mst089

PubMed Abstract | CrossRef Full Text | Google Scholar

Kljuic, A., Bazzi, H., Sundberg, J. P., Martinez-Mir, A., O’Shaughnessy, R., Mahoney, M. G., et al. (2003). Desmoglein 4 in hair follicle differentiation and epidermal adhesion. Cell 113, 249–260. doi: 10.1016/S0092-8674(03)00273-3

CrossRef Full Text | Google Scholar

Koch, P. J., Mahoney, M. G., Cotsarelis, G., Rothenberger, K., Lavker, R. M., and Stanley, J. R. (1998). Desmoglein 3 anchors telogen hair in the follicle. J. Cell Sci. 111( Pt 17), 2529–2537.

PubMed Abstract | Google Scholar

Koch, P. J., Mahoney, M. G., Ishikawa, H., Pulkkinen, L., Uitto, J., Shultz, L., et al. (1997). Targeted disruption of the pemphigus vulgaris antigen (Desmoglein 3) gene in mice causes loss of keratinocyte cell adhesion with a phenotype similar to pemphigus vulgaris. J. Cell Biol. 137:1091.

Google Scholar

Li, C., Wang, Z., Liu, B., Yang, S., Zhu, Z., Fan, B., et al. (2004a). Evaluation of the genetic relationship among ten chinese indigenous pig breeds with twenty-six microsatellite markers. Asian-Australas J Anim Sci 17, 441–444. doi: 10.5713/ajas.2004.441

CrossRef Full Text | Google Scholar

Li, M. H., Li, K., and Zhao, S. H. (2004b). Diversity of chinese indigenous goat breeds: a conservation perspective - a review. Asian-Australas J. Anim. Sci. 17, 726–732. doi: 10.5713/ajas.2004.726

CrossRef Full Text | Google Scholar

Li, M., Tian, S., Jin, L., Zhou, G., Li, Y., Zhang, Y., et al. (2013). Genomic analyses identify distinct patterns of selection in domesticated pigs and tibetan wild boars. Nat. Genet. 45, 1431–1438. doi: 10.1038/ng.2811

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X.-L., Wu, Z.-L., Liu, Z.-Z., Gong, Y.-F., Zhou, R.-Y., and Zheng, G.-R. (2006). SNP identification and analysis in part of intron 2 of goat MSTN gene and variation within and among species. J. Hered. 97, 285–289. doi: 10.1093/jhered/esj026

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Wu, D. D., Boyko, A. R., Wang, G. D., Wu, S. F., Irwin, D. M., et al. (2014). Population variation revealed high-altitude adaptation of tibetan mastiffs. Mol. Biol. Evol. 31, 1200–1205. doi: 10.1093/molbev/msu070

PubMed Abstract | CrossRef Full Text | Google Scholar

Librado, P., and Rozas, J. (2009). DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452. doi: 10.1093/bioinformatics/btp187

PubMed Abstract | CrossRef Full Text | Google Scholar

Ling, Y. H., Zhang, X. D., Yao, N., Ding, J. P., Chen, H. Q., Zhang, Z. J., et al. (2012). Genetic differentiation of chinese indigenous meat goats ascertained using microsatellite information. Asian-Australas J. Anim. Sci. 25, 177–182. doi: 10.5713/ajas.2011.11308

PubMed Abstract | CrossRef Full Text | Google Scholar

Litch, J. A. (2006). Going higher: oxygen, man, and mountains, 5th Edition. Wilderness Environ. Med. 17:e16–e17. doi: 10.1016/S1080-6032(06)70891-2

CrossRef Full Text

Naderi, S., Rezaei, H. R., Pompanon, F., Blum, M. G., Negrini, R., Naghash, H. R., et al. (2008). The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc. Natl. Acad. Sci. U.S.A. 105, 17659–17664. doi: 10.1073/pnas.0804782105

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazari-Ghadikolaei, A., Mehrabani-Yeganeh, H., Miarei-Aashtiani, S. R., Staiger, E. A., Rashidi, A., and Huson, H. J. (2018). Genome-wide association studies identify candidate genes for coat color and mohair traits in the iranian markhoz goat. Front. Genet. 9:105. doi: 10.3389/fgene.2018.00105

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, R. (2005). Molecular signatures of natural selection. Annu. Rev. Genet. 39, 197–218. doi: 10.1146/annurev.genet.39.073003.112420

CrossRef Full Text | Google Scholar

Niu, L., Chen, X., Xiao, P., Zhao, Q., Zhou, J., Hu, J., et al. (2017). Detecting signatures of selection within the tibetan sheep mitochondrial genome. Mitochondrial DNA A DNA Mapp. Seq. Anal. 28, 801–809. doi: 10.1080/24701394.2016.1192614

PubMed Abstract | CrossRef Full Text | Google Scholar

Nordborg, M., and Tavaré, S. (2002). Linkage disequilibrium: what history has to tell us. Trends Genet. 18, 83–90. doi: 10.1016/S0168-9525(02)02557-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Payne, A. S., Ishii, K., Kacir, S., Lin, C., Li, H., Hanakawa, Y., et al. (2005). Genetic and functional characterization of human pemphigus vulgaris monoclonal autoantibodies isolated by phage display. J. Clin. Invest. 115, 888–899. doi: 10.1172/JCI24185

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Q., Zhang, G., Ma, T., Qian, W., Wang, J., Ye, Z., et al. (2012). The yak genome and adaptation to life at high altitude. Nat. Genet. 44, 946–949. doi: 10.1038/ng.2343

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, Y., Zhao, H., Han, N., Zhou, G., Song, G., Gao, B., et al. (2013). Ground tit genome reveals avian adaptation to living at high altitudes in the tibetan plateau. Nat. Commun. 4:2071. doi: 10.1038/ncomms3071

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, P., Bansal, A., and Sharma, P. C. (2015). RNA-seq-based transcriptome profiling reveals differential gene expression in the lungs of sprague-dawley rats during early-phase acute hypobaric hypoxia. Mol. Genet. Genomics 290, 2225–2240. doi: 10.1007/s00438-015-1064-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sievers, F., Wilm, A., Dineen, D., Gibson, T. J., Karplus, K., Li, W., et al. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7:539. doi: 10.1038/msb.2011.75

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, S., Yao, N., Yang, M., Liu, X., Dong, K., Zhao, Q., et al. (2016). Exome sequencing reveals genetic differentiation due to high-altitude adaptation in the tibetan cashmere goat (Capra hircus). BMC Genomics 17:122. doi: 10.1186/s12864-016-2449-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanley, J. R., and Amagai, M. (2006). Pemphigus, bullous impetigo, and the staphylococcal scalded-skin syndrome. N. Engl. J. Med. 355, 1800–1810. doi: 10.1056/NEJMra061111

PubMed Abstract | CrossRef Full Text | Google Scholar

Stephens, M., Smith, N. J., and Donnelly, P. (2001). A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. 68, 978–989. doi: 10.1086/319501

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, J., Zhong, H., Chen, S. Y., Yao, Y. G., and Liu, Y. P. (2013). Association between MT-CO3 haplotypes and high-altitude adaptation in tibetan chicken. Gene 529, 131–137. doi: 10.1016/j.gene.2013.06.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomason, H. A., Scothern, A., McHarg, S., and Garrod, D. R. (2010). Desmosomes: adhesive strength and signalling in health and disease. Biochem. J. 429, 419–433. doi: 10.1042/bj20100567

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, L. G., Yao, T., Mosley-Thompson, E., Davis, M., Henderson, K., and Lin, P.-N. (2000). A high-resolution millennial record of the south asian monsoon from himalayan ice cores. Science 289, 1916–1919. doi: 10.1126/science.289.5486.1916

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsang, S. M., Liu, L., Teh, M. T., Wheeler, A., Grose, R., Hart, I. R., et al. (2010). Desmoglein 3, via an interaction with E-cadherin, is associated with activation of Src. PLoS One 5:e14211. doi: 10.1371/journal.pone.0014211

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G. D., Fan, R. X., Zhai, W., Liu, F., Wang, L., Zhong, L., et al. (2014). Genetic convergence in the adaptation of dogs and humans to the high-altitude environment of the tibetan plateau. Genome Biol. Evol. 6, 2122–2128. doi: 10.1093/gbe/evu162

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Liu, J., Zhou, G., Guo, J., Yan, H., Niu, Y., et al. (2016). Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits. Sci. Rep. 6:38932. doi: 10.1038/srep38932

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Wang, J., Zi, X. D., Huatai, C. R., Ouyang, X., and Liu, L. S. (2011). Genetic diversity of tibetan goats of plateau type using microsatellite markers. Arch. Anim. Breed 54, 188–197. doi: 10.5194/aab-54-188-2011

CrossRef Full Text | Google Scholar

Wei, C., Wang, H., Liu, G., Zhao, F., Kijas, J. W., Ma, Y., et al. (2016). Genome-wide analysis reveals adaptation to high altitudes in tibetan sheep. Sci. Rep. 6:26770. doi: 10.1038/srep26770

PubMed Abstract | CrossRef Full Text | Google Scholar

Weir, B. S., and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. Evolution 38, 1358–1370. doi: 10.2307/2408641

CrossRef Full Text | Google Scholar

Xiang, K., Ouzhuluobu, Peng, Y., Yang, Z., Zhang, X., Cui, C., et al. (2013). Identification of a tibetan-specific mutation in the hypoxic gene EGLN1 and its contribution to high-altitude adaptation. Mol. Biol. Evol. 30, 1889–1898. doi: 10.1093/molbev/mst090

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang-Long, L., and Valentini, A. (2004). Genetic diversity of chinese indigenous goat breeds based on microsatellite markers. J. Anim. Breed. Genet. 121, 350–355. doi: 10.1111/j.1439-0388.2004.00465.x

CrossRef Full Text | Google Scholar

Xu, S., Luosang, J., Hua, S., He, J., Ciren, A., Wang, W., et al. (2007). High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome. J. Genet. Genomics 34, 720–729. doi: 10.1016/s1673-8527(07)60081-2

PubMed Abstract | CrossRef Full Text | Google Scholar