Association between MPHOSPH6 gene polymorphisms and IgA nephropathy risk in a Chinese Han population

Multiple genetic and environmental factors together contribute to the risk of IgA nephropathy (IgAN). MPHOSPH6 play an important role in the recruitment of the exosome to the pre-rRNA. However, to date, little information is found about the association between MPHOSPH6 polymorphisms and the IgAN risk. In this case-control study, we genotyped five single nucleotide polymorphisms (SNPs) in MPHOSPH6 gene in 416 IgAN cases and 495 controls using Sequenom Mass-ARRAY technology and evaluated their association with IgAN using the χ2 and genetic model analysis. In the allelic model analysis, we determined rs1056654 was associated with a 0.774-fold decrease in the risk of IgAN (95%CI= 0.630-0.952; p = 0.015). In the genetic model analysis, we found that the “C/C” genotype of rs1056675 was associated with an increased risk of IgAN based on the codominant model (OR =1.48; 95% CI=1.03-2.13; p=0.033) and recessive model (OR =1.52; 95% CI=1.11-2.09; p=0.0095). The “G/A-A/A” genotype of rs1056654 was associated with a decreased risk of IgAN based on the dominant model (OR =0.75; 95% CI=0.58-0.98; p=0.032) and log-additvie model (OR =0.78; 95% CI=0.64-0.96; p=0.0188). Our data suggested that gene polymorphisms in the MPHOSPH6 may exert influences IgAN susceptibility in a Chinese Han population.


INTRODUCTION
IgA nephropathy (IgAN) is the most common type of primary glomerulonephritis [1], and more than 40% of IgAN patients will eventually worsen to end stage renal disease [2]. However, the prevalence of IgAN is significantly different among ethnicities. It is most common in individuals of Asian ancestry, with a very high prevalence of 3.7% [3], but is rare among Africans [4] and is of intermediate prevalence among Europeans (up to 1.3%) [5]. It has been widely accepted that multiple genetic and environmental factors together contribute to the risk of IgAN.
In the past decades, several genome-wide association studies (GWAS) have linked a plenty variety of genetic abnormality in the onset and development of IgAN. Among them, complement factor H related 1 and 3 (CFHR1, CFHR3) [6], tumor necrosis factor (TNFSF13) [7] and a-defensin (DEFA), major histocompatibility complex (MHC) [8] genes polymorphisms were the most common susceptibility genes with the risk of IgAN.
Recently, microRNAs such as miR-148b [9], miR-n-9 [10] were associated with the development of IgAN. These results remind us that study the polymorphisms of genes involving in formation of RNA may provide us a better understanding of IgAN.

Research Paper
It has been demonstrated that MPHOSPH6 play an important role in the recruitment of the exosome to the pre-rRNA [11]. MPHOSPH6 also might be involved in regulating shrimp cell cycle and ovary development [12]. However, to date, little information is found about the association between MPHOSPH6 polymorphisms and the risk of IgAN. In the present study, we selected five single nucleotide polymorphisms (SNPs) of MPHOSPH6 to investigate the association between MPHOSPH6 polymorphisms and the risk of IgAN in a Chinese Han population.

RESULTS
A total of 416 IgAN cases (271 men and 145 women; mean age, 33.35 ±12.13 years) and 495 controls (180 men and 315 women; mean age, 54.48 ±9.44 years) were included in the study. The basic characteristics of the cases and controls are shown in Table 1.
The MAFs of the analyzed SNPs in the case and control groups are shown in Table 2. All SNPs were in Hardy-Weinberg equilibrium (HWE) in the controls (p > 0.05) with the exception of rs1056629, which was excluded from subsequent analyses. The MAFs of the SNPs in the control group were similar to those reported for the HapMap Asian population. Using chi-square tests, we determined that rs1056654 was associated with a 0.774-fold decrease in the risk of IgAN (95% confidence interval [CI] = 0.630-0.952; p = 0.015).
The genotype frequencies of the MPHOSPH6 polymorphisms are shown in Table 3. Compared with the TT genotype, the CC frequency of rs1056675 polymorphism among cases was significantly different from controls (CC versus TT: OR =1.48, 95% CI=1.03-2.13, p=0.035), which suggested that the rs1056675 polymorphism had an increased effect on the risk of IgAN. Additionally, we observed that rs1056654 polymorphism have a protective role of IgAN risk (AA versus GG: OR =0.61, 95% CI=0.38-0.98, p=0.045).
Furthermore, we assumed that the minor allele of each SNP as a risk factor compared with the wild-type allele. Four models (codominant, dominant, recessive, and log-additive) were applied to analyze the associations between the SNPs and risk of IgAN using an unconditional

