Indian Journal of Animal Research

  • Chief EditorK.M.L. Pathak

  • Print ISSN 0367-6722

  • Online ISSN 0976-0555

  • NAAS Rating 6.50

  • SJR 0.263

  • Impact Factor 0.5 (2023)

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November and December)
Indexing Services :
Science Citation Index Expanded, BIOSIS Preview, ISI Citation Index, Biological Abstracts, Scopus, AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Indian Journal of Animal Research, volume 56 issue 4 (april 2022) : 400-406

​Dominant-genotype Frequency Analysis of Economic Traits Related to SNP Candidate Markers in Three Yak Populations

Q. Zhang1,2, Y.J. Cidan1,2, D.Z. Luosang1,2, Z.D. Pingcuo1,2, Y.L. Dawa1,2, X.Y. Chen1,2, W.D. Basang1,2,*
1State Key Laboratory of Barley and Yak Germplasm Researces and Genetics Improvement, Tibet Academy of Agricultural and Animal Husbandry Science (TAAAS), Lhasa Tibet 850000, China.
2Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agriculture and Animal Husandry Science, Lhasa 850009, China.
Cite article:- Zhang Q., Cidan Y.J., Luosang D.Z., Pingcuo Z.D., Dawa Y.L., Chen X.Y., Basang W.D. (2022). ​Dominant-genotype Frequency Analysis of Economic Traits Related to SNP Candidate Markers in Three Yak Populations . Indian Journal of Animal Research. 56(4): 400-406. doi: 10.18805/IJAR.B-1358.
Background: Yak as a unique domestic animal that has extremely important social value and influence on the local culture and economy in the Qinghai-Tibet Plateau. The current study aimed to evaluate the genotype distribution of a series of economic traits (growth, meat quality and lactation) related to single nucleotide polymorphisms (SNP) in three Tibetan yak populations.

Methods: A total of 238 yaks from three populations [Sibu (SB), Chawula (CWL) and Jiali (JL) yaks] including 34 SB, 104 JL and 100 CWL yak individuals were collected. All samples were genotyped for 12 SNP markers by using SNaPshot technology.

Result: All loci had abundant polymorphisms in the three populations, except for the Hesx1_G618C locus. The dominant growth-related genotype was MyoD1_C1710T (C/T), which had the highest frequency in the three Tibetan yak populations. However, the growth-related dominant genotypes at UCP2_T1499C (C/T) and CYP4A11_G4806A (G/A) loci were rare genotypes in the three Tibetan yak populations. Results of Hardy-Weinberg equilibrium (HWE) analysis showed that all sites did not deviate from the HWE within the population. This finding indicated that these populations belonged to a natural one without having been subjected to artificial selection on economic traits. Overall, this study provided valuable reference for the future molecular breeding of yak based on the genotype distribution of economic-trait candidate markers in three yak populations in Tibet.
(El-Komy et al., 2021) Domestic yak (Bos grunniens) has perfect plateau adaptability and is one of the mammals living in the Qinghai-Tibet Plateau, which has the highest altitude (>3000 m) in the world. This animal provides milk, meat, wool, labor, fuel and other daily necessities for local herdsmen and is thus an important life and economic resource (Nie et al., 2020).
       
The molecular markers (MMs) of economic traits in cattle and buffalo have been largely identified to be related to their growth traits and milk production (Ujan et al., 2011; Lyu et al., 2020; El-Komy et al., 2021). Conversely, only a few candidate genes related to the economic traits of yak have been reported. For example, a single nucleotide polymorphism (SNP; G1069C) from MC4R is significantly (P<0.01) correlated with the body weight of 18-month-old Maiwa yak (Cai et al., 2015). A new SNP (C19913T) in exon 7 of the LPL gene is reportedly responsible for the Phe-to-Ser substitution of amino acids, which have significant effects on the body weight (P<0.01), carcass body weight (P<0.05), eye area and visceral fat weight of yak (Ding et al., 2012).
       
At present, the yak population in the Tibet Plateau of China subsists mostly on extensive grazing and breeding methods. Considering the slow technological development in the region, the production level of yak breeding is low. Therefore, screening the frequency of dominant genotypes of economic-trait-related markers in Tibetan yak populations is needed to assess the breeding potential of these populations for future Tibetan yak breeding. In the present study, we aimed to identify the dominant-genotype distribution of 12 previously known MMs in three Tibetan yak populations and to evaluate the breeding potential for genetic improvement.
Animals and samples
 
