1
|
Liu X, Chen W, Huang B, Wang X, Peng Y, Zhang X, Chai W, Khan MZ, Wang C. Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits. Front Vet Sci 2024; 10:1334434. [PMID: 38274664 PMCID: PMC10808162 DOI: 10.3389/fvets.2023.1334434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 12/27/2023] [Indexed: 01/27/2024] Open
Abstract
Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Muhammad Zahoor Khan
- Liaocheng Research Institute of Donkey High-Efficiency Breeding, Liaocheng University, Liaocheng, China
| | - Changfa Wang
- Liaocheng Research Institute of Donkey High-Efficiency Breeding, Liaocheng University, Liaocheng, China
| |
Collapse
|
2
|
Kong L, Yue Y, Li J, Yang B, Chen B, Liu J, Lu Z. Transcriptomics and metabolomics reveal improved performance of Hu sheep on hybridization with Southdown sheep. Food Res Int 2023; 173:113240. [PMID: 37803553 DOI: 10.1016/j.foodres.2023.113240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 07/04/2023] [Accepted: 07/05/2023] [Indexed: 10/08/2023]
Abstract
Consumers are increasingly demanding high-quality mutton. Cross breeding can improve meat quality and is widely used in sheep breeding. However, little is known about the molecular mechanism of cross breeding sheep meat quality. In this study, male Southdown and female Hu sheep were hybridized. The slaughter performance and longissimus dorsi quality of the 6-month-old hybrid offspring were measured, and the longissimus dorsi of the hybrid offspring was analyzed by transcriptomics and metabolomics to explore the effect of cross breeding on meat quality. The results showed that the production performance of Southdown × Hu F1 sheep was significantly improved, the carcass fat content was significantly decreased, and the eating quality of Southdown × Hu F1 sheep were better. Compared with the HS group (Hu × Hu), the NH group (Southdown × Hu) had 538 differentially expressed genes and 166 differentially expressed metabolites (P < 0.05), which were significantly enriched in amino acid metabolism and other related pathways. Up-regulated genes METTL21C, PPARGC1A and down-regulated gene WFIKKN2 are related to muscle growth and development. Among them, the METTL21C gene, which is related to muscle development, was highly correlated with carnosine, a metabolite related to meat quality (correlation > 0.6 and P < 0.05). Our results provide further understanding of the molecular mechanism of cross breeding for sheep muscle growth and meat quality optimization.
Collapse
Affiliation(s)
- Lingying Kong
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Yaojing Yue
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Jianye Li
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Bohui Yang
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Bowen Chen
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
| | - Jianbin Liu
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.
| | - Zengkui Lu
- Key Laboratory of Animal Genetics and Breeding on the Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; Sheep Breeding Engineering Technology Research Center of Chinese Academy of Agricultural Sciences, Lanzhou 730050, China.
| |
Collapse
|
3
|
Wang Q, Wei Z, Zhu H, Pan C, Akhatayeva Z, Song X, Lan X. Goat Pleomorphic Adenoma Gene 1 ( PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits. Animals (Basel) 2023; 13:2023. [PMID: 37370533 DOI: 10.3390/ani13122023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Revised: 05/19/2023] [Accepted: 06/12/2023] [Indexed: 06/29/2023] Open
Abstract
The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.
Collapse
Affiliation(s)
- Qian Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Zhenyu Wei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Haijing Zhu
- Shaanxi Provincial Engineering and Technology Research Center of Cashmere Goats, Yulin University, Yulin 719000, China
- Life Science Research Center, Yulin University, Yulin 719000, China
| | - Chuanying Pan
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Zhanerke Akhatayeva
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Xiaoyue Song
- Shaanxi Provincial Engineering and Technology Research Center of Cashmere Goats, Yulin University, Yulin 719000, China
- Life Science Research Center, Yulin University, Yulin 719000, China
| | - Xianyong Lan
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| |
Collapse
|
4
|
Wu H, Zhou W, Liu H, Cui X, Ma W, Wu H, Li G, Wang L, Zhang J, Zhang X, Ji P, Lian Z, Liu G. Whole-genome methylation analysis reveals epigenetic variation between wild-type and nontransgenic cloned, ASMT transgenic cloned dairy goats generated by the somatic cell nuclear transfer. J Anim Sci Biotechnol 2022; 13:145. [PMID: 36434676 PMCID: PMC9701027 DOI: 10.1186/s40104-022-00764-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Accepted: 08/03/2022] [Indexed: 11/26/2022] Open
Abstract
BACKGROUND SCNT (somatic cell nuclear transfer) is of great significance to biological research and also to the livestock breeding. However, the survival rate of the SCNT cloned animals is relatively low compared to other transgenic methods. This indicates the potential epigenetic variations between them. DNA methylation is a key marker of mammalian epigenetics and its alterations will lead to phenotypic differences. In this study, ASMT (acetylserotonin-O-methyltransferase) ovarian overexpression transgenic goat was produced by using SCNT. To investigate whether there are epigenetic differences between cloned and WT (wild type) goats, WGBS (whole-genome bisulfite sequencing) was used to measure the whole-genome methylation of these animals. RESULTS It is observed that the different mCpG sites are mainly present in the intergenic and intronic regions between cloned and WT animals, and their CG-type methylation sites are strongly correlated. DMR (differentially methylated region) lengths are located around 1000 bp, mainly distributed in the exonic, intergenic and intronic functional domains. A total of 56 and 36 DMGs (differentially methylated genes) were identified by GO and KEGG databases, respectively. Functional annotation showed that DMGs were enriched in biological-process, cellular-component, molecular-function and other signaling pathways. A total of 10 identical genes related to growth and development were identified in GO and KEGG databases. CONCLUSION The differences in methylation genes among the tested animals have been identified. A total of 10 DMGs associated with growth and development were identified between cloned and WT animals. The results indicate that the differential patterns of DNA methylation between the cloned and WT goats are probably caused by the SCNT. These novel observations will help us to further identify the unveiled mechanisms of somatic cell cloning technology, particularly in goats.
Collapse
Affiliation(s)
- Hao Wu
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China ,Sany Institute of China Agricultural University, Sanya, 572025 China
| | - Wendi Zhou
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Haijun Liu
- Institute of Animal Husbandry and Veterinary, Academy of Agricultural Sciences of Tianjin, Tianjin, 300192 China
| | - Xudai Cui
- Qingdao Senmiao Industrial Co., Ltd., Qingdao, 266101 China
| | - Wenkui Ma
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Haixin Wu
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Guangdong Li
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Likai Wang
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Jinlong Zhang
- Institute of Animal Husbandry and Veterinary, Academy of Agricultural Sciences of Tianjin, Tianjin, 300192 China
| | - Xiaosheng Zhang
- Institute of Animal Husbandry and Veterinary, Academy of Agricultural Sciences of Tianjin, Tianjin, 300192 China
| | - Pengyun Ji
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Zhengxing Lian
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Guoshi Liu
- grid.22935.3f0000 0004 0530 8290National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China ,Sany Institute of China Agricultural University, Sanya, 572025 China
| |
Collapse
|
5
|
Hu L, Zhang L, Li Q, Liu H, Xu T, Zhao N, Han X, Xu S, Zhao X, Zhang C. Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing. Front Genet 2022; 13:971464. [PMID: 36160022 PMCID: PMC9490000 DOI: 10.3389/fgene.2022.971464] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Accepted: 08/23/2022] [Indexed: 01/29/2023] Open
Abstract
Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.
Collapse
Affiliation(s)
- Linyong Hu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Liangzhi Zhang
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Qi Li
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Hongjin Liu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Tianwei Xu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Na Zhao
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Xueping Han
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
- Technology Extension Service of Animal Husbandry of Qinghai, Xining, China
| | - Shixiao Xu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Xinquan Zhao
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Cunfang Zhang
- State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining, China
- *Correspondence: Cunfang Zhang,
| |
Collapse
|
6
|
Wang X, Wang Y, Cao X, Huang Y, Li P, Lan X, Buren C, Hu L, Chen H. Copy number variations of the KAT6A gene are associated with body measurements of Chinese sheep breeds. Anim Biotechnol 2021:1-8. [PMID: 34842492 DOI: 10.1080/10495398.2021.2005616] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Copy number variation (CNV) is one kind of genomic structure variations and presents as gains and losses of genomic fragments. More recently, we have made an atlas of CNV maps for livestock. In the future, it is a primary focus to determine the phenotypic effects of candidate CNVs. Lysine Acetyltransferase 6 A (KAT6A) is a protein coding gene and plays a critical role in many cellular processes. However, the effects of KAT6A CNVs on sheep body measurements remains unknown. In this study, we performed quantitative polymerase chain reaction (qPCR) to detect the presences and distributions of three CNV regions within KAT6A gene in 672 sheep from four Chinese breeds. Association analysis indicated that the three CNVs of KAT6A gene were significantly associated with body measurement(s) in Small-tailed Han sheep (STH) and Hu sheep (HU) (p < 0.05), while no effects on Large-tailed Han sheep (LTH) were observed (p > 0.05) were observed. Additionally, only one CNV was significantly associated with body measurement (body length) in Chaka sheep (CK) (p < 0.05). Our study provided evidence that the CNV(s) of KAT6A gene could be used as candidate marker(s) for molecular breedings of STH, HU, and CK breeds.
Collapse
Affiliation(s)
- Xiaogang Wang
- Key laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Yiru Wang
- Key laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Xiukai Cao
- Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou, China
| | - Yongzhen Huang
- Key laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Pi Li
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, Qinghai, China
| | - Xianyong Lan
- Key laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Chaogetu Buren
- Animal Disease Control Center of Haixi Mongolian and Tibetan Autonomous Prefecture, Delingha, Qinghai, China
| | - Linyong Hu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, Qinghai, China
| | - Hong Chen
- Key laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China.,College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| |
Collapse
|