DISCUSSION
In this study, we investigated the associations between five selected MPHOSPH6 SNPs and risk of IgAN in the Chinese Han population of Shaanxi Province. We found that rs1056675 is associated with an increased risk of IgAN, while rs1056654 has a protective role for IgAN. Our results suggest that the polymorphisms of MPHOSPH6 may play an important role in the risk of IgAN in the Chinese Han population.
MPHOSPH6 is located on chromosome 16q23.3 and encodes the M-phase phosphoprotein 6 (MPP6) that is important for the maturation of 5.8S rRNA. MPP6 is also a RNA-binding protein, which preferentially binds to pyrimidine homopolymers [13]. In fact, little study is found about the association between MPHOSPH6 and disease. Li et al has attempted to investigate the association between MPHOSPH6 polymorphisms and risk of colorectal cancer, but they have no significant result [14]. In our study, we detected five SNPs of MPHOSPH6, and found that genetic polymorphisms of MPHOSPH6 are associated with IgAN, which may shed a new light on the in-depth study for this gene.
Previous association studies have found many SNPs associated with IgAN risk in some populations; however, little of them were successfully replicated in other populations. Earlier GWAS studies revealed two common variants rs3115573 and rs3130315 in MHC that influence IgAN risk in a British population [15]. Chinese researchers also attempted to investigate SNPs associated with IgAN in the Chinese populations, they identified novel SNPs including rs3803800, rs2738048, rs660895, rs1794275 and rs2523946 reached a genome-wide significance [8]. It has been revealed that four SNPs including rs9275224, rs2856717, rs9275424 and rs9275596 were significantly associated with IgAN risk in a GWAS replication study in a combined Chinese, European and African-American population and followed by a meta-analysis in 85 world populations [16]. In the present study, we found that rs1056675 and rs1056654 were significantly associated with IgAN risk in a Chinese Han population. As far as we know, we are the first to report the association between MPHOSPH6 polymorphisms and IgAN risk, so the results identified here should be confirmed in further studies.
Our study had several intrinsic limitations. Firstly, the subjects in this study were all Han Chinese who lived in Shaanxi Province. The results identified here   need to be confirmed in other ethnicities. Secondly, our sample size was relatively small, so we didn't do further subgroup analyses based on age or gender. Thirdly, IgAN is a heterogeneous disease with many other risk factors. We could not explore the interactions between genetic polymorphisms and environmental factors in IgAN patients due to the limited data. Therefore, the relationship between MPHOSPH6 polymorphisms and IgAN risk must be evaluated in future studies with bigger sample size and different populations.

Subjects
We consecutively recruited 416 IgAN cases and 495 healthy controls from the First Affiliated Hospital of Xi'an Jiaotong University from March 2014 to December 2016. All subjects were of northern Han Chinese ancestry and were recruited among local inhabitants of Shanxi province from northern China. All the patients were diagnosed and histologically confirmed to suffer from IgAN according to the renal biopsy, and they had not received any systemic treatment before the time of examination. Besides, we excluded patients with any type of cancer, infection, secondary IgAN (Secondary IgAN is seen most commonly in patients with liver disease or mucosal inflammation, in particular affecting the gastrointestinal tract), other renal diseases or autoimmune diseases. We also excluded control subjects with any chronic disease, conditions involving vital organs (liver, heart, lung and brain), central nervous system-related disease, or aggressive metabolic and endocrinological disease. Peripheral blood was collected from both cases and controls for DNA extraction.
All of the participants provided written informed consent. The Human Research Committee for Approval of Research Involving Human Subjects, the First Affiliated Hospital of Xi'an Jiaotong University, approved the use of human blood samples in this study.

SNP selection and genotyping
In this study, five SNPs in MPHOSPH6 were selected from DbSNP database (http://www.hapmap.org/ index.html.en) and SNP Consortium database (http:// snp.cshl.org/) for analysis. The lower frequency alleles were coded as the minor allele. All of the SNPs had minor allele frequencies (MAFs) > 5% in the HapMap Chinese Han Beijing population. Genomic DNA was isolated from whole blood samples using the GoldMag-Mini Purification Kit (GoldMagCo. Ltd. Xi'an, China), and DNA concentrations were measured using the NanoDrop2000 (Thermo Scientific, Waltham, MA, USA). Sequenom Massarray Assay Design 3.0 softwarewas used to design a multiplexed SNP Mass EXTENDED assay [17][18][19]. Genotyping was performed on a Sequenom MassARRAY RS1000 platform using the manufacturer's protocol. The PCR primers for each SNP are shown in Table 5. Data management and analysis was performed using the Sequenom Typer 4.0 Software [20,21].

Statistical analysis
We used Microsoft Excel (Microsoft Corporation, Redmond, WA, USA) and the SPSS 18.0 statistical package (SPSS, Chicago, IL, USA) to perform statistical analyses. All p values presented in this study were two sided, and p = 0.05 was considered the cutoff for statistical significance. Differences in the characteristics of the case and control study populations were analyzed using chi-square tests for categorical variables and Welch's t tests for continuous variables. In all analyses, the lower frequency allele was considered to be the 'risk' allele. Control genotype frequencies for each SNP were tested for departure from HWE using Fisher's exact tests. Allele and genotype frequencies in the cases and controls were compared using chi-square tests [22]. Four models (codominant, dominant, recessive, and log-additive) were used to assess the association between each genotype and the risk of IgAN. The effects of the polymorphisms on the risk of IgAN were expressed as odds ratios (ORs) with 95% confidence interval (CIs), which were calculated using unconditional logistic regression analysis after adjusting for age and gender [23].

CONCLUSIONS
The present study provided evidence that SNPs in the MPHOSPH6 are associated with IgAN in a Chinese Han population. It is possible that these variants are IgAN risk factors and these data can provide a theoretical foundation for other researchers to further study the association between the MPHOSPH6 gene and IgAN risk in the Chinese Han or other populations.