A total of 238 yak individuals from three populations were collected, including 34 Sibu (SB), 104 Jiali (JL) and 100 Cawula (CWL) yak individuals (Table 1). Venous blood (5 mL) was collected from each animal and genomic DNA was extracted using a PureLink Genomic DNA Mini Kit (Shanghai, China). Nano Drop 2000 was used to determine DNA concentration and quality.
 

Table 1: Geographic location and sample information of three Tibetan yak populations.


 
Primers and PCR amplification
 
PCR amplification primers were designed using Primer3Plus (http://www.primer3plus.com/) according to the sequence information and previous reports of 12 SNP markers (Table 2).
 

Table 2: PCR amplification and genotyping primers for 12 SNP markers related to the economic traits of yak.


       
The 10 µL PCR amplification system included 5 µL of 2×Taq PCR Mix, 1 µL of Primer Mix, 1 µL of DNA (20 ng/μL) and 3 µL of ddH2O. The PCR amplification-reaction program was as follows: 95°C pre-denaturation for 5 min; 40 cycles each of 94°C denaturation for 30 s, 52°C annealing for 30 s and 72°C extension for 80 s and 72°C extension for 10 min. The PCR product genotyped using SNaPshot technology was referenced as conducted in previous studies (Wang et al., 2021; Cheng et al., 2020).
 
Data analysis
 
A Microsatellite Toolkit (Park, 2012) was used to analyze genotype frequency and genetic-diversity parameters, such as observed heterozygosity (HO), expected heterozygosity (HE) and polymorphism information content (PIC). Arlequin 3.5 software (Excoffier and Lischer, 2010) was used for Hardy-Weinberg equilibrium (HWE) analysis.
Results of genetic-diversity parameter analysis (Table 3) showed that the distribution of HE in the CWL population was 0.000 (Hesx1-G618C and OXGR1-A347G) to 0.4984 (MyoD1-C1710T) and that of HO ranged from 0.0000 (Hesx1-G618C and OXGR1-A347G) to 0.5400 (CAPN_G-1222A). The PIC of the MyoD1_C1710T locus was the highest (0.373). In the JL population, CAPN_G-1222A had the highest HE (0.5023) and HO (0.5098) and CAPN_G-1222A (0.3749) and GHSR_T1387C (0.3731) had the highest PICs. In the SB population, the highest HE was found in GHSR_T1387C (0.5035) and MyoD1_C1710T (0.4965). MyoD1_C1710T and CAPN_G-1222A had the highest HO and five sites, namely, ACSL1_A2079T, GHSR-T1387C, MyoD1_C1710T, CAPN_G-1222A and TMEM-18_C1267T, were medium polymorphic (0.25<PIC<0.5). The HWE results further revealed that no loci considerably deviated from the HWE within the population. This finding implied that the genetic diversity of the three yak populations was well maintained and that the three populations had not undergone artificial selection for economic traits.
 

Table 3: Genetic diversity parameters and Hardy Weinberg equilibrium analysis of 12 SNP markers in three Tibetan yak populations.


       
Table 4 shows the genotype frequencies of 12 SNPs in the three yak populations. The HesxI gene is known to indirectly affect the growth and development of animals (Kaminski et al., 2008). A previous report has found that HesxI_G-618C and HesxI_T226C are significantly correlated (P<0.05) with the weight and height of 12-month-old hornless yak (Zhang et al., 2019). The present study showed that the HesxI_G-618C locus was a homozygous G/G type in the three populations and more than 90% of individuals carried the T/T type of the HesxI_T226C locus. Thus, whether these two markers can be used for the molecular breeding of these three Tibet yaks requires further research.
 

Table 4: Genotypic frequency distribution of 12 SNP markers related to the economic traits of three Tibetan yak populations.


       
Myogenic differentiation 1 (MyoD1) from the muscle regulatory factor (MyoD) family is an important transcription factor that regulates skeletal-muscle production (Talbot et al., 2016). It can activate muscle-gene transcription, transform non-muscle cells into muscle cells, regulate muscle-cell fusion and promote the differentiation of myoblasts into myotubes, subsequently fusing them into muscle fibers and thus regulating embryonic muscle development (Hodge et al., 2019). A large number of genetic variations distributed in the MyoD1 gene region are remarkably related to the growth and meat quality of domestic animals (Liu et al., 2008; Ujan et al., 2011). Specifically, a SNP (MyoD1_C1710T) on exon 3 of the MyoD1 gene has been identified to be significantly correlated with body-size traits (P<0.05) in Shenzha yak and C/T is the dominant genotype (Huang et al., 2019). Herein, we found that the C/T genotype of MyoD1_C1710T had the highest frequency (51.00%-67.64%) within all three yak populations.
       
Cytochrome P450 4A11 (CYP4A11) can be used to diagnose renal cell carcinoma (Kim et al., 2020) and is related to the occurrence of non-alcoholic fatty liver and coronary artery disease (Gao et al., 2020). In particular, a recently discovered SNP marker of the CYP4A11 gene (CYP4A11_G4806A) is significantly correlated with weight and tube circumference in Maiwa yak and G/A genotype is significantly higher in these phenotypes than that of individuals with G/G (P<0.05) (Guan et al., 2019). In the current study, we found that the CYP4A11_G4806A locus had a high frequency of the G/G genotype (73.53%-85.58%) in the three yak populations, the G/A genotype had a certain proportion of distribution, ranging from 14.42% (JL) to 26.47% (SB). Whether the two SNPs of CYP4A11 can be extensively used as MMs of yak growth traits requires further clarification.
       
Transmembrane protein 18 (TMEM18) participates in the migration and regulation of neural stem cells in vivo and in vitro (Jurvansuu et al., 2008) and regulates adipocyte differentiation through the nervous system to affect obesity (Gutierrez-Aguilar et al.,  2011). In previous studies, some mutations from TMEM18 gene have been confirmed to be related to appetite in humans (Larder et al., 2017) and growth traits in cattle (Crispim et al., 2015). Another study has suggested that the TMEM18_C1267T and TMEM18_C44 47T loci are significantly associated with body height and slaughter rate (P<0.05) in Tianzhu yak (Zhang et al., 2017) and that T/T and C/T are the dominant genotype of these loci, respectively. However, the current study found that the two loci from the TMEM-18 gene (C1267T and C4447T) carried the highest frequency of the C/C genotype in the three populations, whereas T/T in TMEM18_C1267T was rare.
       
Calpain 4 (CAPN4) is extensively known to be associated with beef quality traits (Colle et al., 2017). An SNP (CAPN_G-1222A) in the CAPN4 gene’s promoter region is significantly related (P<0.05) to multiple meat quality and growth traits of Gannan yak, including cooking rate and water-loss rate (Niu et al., 2015). This study found that the CAPN_G-1222A locus had abundant polymorphisms in the three populations. Thus, this locus had potential in the molecular breeding of meat-quality traits in the three yak populations.
       
Growth hormone (GH), growth hormone receptor (GHR) and growth hormone secretagogue receptor (GHSR) are important regulatory genes for body growth and development (Wegmann et al., 2017; Wang et al., 2020). GHSR can stimulate the release and secretion of GH and thus participate in the body’s growth and energy metabolism (Zhang et al., 2017; Lv et al., 2018). Particularly, an SNP (T1387C) at the 52  UTR of the GHRS gene is confirmed as significantly related to the body weight in Maiwa yak (P<0.05) (Hai et al., 2017) and the dominant gene has a C/C genotype. In the present study, the frequency of the C/C genotype of GHRS_SNP (T1387C) in the three Tibetan yak populations was 32.69%-45.00%. This finding indicated that the three yak populations carried abundant dominant genotypes in this locus.
       
Uncoupling proteins (UCPs), as an anion carrier member of the inner mitochondrial membrane, maintain energy metabolism and homeostasis in vivo (Krauss et al., 2005). Some studies have shown that UCP2 is associated with body weight in humans and domestic animals (Oliveira et al., 2017; Oguzkan-Balci et al.,  2013). Coincidentally, an SNP (UCP2_T1499C) of UCP2 gene is related to the weight of Maiwa yak and individuals with the C/T genotype have significantly higher weight (P<0.05) than those with the T/T genotype (Hao et al., 2019). In the current work, we found that the dominant genotype C/T was a rare genotype in the three Tibetan yak populations.
       
A series of studies has shown that long-chain acyl-CoA synthetase 1 (ACSL1) is the most important quantitative-trait locus and functional candidate gene that affects fatty acid composition in pork. Moreover, the polymorphism of the ACSL1 gene is significantly related (P<0.05) to the composition and ratio of fatty acid in the skeletal muscle of cattle (Widmann et al., 2011). Some genetic variations in ACSL1 are also substantially related to the milk-production traits of Holstein cows (Liang et al., 2020). In Gannan yaks, the two SNPs (ACSL1_A2079T and ACSL1_G2409A) in the promoter region of the ACSL1 gene are significantly related (P<0.05) to the protein rate and milk-fat rate (Zhao et al., 2019). The present study found that these two loci were rich in polymorphisms in three Tibetan yak populations.
This study analyzed the genotype distribution of 12 candidate SNP genetic markers related to the economic traits of yak in three Tibetan yak populations. Some of them may have value in the breeding selection of Tibet yaks in the future, although extensive and in-depth research is required.
Healthy breeding of Tibetan characteristic livestock (XZ201901NA02), differential expression analysis of several factors related to fat metabolism in Yak (XZNKY-2020-C-007Z03), National Meat cattle and Yak Industrial Technology System (CARS-37), Research on the impact of oregano oil in yak productive performance (XZNKY-2020-C-007Z05), New breed cultivation of Ali yak (XZ202001ZY0021N).

  1. Cai, X., Mipam, T.D., Zhao, F.F., Sun, L. (2015). SNPs detected in the yak MC4R gene and their association with growth traits. Animal. 9(7): 1097-103.

  2. Cheng, S.Z., Guang-Xin, E., Liu, C.L., Basang, W.D., Zhu, Y.B., et al. (2020). SNP of AHSA2 gene in three cattle breeds using snapshot technology. J Genet. 99: 54.

  3. Colle, M.J., Doumit, M.E. (2017). Effect of extended aging on calpain-1 and -2 activity in beef longissimus lumborum and semimembranosus muscles. Meat Science. 131: 142-145. 

  4. Crispim, A.C., Kelly, M.J., Guimarães, S.E., Fonseca e Silva, F., Fortes, M.R., Wenceslau, R.R., Moore, S. (2015). Multi- Trait GWAS and New Candidate Genes Annotation for Growth Curve Parameters in Brahman Cattle. PLoS One. 10(10): e0139906. 

  5. Ding, X.Z., Liang, C.N., Guo, X., Xing, C.F., Bao, P.J., Chu, M., Pei, J., Zhu, X.S., Yan, P. (2012). A novel single nucleotide polymorphism in exon 7 of LPL gene and its association with carcass traits and visceral fat deposition in yak (Bos grunniens) steers. Molecular Biology Reports. 39(1): 669-673. 

  6. El-Komy, S.M., Saleh, A.A., Abd El-Aziz, R.M., El-Magd, M.A. (2021) Association of GH polymorphisms with growth traits in buffaloes. Domestic Animal Endocrinology. 74: 106541. 

  7. Excoffier, L., Lischer, H.E. (2021) Arlequin 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources. 10(3): 564-567. 

  8. Gao, H., Cao, Y., Xia, H., Zhu, X., Jin, Y. (2020). CYP4A11 is involved in the development of nonalcoholic fatty liver disease via ROS-induced lipid peroxidation and inflammation. International Journal of Molecular Medicine. 45(4): 1121-1129.

  9. Guan, C.Y., Wang, J.K., Chai, Z.X., Wang, H., Xin, J.W., Ji, Q.M. and Zhong, J.C. (2019). Correlation between CYP4A11 gene polymorphisms and growth traits in Maiwa Yak. Genomics and Applied Biology. 39(8): 3475-3481.

  10. Gutierrez-Aguilar, R., Kim, D.H., Woods, S.C., Seeley, R.J. (2011). Expression of new loci associated with obesity in diet- induced obese rats: From genetics to physiology. Obesity (Silver Spring). 20: 306-312.

  11. Hai, D., Chai, Z.X., Zhong, J.C. (2017). Identification of SNPs in GH, GHR and GHRS genes and their association with body measurement traits in Maiwa Yak. Acta Veterinaria et Zootechnica Sinica. 48(4): 605-617.

  12. Hao, J., Chai, Z.X., Wang, J.K., Wang, H., Xin, J.W., Ji, Q.M., Zhong, J.C. (2019). Analysis of UCP gene polymorphism and its genetic effects associated with growth traits in Maiwa yak. Heilongjiang Animal Science and veterinary Medicine. 2019(21): 149-153.

  13. Hodge, B.A., Zhang, X., Gutierrez-Monreal, M.A., Cao, Y., Hammers, D.W., Yao, Z., Wolff, C.A., Du, P., Kemler, D., Judge, A.R., Esser, K.A. (2019). MYOD1 functions as a clock amplifier as well as a critical co-factor for downstream circadian gene expression in muscle. Elife. 8: e43017. 

  14. Huang, X., Chai, Z.X., Wang, H., Ji, Q.M., Xin, J.W., Zhong, J.C. (2019). Association analysis between polymorphism of MyoD1 gene and growth traits in Tibetan Yak. Southwest China Journal of Agricultural Sciences. 32(2): 439-445.

  15. Jurvansuu, J., Zhao, Y., Leung, D.S., Boulaire, J., Yu, Y.H., Ahmed, S., Wang, S. (2008). Transmembrane protein 18 enhances the tropism of neural stem cells for glioma cells. Cancer Research. 68(12): 4614-4622.

  16. Kaminski, B.A., Palmert, M.R. (2008). Genetic control of pubertal timing. Current Opinion in Pediatrics. 20(4): 458-464.

  17. Kim, S., Kim, J.M., Lee, H.J., Lim, J.S., Seong, I.O., Kim, K.H. (2020). Alteration of CYP4A11 expression in renal cell carcinoma: Diagnostic and prognostic implications. Journal of Cancer. 11(6): 1478-1485. 

  18. Krauss, S., Zhang, C.Y., Lowell, B.B. (2005). The mitochondrial uncoupling-protein homologues. Nature Reviews Molecular Cell Biology. 6(3): 248-261.

  19. Liang, Y., Gao, Q., Zhang, Q., Arbab, A.A.I., Li, M., Yang, Z., Karrow, N.A., Mao, Y. (2020). Polymorphisms of the ACSL1 gene influence milk production traits and somatic cell score in Chinese Holstein cows. Animals (Basel). 10(12): 2282. 

  20. Liu, M., Peng, J., Xu, D.Q., Zheng, R., Li, F.E., Li, J.L., Zuo, B., Lei, M.G., Xiong, Y.Z., Deng, C.Y., Jiang, S.W. (2008). Association of MYF5 and MYOD1 gene polymorphisms and meat quality traits in Large White x Meishan F2 pig populations. Biochemical Genetics. 46(11-12): 720-732.

  21. Lv, Y., Liang, T., Wang, G., Li, Z. (2018). Ghrelin, a gastrointestinal hormone, regulates energy balance and lipid metabolism. Bioscience Reports. 38(5): BSR20181061. 

  22. Lyu, S., Yang, P., Liu, Y., Song, T., Zhang, Z., Shi, Q., Chen, F., Liu, X., Li, Z., Ru, B., Cai, C., Xie, J., Lei, C., Chen, H., Xu, Z., Wang, E., Huang, Y. (2020). Genetic effects of MOGAT1 gene SNP in growth traits of Chinese cattle. Gene. 769: 145201.

  23. Nie, J., Shao, S., Xia, W., Liu, Z., Yu, C., Li, R., Wang, W., Li, J., Yuan, Y., Rogers, K.M. (2020). Stable isotopes verify geographical origin of yak meat from Qinghai-Tibet plateau. Meat Science.165: 108113. doi: 10.1016/j.meatsci. 2020. 108113.

  24. Niu, X.L., Li, Y.Q., Hu, J., Luo, Y.Z., Guo, S.Z., Yan, W., Yang, S.M. (2015). The mutations in the promoter region of CAPN4 and its effects on carcass and meat quality traits in Yak. Acta Agriculture Boreali-Sinica. 30(2): 28-34.

  25. Oguzkan-Balci, S., Col-Araz, N., Nacak, M., Araz, M., Sabanci, H., Balat, A., Pehlivan, S. (2013). Mitochondrial uncoupling protein 2 (ucp2) gene polymorphisms are associated with childhood obesity and related metabolic disorders. Journal of Pediatric Endocrinology and Metabolism Jpem. 26: 277-283.

  26. Oliveira, B.A.P.D., Pinhel, M.A.D.S., Nicoletti, C.F., Oliveira, C.C.D., Nonino, C.B. (2017). Ucp2 and plin1 expression affects the resting metabolic rate and weight loss on obese patients. Obesity Surgery. 27(2): 343-348.

  27. Park, S. (2012). Excel Microsatellite Toolkit, version 3.1.1. Animal Genomics Lab website. University College, Dublin. 

  28. Talbot, J., Maves, L. (2016). Skeletal muscle fiber type: Using insights from muscle developmental biology to dissect targets for susceptibility and resistance to muscle disease. Wiley Interdisciplinary Reviews-Developmental Biology. 5(4): 518-34.

  29. Ujan, J.A., Zan, L.S., Wang, H.B., Ujan, S.A., Adoligbe, C., Wang, H.C., Biao, S.F. (2011). Lack of an association between a single nucleotide polymorphism in the bovine myogenic determination 1 (MyoD1) gene and meat quality traits in indigenous Chinese cattle breeds. Genetics and Molecular Research. 10(3): 2213-22. 

  30. Wang, M., Qian, L., Li, J., Ming, H., Fang, L., Li, Y., Zhang, M., Xu, Y., Ban, Y., Zhang, W., Zhang, Y., Liu, Y., Wang, N. (2020). GHSR deficiency exacerbates cardiac fibrosis: Role in macrophage inflammasome activation and myofibroblast differentiation. Cardiovascular Research. 116(13): 2091- 2102.

  31. Wang, X., E, G.X., Na, R.S., Liu, C.L., Huang, Y.F. (2021). Single nucleotide polymorphism of cdc37, ahsa1 and stip1 gene in three cattle breeds using snapshot technology. Indian Journal of Animal Research. DOI: 10.18805/IJAR.B-1266.

  32. Wegmann, M.G., Thankamony, A., Roche, E., Hoey, H., Kirk, J., Shaikh, G., Ivarsson, S.A., Söder, O., Dunger, D.B., Juul, A., Jensen, R.B. (2017). The exon3-deleted growth hormone receptor gene polymorphism (d3-GHR) is associated with insulin and spontaneous growth in short SGA children (NESGAS). Growth Hormone and Igf Research. 35: 45-51. 

  33. Widmann, P., Nuernberg, K., Kuehn, C., Weikard, R. (2011). Association of an ACSL1 gene variant with polyunsaturated fatty acids in bovine skeletal muscle. BMC Genetics. 12: 96.

  34. Zhang, B., Shang, P., Tao, Z., Qiangba, Y., Wang, Z., Zhang, H. (2017). Effect of a single nucleotide polymorphism in the growth hormone secretagogue receptor (GHSR) gene on growth rate in pigs. Gene. 634: 68-73. 

  35. Zhang, H., Wang, H.B., Ma, W., Jia, C.J., Ma, X.M., Wu, X.Y., Chu, M., Yan, P., Cheng, S.R., Liang, C.N. (2019). Association analysis of HesxI gene polymorphism and growth traits in polled yak. Journal of Henan Agricultural Science. 48(8): 134-139.

  36. Zhang, H., Zhang, Q.W., Wang, Q., Ma, Y.J., Zhang, Y., Zhao, X.X. (2017). Polymorphism of TMEM-18 gene and its correlation with production traits in Yak. Biotechnology Bulletin. 33(2): 89-96.

  37. Zhao, Z.D., Tian, H.S., Jiang, Y.Y., Shi, B.G., Liu, X., Li, X.P., Wang, D.Z., Chen, J.L., Hu, J. (2019). Polymorphisms of ACSL1 gene promoter and their association analysis with milk quality traits in Yak (Bos grunniens). Journal of Agricultural Biotechnology. 27(9): 1596-1630.

Editorial Board

View all (0)