- 1Department of Medical College, Shihezi University, Shihezi, China
- 2Department of Respiratory and Critical Care Medicine, First People’s Hospital of Kashi, Kashi, China
Genetic factors are important factors in chronic obstructive pulmonary disease (COPD) onset. Plenty of risk and new causative genes for COPD have been identified in patients of the Chinese Han population. In contrast, we know considerably little concerning the genetics in the Kashi COPD population (Uyghur). This study aims at clarifying the genetic maps regarding COPD susceptibility in Kashi (China). Whole-exome sequencing (WES) was used to analyze three Uyghur families with COPD in Kashi (eight patients and one healthy control). Sanger sequencing was also used to verify the WES results in 541 unrelated Uyghur COPD patients and 534 Uyghur healthy controls. WES showed 72 single nucleotide variants (SNVs), two deletions, and small insertions (InDels), 26 copy number variants (CNVs), and 34 structural variants (SVs), including g.71230620T > A (rs12449210T > A, NC_000,016.10) in the HYDIN axonemal central pair apparatus protein (HYDIN) gene and g.61190482A > G (rs777591A > G, NC_000002.12) in the ubiquitin-specific protease 34 (USP34) gene. After Sanger sequencing, we found that rs777591“AA” under different genetic models except for the dominant model (adjusted OR = 0.8559, 95%CI 0.6568–1.115, p > .05), could significantly reduce COPD risk, but rs12449210T > A was not related to COPD. In stratified analysis of smoking status, rs777591“AA” reduced COPD risk significantly among the nonsmoker group. Protein and mRNA expression of USP34 in cigarette smoke extract-treated BEAS-2b cells increased significantly compared with those in the control group. Our findings associate the USP34 rs777591“AA” genotype as a protector factor in COPD.
Introduction
Chronic obstructive pulmonary disease (COPD) is a common and complicated disease of the lungs. It is characterized by continuous and irreversible airflow destruction due to chronic inflammation (Lozano et al., 2012). As of 2020, COPD is expected to become the fifth most prevalent burden of disease and the third leading cause of death worldwide (Vestbo et al., 2013).
Our epidemiological investigation shows the COPD prevalence in Kashi (Xinjiang, China) to be 17.01% (Li et al., 2021), which is higher than that in other parts of China. COPD development involves environmental factors (e.g., smoking, air pollution), genetic susceptibility, and infection (Zhang et al., 2014). A trend of familial aggregation in COPD patients is also documented (Wang et al., 2020). Scholars have found that some genes might be related to COPD genetic susceptibility of the Chinese Han population (Wang C. et al., 2018; Zhang et al., 2020). In contrast, the study for COPD is still lacking in Kashi.
Kashi city is located in the northwestern China, and more than 90% of the total local population in Kashi are Uyghurs (Zhang et al., 2018). The Uyghurs demonstrate a range of mixed Asian and European anthropological features (Xu and Jin 2008). At the same time, they usually do not marry other ethnic groups. Hence, their genomes are significantly different compared with the Han populations. Furthermore, the history, geographic location, and local customs of Kashi are quite different from other regions of China (Abuzhalihan et al., 2019; Abudureheman et al., 2021). These features make the Kashi Uyghurs a resourceful population for describing the ethnicity-specific variants associated with COPD.
The HYDIN axonemal central pair apparatus protein (HYDIN) gene is located on human chromosome 16 with a length of 423 kb (Laske et al., 2013). Population studies show that HYDIN mutations can cause primary ciliary dyskinesia (PCD) (Shapiro et al., 2018). PCD is a genetic disease with abnormal ciliary motility and is characterized by chronic respiratory infections (Olbrich et al., 2012; Kurkowiak et al., 2015). Respiratory infection is known to play an important role in the pathogenesis and progression of COPD (Leung et al., 2017).
The ubiquitin-specific protease 34 (USP34) gene is a member of the ubiquitin-specific protease family. USPs can regulate cell growth (Liu et al., 2019) and inhibit apoptosis (de Castro et al., 2019). Moreover, USP34 regulates the Wnt pathway (Lui et al., 2011) and plays an important role in DNA damage (Sy et al., 2013). Studies confirm that USP34 also has an influence in the NF-κB pathway (Truong et al., 2021). DNA damage, excessive apoptosis, and the NF-κB pathway are known to be associated with the development of COPD (Neofytou et al., 2012; Sauler et al., 2018; Canadas et al., 2020; Wang et al., 2021).
In the present study, eight people with COPD and one healthy person from three Uyghur families with COPD in Kashi were subjected to whole-exome sequencing (WES) to screen for the susceptibility genes and polymorphisms related to COPD. The two single nucleotide variants (SNVs) (rs12449210T > A in HYDIN and rs777591A > G in USP34) were determined to be studied because they were found for the first time in the Kashi COPD population. Furthermore, the association between the two SNVs and COPD risk have not been previously described. Based on the WES data sets, we continued to evaluate their relationship with clinical characteristics in a case-control study of 1075 people. Then, we hypothesized that the presence of SNVs combined with environmental factors (e.g., smoking) might regulate gene expression, thereby increasing the risk of COPD. Also, protein and mRNA expression of them were determined in cigarette smoke extract (CSE)-treated bronchial epithelial beas-2b (BEAS-2b) cells.
Materials and Methods
Collection of Information of COPD Families for WES
Eight COPD patients and one healthy control (HC) were selected for WES. All nine people were permanent residents of Kashi and aged >40 years. These Uyghur families were derived from a previous epidemiological investigation of COPD in Kashi. Families had ≥2 generations with COPD in three generations. The inclusion and exclusion criteria for the subjects were described by Gong et al. (2022).
Peripheral blood samples (5 ml) were obtained from each participant and transferred to EDTA-K2 Vacutainer tubes for DNA extraction. Apart from blood samples, the basic information of patients and HC were also collected: age, sex, body mass index (BMI), pulmonary function, smoking history, and other data. All pulmonary function tests were undertaken using a spirometer (Cosmed, Rome, Italy).
Collection of Samples for a Case-Control Study
Recruited individuals were aged >40 years and from Kashi during 2018–2019. The study cohort was 541 Uyghur unrelated COPD patients and 534 Uyghur HCs from the First People’s Hospital of Kashi. For the COPD group, the inclusion criteria were people whose forced expiratory volume in the first second (FEV1)/forced vital capacity (FVC) < 0.70 after bronchodilator inhalation denotes airflow limitation. The exclusion criteria are identical to COPD family patients. For the HC group, the exclusion criteria were the same as HC of COPD family.
After providing written informed consent, all individuals were required to provide the same basic information as the family subjects. Similarly, 5 ml peripheral blood samples were obtained from each of the 1075 subjects for DNA extraction.
Genetic Analyses
Genetic analyses, including WES, data analyses, variants selection, filtering strategy, Sanger sequencing, enrichment analyses, chemicals and reagents, Western blotting, and real-time RT-qPCR, had been described in the supplemental methods.
Statistical Analyses
SPSS 18 (IBM, Armonk, NY, United States) and PLINK v1.07 (pngu.mgh.harvard.edu/∼purcell/plink/index.shtml) were employed for statistical analyses. Quantitative data are the mean ± SD or median (interquartile range). The independent t-test was used to compare the difference between age and BMI. We used the chi-square test to ascertain the influence of sex, smoking status, coal consumption, wood consumption, cigarettes per day, and cumulative quantity of active smoking. The difference in quitting smoking years was assessed by Fisher’s exact test. For parameters that did not have a normal distribution (e.g., annual household income, FEV1%, and FEV1/FVC), we used the Mann–Whitney U-test to evaluate differences. We tested for the Hardy–Weinberg equilibrium on each SNV among samples in the case-control study. Akaike’s information criteria (AIC) were used to calculate the genetic model for rs12449210 T > A and rs777591A > G (including genotype, dominant, recessive, allele, and additive). We used odds ratios (ORs) and corresponding 95% confidence intervals (CIs) by logistic regression analysis after correcting for sex, age, BMI, and smoking status to assess the relationship between SNVs and case–control groups. R function genpwr in package devtools was used to perform post hoc power estimation. The homogeneity test of stratified analyses was evaluated by the chi-square test. In the enrichment analyses, verification criteria of multiple SNVs were evaluated using the false discovery rate (FDR) with multiple corrections. The experimental data of Western blotting and real-time RT-qPCR were the mean ± SEM, and differences between groups were evaluated by one-way ANOVA. The replicates for the experiments were at least three times. p < .05 was considered significant.
Results
Family Members With COPD
WES
WES was undertaken for Uyghur family members with COPD to identify common and significant variants. The genealogical trees of the three Uyghur families are described in Figure 1A. Eight patients and one healthy subject were included in Supplementary.Table S1. The family members with COPD were aged 48–90 (mean age = 63.63) and had always lived in Kashi.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792027/fcell-09-792027-HTML/image_m/fcell-09-792027-g001.jpg)
FIGURE 1. Family information and variants. (A): Genealogical trees of three Uyghur families with COPD. Square = male, and circle = female. A square with a slash represents a dead male, and a circle with a slash denotes a dead female. A solid square represents a male with COPD, and the solid circle denotes a female with COPD. A solid square with “R” represents a COPD male who refused to join our study, and the solid circle with “R” denotes a COPD female who refused to join our study. (B): Summary of all variations among exomes of family members with COPD. There are nine circles of information represented from the outside to the inside of this. Panel 1. Chromosome. 2. Map of SNVs. The density calculation is based on the number of SNVs in each window using log10. The color changes from red, yellow, and blue. 3. Density graph of inserts. The density calculation is based on the number of inserts in each window using log10. 4. Density of deletions. The density calculation is based on the number of deletions in each window using log10. 5. Density map of the variation sites in coding regions, including SNVs and InDels. The density calculation is based on the various sites in each window using log10. 6. Density map of the mutation sites in the noncoding area, including SNVs and InDels. The density calculation is based on the number of mutation sites in each window using log10. 7. Map of CNV locations. The area size represents the size of CNV. Red = gain, blue = loss (when the number of samples >1, there is no circle). 8. Map of SV position. The area size represents the size of SVs. Orange = DEL, green = INS. 9. Association of SV type diagram. Blue = INV, red = CTX, green = ITX.
After WES, there were ∼0.26 Tb of sequencing data. The minimum value of sequencing depth was 30×, the average sequencing depth was 116.29 ×. Q20 of all samples was >95%, Q30 was >90%, and the average coverage rate was >99.5%. We detected 100,193 SNVs, 11,009 deletions, and small insertions (InDels), 2915 copy number variants (CNVs), and 29,016 structural variants (SVs). To facilitate observation of all variants in family members with COPD, we used Circos to display all variants (Figure 1B).
Analysis for SNVs
After WES, we obtained 100,193 SNVs. We undertook statistical analyses of all SNVs (Supplementary Figure S1). We found that the top three regions with the largest number of SNVs were the intron (43.47%), exon (42.52%), and noncoding RNA (5.76%) (Supplementary Figures S1A–C).
Then, 3096 SNVs were obtained after calculating the corresponding scores, SIFT and PolyPhen-2 (Supplementary Figure S1G). Finally, we obtained 72 SNVs after the filtering strategy (Supplementary Table S2).
Analysis for Small InDels, CNVs, and SVs
The top three regions with the largest number of InDels were the intron (71.90%), exon (8.28%), and noncoding RNA (8.28%) (Supplementary Figures S1D,E). Finally, we obtained 26 InDels (Supplementary Table S3).
After drawing a heat map for all 2915 CNVs, we detected the gain and loss of each chromosome (Supplementary Figure S2). Finally, we obtained two CNVs (Supplementary Table S4).
After analyzing the distribution of SVs in different regions of the genome, we found that the top three regions with the largest number of SVs were the intron (55.99%), exon (18.65%), and 3′ UTR (8.24%) (Supplementary Figure S1F). Finally, we obtained 34 SVs (Supplementary Table S5).
In conclusion, we obtained 72 SNVs, two CNVs, 26 InDels, and 34 SVs by the aforementioned filtering strategy.
Enrichment Analyses of Common Mutant Genes in Family Members with COPD
We selected the top 30 q-value GO terms, pathways, and diseases to draw a bubble chart (Figures 2A–C). We also selected the top 30 q-value GO terms among “cell component,” “molecular function,” and “biological process” to draw bubble charts (Figures 2D–F). The enrichment analyses indicate common mutant genes were related to cilium- or flagellum-dependent cell motility, ciliary cytoplasm, axoneme, calcium-ion binding, olfactory transduction, ABC transporters, obesity-related traits, and congenital disorders of metabolism.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792027/fcell-09-792027-HTML/image_m/fcell-09-792027-g002.jpg)
FIGURE 2. Bubble charts. (A–C): The top 30 q-values for functional (GO database), pathway (KEGG database), and disease enrichment. (D–F): The top 30 q-value GO terms among the categories of “cell component,” “molecular function,” and “biological process.” The abscissa “rich factor” represents the ratio of input or background frequency in the enrichment-analysis result. The bubble size represents the number of genes annotated to this functional entry for the mutant gene, and the color corresponds to the q-value in the enrichment-analysis result.
Case-Control Study
Clinical Characteristics of the COPD and Non-COPD Groups
A total of 1075 individuals were invited to participate (Table 1). The study cohort comprised 541 Uyghur COPD patients and 534 Uyghur HCs. Significant differences between COPD cases and HCs were observed with regard to age, sex, BMI, smoking status, FEV1%, and FEV1/FVC (p < .05 for all) but not for annual household income, coal consumption, wood consumption, cigarettes per day, cumulative quantity of active smoking, and quitting smoking years (p > .05 for all) (Supplementary Table S6).
Hardy–Weinberg Equilibrium of rs12449210T > A and rs777591A > G
Rs12449210T > A and rs777591A > G both met the criteria for the Hardy–Weinberg equilibrium (p > .05) (Table 2). Therefore, rs12449210T > A and rs777591A > G could be analyzed further.
Analysis of Genotypes of rs12449210T > A and rs777591A > G in the Case-Control Study
We undertook genetic model analysis (genotype, dominant, recessive, allele, and additive) on the SNVs of rs12449210T > A and rs777591A > G (Table 2). Because a few subjects might have de novo mutations (Francioli et al., 2015; Besenbacher et al., 2016; Garimella et al., 2020) in the primer binding region, the call rate of rs12449210 of HYDIN and rs777591 of USP34 were 98.88% (1,063/1,075) and 99.44% (1,069/1,075).
For rs777591A > G located on USP34, after correcting for sex, age, smoking status, and BMI, this site significantly decreased the risk of COPD with the “AA” genotype based on the genotype model (OR: 0.5767, 95%CI: 0.3716–0.8948, p < .05) and recessive model (OR: 0.5942, 95%CI: 0.3927–0.8993, p < 0.05). In the allele model, the “A” genotype was related to a reduced risk of COPD (OR: 0.8188, 95%CI: 0.6746–0.9938, p < .05). For the additive model, a reduced tendency of COPD risk was also present (OR: 0.8146, 95%CI: 0.6693–0.9915, p < .05).
For rs12449210T > A located on HYDIN before correction “TA” or “AA + TA” increased the COPD risk under the genotype and dominant models (genotype model: OR: 1.345, 95%CI: 1.018–1.779, p < .05; dominant model: OR: 1.304, 95%CI: 1.004–1.695, p < .05). After correction, a significant effect was not present in the five models (p > .05).
Stratified Analysis of rs12449210T > A and rs777591A > G in the Case–Control Study
We undertook stratified analysis for the relationship between rs12449210T > A and rs777591A > G and COPD risk based on smoking status and FEV1%.
The beneficial effects of “AA” in rs777591A > G was more evident in nonsmokers in the genotype model [adjusted odds ratio (aOR): 0.511, 95%CI: 0.3102–0.8417], recessive model (aOR: 0.514, 95%CI: 0.3208–0.8237), allele model (aOR: 0.7978, 95%CI: 0.6415–0.9922), and additive model (aOR: 0.8053, 95%CI: 0.6507–0.9964). However, in the smoker group, the protective effects were not significant (Table 3).
For rs12449210T > A (Supplementary Table S7), no differences were observed in smokers and nonsmokers before or after correction of potentially influencing factors.
Rs12449210T > A and rs777591A > G did not show an obvious difference in FEV1% (Supplementary Table S8).
Expression of USP34 and Iκbα in CSE-Treated BEAS-2b Cells
We wished to determine the impact of CSE on expression of USP34 and Iκbα. We measured expression in BEAS-2b cells treated with 0%–3% CSE for 24 h. Compared with the control group, 3% CSE obviously increased expression of USP34 and Iκbα protein (Figures 3E,F) and mRNA (Figures 3A,B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792027/fcell-09-792027-HTML/image_m/fcell-09-792027-g003.jpg)
FIGURE 3. Expression of USP34 and Iκbα in CSE-treated BEAS-2b cells. mRNA expression assessed by RT-qPCR for USP34 (A,C) or Iκbα (B,D) after treatment with (A,B) 0%–3% CSE for 24 h (n = 3) or (C,D) 3% CSE for different times (n = 3). Protein expression assessed by Western blotting for USP34 (E,G) or Iκbα (F,H) after treatment with (E,F) 0%–3% CSE (n = 3) or (G,H) 3% CSE for different times (n = 3). Anti-β-actin antibody was used as a control in Western blotting. *p < .05, **p < .01, and ***p < .001, versus control or 24-h group.
To determine the impact of different treatment times on expression of USP34 and Iκbα, we measured expression in BEAS-2b cells treated with 3% CSE for 24, 48, and 72 h. Compared with the 24 h group, treatment for 48 and 72 h significantly increased expression of USP34 and Iκbα protein (Figures 3G,H) and mRNA (Figures 3C,D).
Discussion
A total of 72 SNVs, two CNVs, 26 InDels, and 34 SVs were obtained by WES from three Uyghur families with COPD. Two SNVs, rs12449210T > A of HYDIN and rs777591A > G of USP34, were chosen to verify in a large Uyghur population. We found that “AA” in rs777591A > G was a protective factor against COPD, whereas rs12449210T > A was not related to COPD susceptibility in Kashi COPD population. mRNA and protein expression of USP34 and Iκbα were obviously increased in CSE-treated BEAS-2b cells in vitro study.
WES has been used widely in the discovery of pathogenic genes (Adalsteinsson et al., 2017; Epi25 Collaborative, 2019; Petrovski et al., 2019). WES could be used to detect high-frequency, low-frequency, and even rare variants that have important roles in the occurrence of complex diseases and could be employed to discover new variants (Meienberg et al., 2015).
An approach, used in Mendelian diseases of high penetrance, is to use a direct filtering strategy in a few families to identify rare or novel variants (Qiao et al., 2016). Many studies identify and validate the candidate variants that were discovered using this approach. They performed WES in a small sample size of family patients to identify potentially related causative genes. The candidate variants were then assessed in a large sample size of case-control population by Sanger sequencing (Mescheriakova et al., 2019; Cetani et al., 2020; Engelbrecht et al., 2020; Froukh et al., 2020).
We undertook WES on eight Uyghur patients from three families with COPD. We showed a DNA-variation map of COPD patients (SNVs, InDels, CNVs, and SVs) at the genetic level, and found 72 SNVs of 55 genes that might be associated with COPD. Scholars have found four genes that might be related to COPD using WES in 49 families with COPD: DNAH8, ALCAM, RARS, and GBF1. The main populations they studied were Caucasian and African Americans (Regan et al., 2010; Qiao et al., 2016). Also, the patients they studied were aged <53 years with FEV1% ≤40% before bronchodilator inhalation. In our study, we studied a Uyghur population aged >40 years with FEV1/FVC <0.70 after bronchodilator inhalation. Those were the differences between ours and previous studies. Furthermore, the final 72 SNVs of 55 genes for COPD identified in our study were different from those determined in previous studies. On the one hand, previous studies selected rare variants, according to the common disease–rare mutation hypothesis (Khera et al., 2018) and most candidate gene methods (Gibson 2012). We did not select rare variants as the research targets. On the other hand, the genetic backgrounds of different ethnicities are different, and even though the same filtering strategy was used, the results might be different (An et al., 2016). Consequently, rs12449210T > A of HYDIN and rs777591A > G of USP34 in our study were found for the first time in a Uyghur COPD population in Kashi.
We evaluated only one HC as a reference; the selected 72 SNVs might contain the SNVs of healthy people. Hence, it was necessary to verify the relationship between these 72 SNVs and COPD in a case-control study.
In the case-control study of two SNVs (rs12449210T > A and rs777591A > G), they were analyzed through genetic models, tobacco smoking, and pulmonary function. Because the PLINK software defaulted a higher mutation frequency to a wild type, in Asia, the “G” mutant frequency of rs777591A > G was higher than “A” (via NCBI database: https://www.ncbi.nlm.nih.gov/). Therefore, our data determined “AA” in rs777591A > G (USP34) to be a protective factor.
Smoking is the main risk factor for COPD (Vogelmeier et al., 2017); rs777591“AA” showed weak protection. However, in the stratified analysis of smokers and nonsmokers, the protective effect of rs777591“AA” was meaningful only in the nonsmoker group. Based on this analysis, we think rs777591“AA” of USP34 had only a weak protective effect against COPD in the nonsmoker group.
In the present study for the COPD population in Kashi, the allele “A” frequency of rs777591A > G is 35.13%, and the COPD prevalence is 17.01% (Li et al., 2021). In East Asia, the frequency of allele “A” is 32%, but the prevalence of COPD in China (East Asia), Korea (East Asia), and Japan (East Asia) is 13.7% (Wang Y. et al., 2018), 7.82% (Kwon and Kim 2016), and 10.9% (Tan and Ng 2008), respectively. Although South Asia has a higher allele “A” frequency (45%), the prevalence of COPD in countries such as India and Nepal is 4.1% and 18% (Tan and Ng 2008), respectively. Thus, despite the similar frequency of allele “A”, the incidence of COPD varied considerably. This might be because the pathogenesis of COPD was not caused by genes alone, but by the interaction of genetic factors, environmental factors (e.g., tobacco smoking, air pollution, wood and coal consumption) and economic status (Shetty et al., 2021). The protective effect of “A” in rs777591A > G might only be demonstrated in the Kashi COPD population at present, and further validation in different ethnic groups or different regions is required.
HYDIN mainly encodes the C2b protein on the central microtubules of motile cilia to adjust the amplitude of the cilia swing and the synergy between cilia (Olbrich et al., 2012). The influence of HYDIN on the cilia ultrastructure could lead to abnormal cilia function and, ultimately, reduce the ability to remove foreign bodies (Lechtreck et al., 2008). Some studies show that, in COPD patients, the cilia on the mucosal surface of the respiratory tract are shortened, and cilia motility is reduced, resulting in significant impairment of foreign-body excretion (Hedstrom et al., 2021). Therefore, we speculated that HYDIN might have a key role in the occurrence and development of COPD. In this study, the rs12449210T > A of HYDIN located in 5′ UTR. HYDIN might be related to the abnormal function of respiratory-tract mucosal cilia in COPD, but we showed in a case-control study that HYDIN had little effect on the occurrence and development of COPD. This did not mean that HYDIN was not associated with COPD in other ethnicities.
Studies show that USP34 inhibits nuclear factor-kappa B (NF-κB) signal transduction by deubiquitinating and stabilizing the NF-κB inhibitor Iκbα (Li et al., 2020). Studies on the inflammatory cell surface and bronchial epithelial biopsies in COPD cases show that the NF-κB pathway is highly activated (Di Stefano et al., 2002). NF-κB expression has also been found to be increased significantly in the lung tissues of animals with COPD (Yang et al., 2009). Therefore, USP34 might have a key role in the occurrence and development of COPD. rs777591A > G is located in the intron region at HaploReg v4.1 (ANNOVAR annotated it in 3′UTR of USP34). Although rs777591A > G is located in the noncoding regions, genome-wide association studies (GWASs) indicate that only 7% single nucleotide polymorphisms (SNPs) related to complex diseases were located in the coding regions, whereas 93% were in the noncoding regions (Yang et al., 2019). We showed that rs777591“AA” could reduce COPD, so it is a protective factor against COPD. In our in vitro study, mRNA and protein expression of USP34 and Iκbα in CSE-treated BEAS-2b cells were obviously higher than those in controls.
We note important limitations of our experimental setting. First, only one non-COPD sample of a COPD family was used to be a control in WES. This may result in containing family-specific variants that are not associated with COPD. It would increase costs and reduce experimental efficiency for case-control validation. Second, the recruited participants were all Uyghurs, and the results might not be representative for ethnic populations. Also, the sample size of the smoking population in the case-control group was small. In addition, although USP34 mRNA and protein expression levels were increased in CSE-treated BEAS-2b cells, the specific role of USP34 in COPD pathogenesis is unclear.
Due to the above limitations, we will conduct multicenter studies and expand the sample size to further validate our data. Further experiments (including cellular experiments in genetic models, knockout or overexpressed animal experiments) are also required.
Conclusion
In the present study, WES revealed 72 SNVs, two CNVs, 26 InDels, and 34 SVs in three families with COPD. rs777591 “AA” of USP34 was a protective factor against COPD, especially in the nonsmoking population. USP34 mRNA and protein expression levels were increased in CSE-treated BEAS-2b cells. Therefore, our data provide new clues for the relationship between USP34 polymorphisms and COPD susceptibility in a Chinese Uyghur population.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: PRJNA785331
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethics Committee of the First People’s Hospital of Kashi (2019-95). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
XZ and LL conceived and designed the study. JX and JR carried out the experiments and wrote the first version of the manuscript. XZ, CX, and AZ conducted the variants analyses. AA, MT, and SZ recruited individuals and undertook clinical work. LT and DH carried out statistical analyses. All authors read and approved this version of the manuscript.
Funding
This study was supported by the Xinjiang Uygur Autonomous Natural Science Foundation of China (Grant numbers: 2019D01C006 and 2017D01C009).
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all the participants in this study and Hui Gong for statistical analyses as well as general assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.792027/full#supplementary-material
References
Abudureheman, Z., Li, L., Zhong, X., Xu, J., Gong, H., Yilamujiang, S., et al. (2021). The rs74794265 SNP of the SREK1 Gene is Associated with COPD in Kashi, China. Copd 16, 2631–2636. doi:10.2147/COPD.S321150
Abuzhalihan, J., Wang, Y.-T., Adi, D., Ma, Y.-T., Fu, Z.-Y., Yang, Y.-N., et al. (2019). Prevalence of Dyslipidemia in Students from Han, Uygur, and Kazakh Ethnic Groups in a Medical University in Xinjiang, China. Sci. Rep. 9, 19475. doi:10.1038/s41598-019-55480-5
Adalsteinsson, V. A., Ha, G., Freeman, S. S., Choudhury, A. D., Stover, D. G., Parsons, H. A., et al. (2017). Scalable Whole-Exome Sequencing of Cell-free DNA Reveals High Concordance with Metastatic Tumors. Nat. Commun. 8, 1324. doi:10.1038/s41467-017-00965-y
An, L., Lin, Y., Yang, T., and Hua, L. (2016). Exploring the Interaction Among EPHX1, GSTP1, SERPINE2, and TGFB1 Contributing to the Quantitative Traits of Chronic Obstructive Pulmonary Disease in Chinese Han Population. Hum. Genomics. 10, 13. doi:10.1186/s40246-016-0076-0
Besenbacher, S., Sulem, P., Helgason, A., Helgason, H., Kristjansson, H., Jonasdottir, A., et al. (2016). Multi-nucleotide de novo Mutations in Humans. Plos. Genet. 12, e1006315. doi:10.1371/journal.pgen.1006315
Cañadas, O., Olmeda, B., Alonso, A., and Pérez-Gil, J. (2020). Lipid-Protein and Protein-Protein Interactions in the Pulmonary Surfactant System and Their Role in Lung Homeostasis. Int. J. Mol. Sci. 21, 3708. doi:10.3390/ijms21103708
Cetani, F., Pardi, E., Aretini, P., Saponaro, F., Borsari, S., Mazoni, L., et al. (2020). Whole Exome Sequencing in Familial Isolated Primary Hyperparathyroidism. J. Endocrinol. Invest. 43, 231–245. doi:10.1007/s40618-019-01107-5
de Castro, G. S., Simoes, E., Lima, J. D. C. C., Ortiz-Silva, M., Festuccia, W. T., Tokeshi, F., et al. (2019). Human Cachexia Induces Changes in Mitochondria, Autophagy and Apoptosis in the Skeletal Muscle. Cancers 11, 1264. doi:10.3390/cancers11091264
Di Stefano, A., Caramori, G., Oates, T., Capelli, A., Lusuardi, M., Gnemmi, I., et al. (2002). Increased Expression of Nuclear Factor- B in Bronchial Biopsies from Smokers and Patients with COPD. Eur. Respir. J. 20, 556–563. doi:10.1183/09031936.02.00272002
Engelbrecht, H.-R., Dalvie, S., Agenbag, G., Stein, D. J., and Ramesar, R. S. (2020). Whole-exome Sequencing in an Afrikaner Family with Bipolar Disorder. J. Affective Disord. 276, 69–75. doi:10.1016/j.jad.2020.06.045
Epi25 Collaborative (2019). Ultra-Rare Genetic Variation in the Epilepsies: A Whole-Exome Sequencing Study of 17,606 Individuals. Am. J. Hum. Genet. 105, 267–282. doi:10.1016/j.ajhg.2019.05.020
Francioli, L. C., Polak, P. P., Polak, P. P., Koren, A., Menelaou, A., Chun, S., et al. (2015). Genome-wide Patterns and Properties of de novo Mutations in Humans. Nat. Genet. 47, 822–826. doi:10.1038/ng.3292
Froukh, T., Hawwari, A., and Al Zubi, K. (2020). Whole Exome Sequencing Highlights Variants in Association with Keratoconus in Jordanian Families. BMC. Med. Genet. 21, 177. doi:10.1186/s12881-020-01112-z
Garimella, K. V., Iqbal, Z., Krause, M. A., Campino, S., Kekre, M., Drury, E., et al. (2020). Detection of Simple and Complex de novo Mutations with Multiple Reference Sequences. Genome Res. 30, 1154–1169. doi:10.1101/gr.255505.119
Gibson, G. (2012). Rare and Common Variants: Twenty Arguments. Nat. Rev. Genet. 13, 135–145. doi:10.1038/nrg3118
Gong, H., Ren, J., Xu, J., Zhong, X., Abudureheman, Z., Yilamujiang, S., et al. (2022). SMAD3 rs36221701 T>C Polymorphism Impacts COPD Susceptibility in the Kashi Population. Gene 808, 145970. doi:10.1016/j.gene.2021.145970
Hedström, U., Öberg, L., Vaarala, O., Dellgren, G., Silverborn, M., Bjermer, L., et al. (2021). Impaired Differentiation of COPD Bronchial Epithelial Cells Grown on Bronchial Scaffolds. Am. J. Respir. Cell. Mol. Biol. 65, 201–213. doi:10.1165/rcmb.2019-0395OC
Khera, A. V., Chaffin, M., Aragam, K. G., Haas, M. E., Roselli, C., Choi, S. H., et al. (2018). Genome-wide Polygenic Scores for Common Diseases Identify Individuals with Risk Equivalent to Monogenic Mutations. Nat. Genet. 50, 1219–1224. doi:10.1038/s41588-018-0183-z
Kurkowiak, M., Ziętkiewicz, E., and Witt, M. (2015). Recent Advances in Primary Ciliary Dyskinesia Genetics. J. Med. Genet. 52, 1–9. doi:10.1136/jmedgenet-2014-102755
Kwon, H., and Kim, E. (2016). Factors Contributing to Quality of Life in COPD Patients in South Korea. Copd 11, 103–109. doi:10.2147/COPD.S90566
Laske, K., Shebzukhov, Y. V., Grosse-Hovest, L., Kuprash, D. V., Khlgatian, S. V., Koroleva, E. P., et al. (2013). Alternative Variants of Human HYDIN Are Novel Cancer-Associated Antigens Recognized by Adaptive Immunity. Cancer Immunol. Res. 1, 190–200. doi:10.1158/2326-6066.CIR-13-0079
Lechtreck, K.-F., Delmotte, P., Robinson, M. L., Sanderson, M. J., and Witman, G. B. (2008). Mutations in Hydin Impair Ciliary Motility in Mice. J. Cell. Biol. 180, 633–643. doi:10.1083/jcb.200710162
Leung, J. M., Tiew, P. Y., Mac Aogáin, M., Budden, K. F., Yong, V. F. L., Thomas, S. S., et al. (2017). The Role of Acute and Chronic Respiratory Colonization and Infections in the Pathogenesis of COPD. Respirology 22, 634–650. doi:10.1111/resp.13032
Li, L., Zhong, X., Zheng, A., JianKun, C., Budukadeer, A. A., Aini, P., et al. (2021). Prevalence and Risk Factors of Chronic Obstructive Pulmonary Disease in Kashi Region, Northwestern China. Copd 16, 655–663. doi:10.2147/COPD.S289620
Li, Q., Wang, M., Xue, H., Liu, W., Guo, Y., Xu, R., et al. (2020). Ubiquitin‐Specific Protease 34 Inhibits Osteoclast Differentiation by RegulatingNF‐κBSignaling. J. Bone. Miner. Res. 35, 1597–1608. doi:10.1002/jbmr.4015
Liu, C., Yao, X., Li, M., Xi, Y., and Zhao, L. (2019). USP39 Regulates the Cell Cycle, Survival, and Growth of Human Leukemia Cells. Biosci. Rep. 39. doi:10.1042/BSR20190040
Lozano, R., Naghavi, M., Foreman, K., Lim, S., Shibuya, K., Aboyans, V., et al. (2012). Global and Regional Mortality from 235 Causes of Death for 20 Age Groups in 1990 and 2010: a Systematic Analysis for the Global Burden of Disease Study 2010. Lancet 380, 2095–2128. doi:10.1016/S0140-6736(12)61728-0
Lui, T. T. H., Lacroix, C., Ahmed, S. M., Goldenberg, S. J., Leach, C. A., Daulat, A. M., et al. (2011). The Ubiquitin-specific Protease USP34 Regulates Axin Stability and Wnt/-Catenin Signaling. Mol. Cell Biol. 31, 2053–2065. doi:10.1128/MCB.01094-10
Meienberg, J., Zerjavic, K., Keller, I., Okoniewski, M., Patrignani, A., Ludin, K., et al. (2015). New Insights into the Performance of Human Whole-Exome Capture Platforms. Nucleic Acids Res. 43, e76. doi:10.1093/nar/gkv216
Mescheriakova, J. Y., Verkerk, A. J., Amin, N., Uitterlinden, A. G., van Duijn, C. M., and Hintzen, R. Q. (2019). Linkage Analysis and Whole Exome Sequencing Identify a Novel Candidate Gene in a Dutch Multiple Sclerosis Family. Mult. Scler. 25, 909–917. doi:10.1177/1352458518777202
Neofytou, E., Tzortzaki, E., Chatziantoniou, A., and Siafakas, N. (2012). DNA Damage Due to Oxidative Stress in Chronic Obstructive Pulmonary Disease (COPD). Int. J. Mol. Sci. 13, 16853–16864. doi:10.3390/ijms131216853
Olbrich, H., Schmidts, M., Werner, C., Onoufriadis, A., Loges, N. T., Raidt, J., et al. (2012). Recessive HYDIN Mutations Cause Primary Ciliary Dyskinesia without Randomization of Left-Right Body Asymmetry. Am. J. Hum. Genet. 91, 672–684. doi:10.1016/j.ajhg.2012.08.016
Petrovski, S., Aggarwal, V., Giordano, J. L., Stosic, M., Wou, K., Bier, L., et al. (2019). Whole-Exome Sequencing in the Evaluation of Fetal Structural Anomalies: a Prospective Cohort Study. The Lancet 393, 758–767. doi:10.1016/S0140-6736(18)32042-7
Qiao, D., Lange, C., Beaty, T. H., Crapo, J. D., Barnes, K. C., Bamshad, M., et al. (2016). Exome Sequencing Analysis in Severe, Early-Onset Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 193, 1353–1363. doi:10.1164/rccm.201506-1223OC
Regan, E. A., Hokanson, J. E., Murphy, J. R., Make, B., Lynch, D. A., Beaty, T. H., et al. (2010). Genetic Epidemiology of COPD (COPDGene) Study Design. COPD: J. Chronic Obstructive Pulm. Dis. 7, 32–43. doi:10.3109/15412550903499522
Sauler, M., Lamontagne, M., Finnemore, E., Herazo-Maya, J. D., Tedrow, J., Zhang, X., et al. (2018). The DNA Repair Transcriptome in Severe COPD. Eur. Respir. J. 52, 1701994. doi:10.1183/13993003.01994-2017
Shapiro, A., Davis, S., Manion, M., and Briones, K. (2018). Primary Ciliary Dyskinesia (PCD). Am. J. Respir. Crit. Care Med. 198, P3–P4. doi:10.1164/rccm.1982P3
Shetty, B. S. P., D'Souza, G., and andPadukudru Anand, M. (2021). Effect of Indoor Air Pollution on Chronic Obstructive Pulmonary Disease (COPD) Deaths in Southern Asia-A Systematic Review and Meta-Analysis. Toxics 9. doi:10.3390/toxics90400810.3390/toxics9040085
Sy, S. M. H., Jiang, J., O, W. S., Deng, Y., and Huen, M. S. Y. (2013). The Ubiquitin Specific Protease USP34 Promotes Ubiquitin Signaling at DNA Double-Strand Breaks. Nucleic Acids Res. 41, 8572–8580. doi:10.1093/nar/gkt622
Tan, W. C., and Ng, T. P. (2008). COPD in Asia. Chest 133, 517–527. doi:10.1378/chest.07-11310.1378/chest.07-1131
Truong, L., Zheng, Y.-M., Kandhi, S., and Wang, Y.-X. (2021). Overview on Interactive Role of Inflammation, Reactive Oxygen Species, and Calcium Signaling in Asthma, COPD, and Pulmonary Hypertension. Adv. Exp. Med. Biol. 1304, 147–164. doi:10.1007/978-3-030-68748-9_9
Vestbo, J., Hurd, S. S., Agustí, A. G., Jones, P. W., Vogelmeier, C., Anzueto, A., et al. (2013). Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 187, 347–365. doi:10.1164/rccm.201204-0596PP
Vogelmeier, C. F., Criner, G. J., Martinez, F. J., Anzueto, A., Barnes, P. J., Bourbeau, J., et al. (2017). Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Lung Disease 2017 Report. GOLD Executive Summary. Am. J. Respir. Crit. Care Med. 195, 557–582. doi:10.1164/rccm.201701-0218PP
Wang C., C., Xu, J., Yang, L., Xu, Y., Zhang, X., Bai, C., et al. (2018). Prevalence and Risk Factors of Chronic Obstructive Pulmonary Disease in China (The China Pulmonary Health [CPH] Study): A National Cross-Sectional Study. The Lancet 391, 1706–1717. doi:10.1016/S0140-6736(18)30841-9
Wang, J., Sun, D., Lu, W., Zhang, Z., Zhang, C., and Hu, K. (2020). BRIP1 rs10744996C>A Variant Increases the Risk of Chronic Obstructive Pulmonary Disease in the Mongolian Population of Northern China. Exp. Physiol. 105, 1950–1959. doi:10.1113/EP088210
Wang, X., Huang, L., Jiang, S., Cheng, K., Wang, D., Luo, Q., et al. (2021). Testosterone Attenuates Pulmonary Epithelial Inflammation in Male Rats of COPD Model through Preventing NRF1-Derived NF-κB Signaling. J. Mol. Cell. Biol. 13, 128–140. doi:10.1093/jmcb/mjaa079
Wang, Y., Zhou, Q., Dong, L., Xiong, M., Jiang, H., Guo, M., et al. (2018). The Effects of CXCL10 Polymorphisms on COPD Susceptibility. Mol. Genet. Genomics. 293, 649–655. doi:10.1007/s00438-017-1408-z
Xu, S., and Jin, L. (2008). A Genome-wide Analysis of Admixture in Uyghurs and a High-Density Admixture Map for Disease-Gene Discovery. Am. J. Hum. Genet. 83, 322–336. doi:10.1016/j.ajhg.2008.08.001
Yang, S.-R., Yao, H., Rajendrasozhan, S., Chung, S., Edirisinghe, I., Valvo, S., et al. (2009). RelB Is Differentially Regulated by IκB Kinase-α in B Cells and Mouse Lung by Cigarette Smoke. Am. J. Respir. Cell. Mol. Biol. 40, 147–158. doi:10.1165/rcmb.2008-0207OC
Yang, T., Li, J., Wen, Y., Tan, T., Yang, J., Pan, J., et al. (2019). LINC00673 rs11655237 C>T Polymorphism Impacts Hepatoblastoma Susceptibility in Chinese Children. Front. Genet. 10, 506. doi:10.3389/fgene.2019.00506
Zhang, W., Li, Y., Lin, J., Abduryim, A., and Zhao, J. (2018). Cariogenicity of Candida Albicans of Distinct Genotypes Among 3-5-Year-Old Uygur Children in Kashgar, China- a Case-Control Study. BMC. Oral Health 18, 203. doi:10.1186/s12903-018-0658-4
Zhang, Z., Wang, J., Chen, L., Chen, X., Sun, G., Zhong, N., et al. (2014). Impact of Haze and Air Pollution-Related Hazards on Hospital Admissions in Guangzhou, China. Environ. Sci. Pollut. Res. 21, 4236–4244. doi:10.1007/s11356-013-2374-6
Keywords: whole-exome sequencing, HYDIN, USP34, rs12449210, rs777591, COPD
Citation: Xu J, Li L, Ren J, Zhong X, Xie C, Zheng A, Abudukadier A, Tuerxun M, Zhang S, Tang L, Hairoula D and Zou X (2022) Whole-Exome Sequencing Implicates the USP34 rs777591A > G Intron Variant in Chronic Obstructive Pulmonary Disease in a Kashi Cohort. Front. Cell Dev. Biol. 9:792027. doi: 10.3389/fcell.2021.792027
Received: 09 October 2021; Accepted: 08 December 2021;
Published: 07 February 2022.
Edited by:
Enrique Medina-Acosta, State University of the North Fluminense Darcy Ribeiro, BrazilReviewed by:
Cíntia Barros Santos-Rebouças, Rio de Janeiro State University, BrazilGlauber Dias, State University of the North Fluminense Darcy Ribeiro, Brazil
Copyright © 2022 Xu, Li, Ren, Zhong, Xie, Zheng, Abudukadier, Tuerxun, Zhang, Tang, Hairoula and Zou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaoguang Zou, 3521600081@qq.com
†These authors have contributed equally to this work