Discovery of Pathogenic Variants Associated with Idiopathic Recurrent Pregnancy Loss Using Whole-Exome Sequencing.

Idiopathic recurrent pregnancy loss (RPL) is defined as at least two pregnancy losses before 20 weeks of gestation. Approximately 5% of pregnant couples experience idiopathic RPL, which is a heterogeneous disease with various causes including hormonal, chromosomal, and intrauterine abnormalities. Although how pregnancy loss occurs is still unknown, numerous biological factors are associated with the incidence of pregnancy loss, including genetic variants. Whole-exome sequencing (WES) was conducted on blood samples from 56 Korean patients with RPL and 40 healthy controls. The WES data were aligned by means of bioinformatic analysis, and the detected variants were annotated using machine learning tools to predict the pathogenicity of protein alterations. Each indicated variant was confirmed using Sanger sequencing. A replication study was also conducted in 112 patients and 114 controls. The Variant Effect Scoring Tool, Combined Annotation Dependent Depletion tool, Sorting Intolerant from Tolerant annotation tool, and various databases detected 10 potential variants previously associated with spontaneous abortion genes in patients by means of a bioinformatic analysis of WES data. Several variants were detected in more than one patient. Interestingly, several of the detected genes were functionally clustered, including some with a secretory function (mucin 4; MUC4; rs200737893 G>A and hyaluronan-binding protein 2; HABP2; rs542838125 G>T), in which growth arrest-specific 2 Like 2 (GAS2L2; rs140842796 C>T) and dynamin 2 (DNM2; rs763894364 G>A) are functionally associated with cell protrusion and the cytoskeleton. ATP Binding Cassette Subfamily C Member 6 (ABCC6) was the only gene with two variants. HABP2 (rs542838125 G>T), MUC4 (rs200737893 G>A), and GAS2L2 (rs140842796 C>T) were detected in only the patient group in the replication study. The combination of WES and machine learning tools is a useful method to detect potential variants associated with RPL. Using bioinformatic tools, we found 10 potential variants in 9 genes. WES data from patients are needed to better understand the causes of RPL.

The Effects of Fermented Rice Drink With Lactiplantibacillus plantarum JSA22 in Overweight Irritable Bowel Syndrome Patients: A Randomized, Double-blind, Placebo-controlled Study.

Background/Aims: This study aims to investigate the effect of a fermented rice drink with Lactiplantibacillus plantarum JSA22 on symptoms, blood tests, microbiomes, and fecal metabolites in patients with irritable bowel syndrome (IBS) who were overweight. Methods: Sixty overweight (body mass index ≥ 23 kg/m2) patients aged between 20 and 65 with IBS were enrolled. Patients were divided into 2 groups and administered either a fermented rice drink or an nonfermented rice drink for a month. The symptom questionnaire, blood samples, and stool samples for microbiome and metabolite were collected before and after the month of rice drink administration. The primary efficacy variable was the subject's global assessment of IBS symptoms. Results: In both groups, global IBS symptoms, including abdominal pain, bowel habit, urgency, and abdominal distension, improved significantly (P < 0.01). The abdominal bloating was more significantly improved in the fermented rice drink group than in the nonfermented rice drink group (P < 0.05). Significant changes were not observed in metabolic syndrome-related blood tests or fecal metabolites in either group. However, microbiome analysis showed significant differences in genus levels before and after consuming fermented rice drink, such as in Blautia in stool (P = 0.020) and Prevotella (P = 0.017) and Oribacterium (P = 0.018) in saliva. Conclusions: The fermented rice drink with L. plantarum JSA22 showed a beneficial effect in reducing abdominal distension in IBS patients. Bacteria that reduce visceral fat accumulation increased in the stool and saliva of patients who consumed fermented rice drinks.

Safety tests and clinical research on buccal and nasal microneedle swabs for genomic analysis.

Conventional swabs have been used as a non-invasive method to obtain samples for DNA analysis from the buccal and the nasal mucosa. However, swabs may not always collect pure enough genetic material. In this study, buccal and nasal microneedle swab is developed to improve the accuracy and reliability of genomic analysis. A cytotoxicity test, a skin sensitivity test, and a skin irritation test are conducted with microneedle swabs. Polymer microneedle swabs meet the safety requirements for clinical research and commercial use. When buccal and nasal microneedle swabs are used, the amount of genetic material obtained is greater than that from commercially available swabs, and DNA purity is also high. The comparatively short microneedle swab (250 µm long) cause almost no pain to all 25 participants. All participants also report that the microneedle swabs are very easy to use. When genotypes are compared at five SNP loci from blood of a participant and from that person's buccal or nasal microneedle swab, the buccal and nasal microneedle swabs show 100% concordance for all five SNP genotypes. Microneedle swabs can be effectively used for genomic analysis and prevention through genomic analysis, so the utilization of microneedle swabs is expected to be high.

Corrigendum: Development of a microneedle swab for acquisition of genomic DNA from buccal cells.

[This corrects the article DOI: 10.3389/fbioe.2022.829648.].

Differences in the gut microbiome composition of Korean children and adult samples based on different DNA isolation kits.

Recent studies have revealed that the composition of human gut microbiota varies according to region, race, age, diet, living environment, and sampling and DNA extraction method. The purpose of this study was to broaden our understanding of the intestinal microbial composition of Koreans by conducting a 16S rRNA amplicon sequencing on 78 Korean samples composed of adults, children, normal and obese groups. We compared the microbiome composition and diversity of these groups at different levels including the phylum and genus level using two different stool DNA extraction kits of QIAamp® PowerFecal® DNA Kit (Qiagen, Hilden, Germany) and CT Max Fecal DNA Kit (Ct bio, Korea). We found that Ct bio (Ct) kit recovered higher DNA yields and OTUs than QIAamp® PowerFecal® DNA Kit (Qia). The Ct kit, which adopted more rigorous bead beating method, detected the most Gram-positive (G+) bacteria, Firmicutes, at the Phylum level, whereas the Qia kit, which used a less rigorous cell lysis method, found the most Gram-negative (G-) bacteria, Bacteroidetes. The Firmicutes-to-Bacteroidetes (F/B) ratio showed no significant difference between the obese and the normal groups of same kit; however, they were significantly different with two different kits. There was a difference in the intestinal flora between healthy Korean adults and children. The taxa that differed significantly between the adults and children were Bacteroides, Bifidobacterium, Prevotella, and Subdoligranulum. There was no significant difference in the intestinal flora between the normal weight group and the obese group in adults and children, respectively. This is probably because the difference in body mass index (BMI) between the sample groups collected in this study is statistically significant, but it is not large enough to show a clear difference in the flora. Therefore, these results should be interpreted with caution while considering the BMI values and Korean obesity criterion together.

Development of a Microneedle Swab for Acquisition of Genomic DNA From Buccal Cells.

A swab is a tool for obtaining buccal DNA from buccal mucus for biological analysis. The acquisition of a sufficient amount and high quality of DNA is an important factor in determining the accuracy of a diagnosis. A microneedle swab (MN swab) was developed to obtain more oral mucosal tissues non-invasively. Eight types of MN swabs were prepared with varying combinations of patterns (zigzag or straight), number of MNs, intervals of MNs, and sharpness of tips. When MN swab was applied up to 10 times, the tissue amount and DNA yield increased compared to commercial swabs. A zigzag pattern of microneedles was found to be more efficient than a straight pattern and increasing the number of microneedles in an array increased the DNA yield. The MN swab collected about twice the DNA compared to the commercial swab. In an in vivo test using mini pigs, the lower cycle threshold values of mucosal samples collected with MN swabs compared to samples collected with commercial swabs indicated that a greater amount of DNA was collected for SNP genotyping. A polymer MN swab is easy to manufacture by a single molding process, and it has a greater sampling capacity than existing commercial swabs.

Prospective clinical evaluation of Momguard non-invasive prenatal test in 1011 Korean high-risk pregnant women.

Clinical performance of the Momguard non-invasive prenatal test (NIPT) was evaluated in a cohort of Korean pregnant women. The foetal trisomies 21, 18 and 13 (T21, T18 and T13) were screened by low-coverage massive parallel sequencing in the maternal blood. Among the 1011 confirmed samples, 32 cases (3.2%) had positive NIPT results. Of these positive cases, 20 cases of T21, all cases of T18 and two cases of T13 had concordant karyotype findings. Only one case out of the remaining 979 negative NIPT samples showed a false negative result. The overall sensitivity and specificity of Momguard to detect the three chromosomal aneuploidies were 96.8% and 99.8%, respectively. Momguard is a clinically useful tool for the detection of T21, T18 and T13 in singleton pregnancy. However, as other NIPT tests, it carries the risk of false positive and false negative results. Hence, the genetic counsellors should provide these limitations to the examinees.Impact StatementWhat is already known on this subject? The NIPT approach using massive parallel sequencing (MPS) showed high sensitivity and specificity in various clinical studies. These results are based on analysis systems using their own bioinformatics algorithms.What the results of this study add? When this NIPT technology was introduced in Korea, the first biological specimens collected in Korea were transported overseas for processing in overseas laboratories and analysed by other country's analysis methods. We needed our own NIPT algorithm and developed Momguard NIPT for the first time in Korea. This study attempted to evaluate this Momguard NIPT protocol prospectively in a large number of samples obtained from three Korean hospitals.What the implications are of these findings for clinical practice and/or further research? The overall sensitivity and specificity to identify T13, T18 and T21 were 96.8% and 99.8%, respectively. These accuracy values were comparable to that of other studies. From this study, we found that Momguard is a clinically useful tool for the detection of three chromosomal aneuploidies. However, as other NIPT tests, it carries the risk of false positive and false negative results. Hence, the genetic counsellors should provide these limitations to the examinees.

Obesity Alters the Microbial Community Profile in Korean Adolescents.

Obesity is an increasing public health concern worldwide. According to the latest Organization for Economic Co-operation and Development (OECD) report (2014), the incidence of child obesity in Korea has exceeded the OECD average. To better understand and control this condition, the present study examined the composition of the gut microbial community in normal and obese adolescents. Fecal samples were collected from 67 obese (body mass index [BMI] ≥ 30 kg/m2, or ≥ 99th BMI percentile) and 67 normal (BMI < 25 kg/m2 or < 85th BMI percentile) Korean adolescents aged 13-16 years and subjected to 16S rRNA gene sequencing. Analysis of bacterial composition according to taxonomic rank (genus, family, and phylum) revealed marked differences in the Bacteroides and Prevotella populations in normal and obese samples (p < 0.005) at the genus and family levels; however, there was no difference in the Firmicutes-to-Bacteroidetes (F/B) ratio between normal and obese adolescents samples at the phylum level (F/B normal = 0.50 ± 0.53; F/B obese = 0.56 ± 0.86; p = 0.384). Statistical analysis revealed a significant association between the compositions of several bacterial taxa and child obesity. Among these, Bacteroides and Prevotella showed the most significant association with BMI (p < 0.0001 and 0.0001, respectively). We also found that the composition of Bacteroides was negatively associated with triglycerides (TG), total cholesterol, and high-sensitive C-reactive protein (hs-crp) (p = 0.0049, 0.0023, and 0.0038, respectively) levels, whereas that of Prevotella was positively associated with TG and hs-crp levels (p = 0.0394 and 0.0150, respectively). We then applied the association rule mining algorithm to generate "rules" to identify the association between the populations of multiple bacterial taxa and obesity; these rules were able to discriminate obese from normal states. Therefore, the present study describes a systemic approach to identify the association between bacterial populations in the gut and childhood obesity.

Polymorphisms in PDE4D are associated with a risk of COPD in non-emphysematous Koreans.

Despite extensive effort, only a few chronic obstructive pulmonary disease (COPD)-associated genes have been suggested, indicating that there must be additional risk-associated loci. Here we aimed to identify additional COPD-associated SNPs and to explore the potential relationship between COPD subgroups and the SNPs in the Korean population. We performed a genome-wide association study (GWAS) with 990 Korean individuals; 102 COPD cases and 544 controls for GWAS using Affymetrix SNP array 5.0, and 173 COPD cases and 171 controls for replication. After validating the candidate single nucleotide polymorphisms (SNP), we performed subgroup analysis by disease phenotype. Through GWAS, we identified a novel SNP in the phosphodiesterase-4D (PDE4D) gene [rs16878037 (C>T), p = 1.66 ◊ 10(-6)] that was significantly associated with COPD. This signal in PDE4D was successfully replicated in the independent set (p = 0.041). When we combined the discovery and replication data, the association signal became more significant (p = 5.69 ◊ 10(-7)). In the COPD subgroup analysis, the T allele of rs16878037 was significantly more frequent in COPD patients without severe diffusion capacity impairment (mild mixed and obstruction-dominant group) than in patients with severe impairment (severe mixed and emphysema-dominant groups). This result supports that PDE4D polymorphisms might be involved in the susceptibility to COPD especially in non-emphysematous individuals and that they could also affect the responsiveness of the PDE4 inhibitor treatment.

Genome-wide copy number variation analysis identifies deletion variants associated with ankylosing spondylitis.

OBJECTIVE: To identify ankylosing spondylitis (AS)-associated copy number variations (CNVs) in Korean subjects and their synergistic roles in the development of AS. METHODS: A genome-wide association study (GWAS) was performed in 309 patients with AS and 309 control subjects, using a copy number variant (CNV) microarray. AS-associated CNV regions were replicated in 2 independent sets (625 patients and 891 control subjects) by quantitative polymerase chain reaction (PCR) and deletion-typing PCR. RESULTS: In the CNV GWAS, 227 CNV regions were shown to be significantly associated with the risk of AS. Of the candidate CNV regions, 9 were successfully replicated in the first replication analysis: 1q32.2 (HHAT), 1p34.2 (BMP8A), 2q31.2 (PRKRA), 6p21.32 (HLA-DPB1), 11q22.1 (CNTN5), 13q13.1 (EEF1DP3), 14q24.2 (RGS6), 16p13.3, and 22q11.1 (IL17RA). The 5 deletion-type CNV regions, in 1q32.2, 2q31.2, 6p21.32, 13q13.1, and 16p13.3, were associated with an increased risk of AS, and the other 4 CNV regions were protective. In the second replication analysis, 4 CNV regions in 1q32.2, 2q31.2, 6p21.32, and 16p13.3 were replicated. Among patients with CNV regions in ≥4 risk-increasing loci, the risk was 18.0 times higher than that in patients without any deletions (odds ratio [OR] 17.98, P = 2.3 × 10(-7) ). Among patients with CNV regions in ≥2 protective loci, the risk was 5.2 times lower than that in those without any deletions (OR 0.19, P = 4.0 × 10(-10) ). The additive effects of simultaneous events were shown to be dependent on the frequency of CNV regions. Through deletion-typing PCR and sequencing, the exact sizes and breakpoint sequences were defined in 4 CNV regions. The mechanism of all 3 deletions was shown to be microhomology-based nonhomologous end joining. CONCLUSION: The results of this study can help to identify pathogenic mechanisms of AS and can easily be applied in the development of algorithms estimating the risk of AS.

New variants including ARG1 polymorphisms associated with C-reactive protein levels identified by genome-wide association and pathway analysis.

C-reactive protein (CRP) is a general marker of systemic inflammation and cardiovascular disease (CVD). The genetic contribution to differences in CRP levels remains to be explained, especially in non-European populations. Thus, the aim of this study was to identify genetic loci associated with CRP levels in Korean population. We performed genome-wide association studies (GWAS) using SNPs from 8,529 Korean individuals (7,626 for stage 1 and 903 for stage 2). We also performed pathway analysis. We identified a new genetic locus associated with CRP levels upstream of ARG1 gene (top significant SNP: rs9375813, Pmeta = 2.85×10(-8)), which encodes a key enzyme of the urea cycle counteract the effects of nitric oxide, in addition to known CRP (rs7553007, Pmeta = 1.72×10(-16)) and HNF1A loci (rs2259816, Pmeta = 2.90×10(-10)). When we evaluated the associations between the CRP-related SNPs with cardiovascular disease phenotypes, rs9375813 (ARG1) showed a marginal association with hypertension (P = 0.0440). To identify more variants and pathways, we performed pathway analysis and identified six candidate pathways comprised of genes related to inflammatory processes and CVDs (CRP, HNF1A, PCSK6, CD36, and ABCA1). In addition to the previously reported loci (CRP, HNF1A, and IL6) in diverse ethnic groups, we identified novel variants in the ARG1 locus associated with CRP levels in Korean population and a number of interesting genes related to inflammatory processes and CVD through pathway analysis.

Genome-wide association study of lung cancer in Korean non-smoking women.

Lung cancer in never-smokers ranks as the seventh most common cause of cancer death worldwide, and the incidence of lung cancer in non-smoking Korean women appears to be steadily increasing. To identify the effect of genetic polymorphisms on lung cancer risk in non-smoking Korean women, we conducted a genome-wide association study of Korean female non-smokers with lung cancer. We analyzed 440,794 genotype data of 285 cases and 1,455 controls, and nineteen SNPs were associated with lung cancer development (P < 0.001). For external validation, nineteen SNPs were replicated in another sample set composed of 293 cases and 495 controls, and only rs10187911 on 2p16.3 was significantly associated with lung cancer development (dominant model, OR of TG or GG, 1.58, P = 0.025). We confirmed this SNP again in another replication set composed of 546 cases and 744 controls (recessive model, OR of GG, 1.32, P = 0.027). OR and P value in combined set were 1.37 and < 0.001 in additive model, 1.51 and < 0.001 in dominant model, and 1.54 and < 0.001 in recessive model. The effect of this SNP was found to be consistent only in adenocarcinoma patients (1.36 and < 0.001 in additive model, 1.49 and < 0.001 in dominant model, and 1.54 and < 0.001 in recessive model). Furthermore, after imputation with HapMap data, we found regional significance near rs10187911, and five SNPs showed P value less than that of rs10187911 (rs12478012, rs4377361, rs13005521, rs12475464, and rs7564130). Therefore, we concluded that a region on chromosome 2 is significantly associated with lung cancer risk in Korean non-smoking women.

Deletion variants of RABGAP1L, 10q21.3, and C4 are associated with the risk of systemic lupus erythematosus in Korean women.

OBJECTIVE: Several copy number variations (CNVs) have been found to be associated with systemic lupus erythematosus (SLE) through the target gene approach. However, genome-wide features of CNVs and their role in the risk of SLE remain unknown. The aim of this study was to identify SLE-associated CNVs in Korean women. METHODS: Genome-wide assessments of CNVs were performed in 382 SLE patients and 191 control subjects, using an Illumina HumanHap610 BeadChip genotyping platform. SLE-associated CNV regions that were identified by genome-wide association study (GWAS) were replicated in quantitative polymerase chain reaction (PCR) and deletion-typing PCR analyses in an independent sample set comprising 564 SLE patients and 511 control subjects. RESULTS: Of 144 common CNV regions, 3 deletion-type CNV regions in 1q25.1, 8q23.3, and 10q21.3 were found to be significantly associated with SLE by GWAS analysis. In the independent replication, the CNV regions in 1q25.1 (RABGAP1L) and 10q21.3 were successfully replicated (odds ratio [OR] 1.30, P=0.038 and OR 1.90, P=3.6×10(-5), respectively), and the associations were confirmed again by deletion-typing PCR. The CNV region in the C4 gene, which showed a potential association in the discovery stage, was included in the replication analysis and was found to be significantly associated with the risk of SLE (OR 1.88, P=0.01). Through deletion-typing PCR, the exact sizes and breakpoint sequences of the deletions were defined. Individuals with the deletions in all 3 loci (RABGAP1L, 10q21.3, and C4) had a much higher risk of SLE than did those without any deletions in the 3 loci (OR 5.52, P=3.9×10(-4)). CONCLUSION: These CNV regions can be useful to identify the pathogenic mechanisms of SLE, and might be used to more accurately predict the risk of SLE by taking into consideration their synergistic effects on disease susceptibility.

Web-based database and viewer of East asian copy number variations.

We have discovered copy number variations (CNVs) in 3,578 Korean individuals with the Affymetrix Genome-Wide SNP array 5.0, and 4,003 copy number variation regions (CNVRs) were defined in a previous study. To explore the details of the variants easily in related studies, we built a database, cataloging the CNVs and related information. This system helps researchers browsing these variants with gene and structure variant annotations. Users can easily find specific regions with search options and verify them from system-integrated genome browsers with annotations.

Identifying Copy Number Variants under Selection in Geographically Structured Populations Based on F-statistics.

Large-scale copy number variants (CNVs) in the human provide the raw material for delineating population differences, as natural selection may have affected at least some of the CNVs thus far discovered. Although the examination of relatively large numbers of specific ethnic groups has recently started in regard to inter-ethnic group differences in CNVs, identifying and understanding particular instances of natural selection have not been performed. The traditional F(ST) measure, obtained from differences in allele frequencies between populations, has been used to identify CNVs loci subject to geographically varying selection. Here, we review advances and the application of multinomial-Dirichlet likelihood methods of inference for identifying genome regions that have been subject to natural selection with the F(ST) estimates. The contents of presentation are not new; however, this review clarifies how the application of the methods to CNV data, which remains largely unexplored, is possible. A hierarchical Bayesian method, which is implemented via Markov Chain Monte Carlo, estimates locus-specific F(ST) and can identify outlying CNVs loci with large values of F(ST). By applying this Bayesian method to the publicly available CNV data, we identified the CNV loci that show signals of natural selection, which may elucidate the genetic basis of human disease and diversity.

Simultaneous copy number gains of NUPR1 and ERBB2 predicting poor prognosis in early-stage breast cancer.

BACKGROUND: The full extent of chromosomal alterations and their biological implications in early breast carcinogenesis has not been well examined. In this study, we aimed to identify chromosomal alterations associated with poor prognosis in early-stage breast cancers (EBC). METHODS: A total of 145 EBCs (stage I and II) were examined in this study. We analyzed copy number alterations in a discovery set of 48 EBCs using oligoarray-comparative genomic hybridization. In addition, the recurrently altered regions (RARs) associated with poor prognosis were validated using an independent set of 97 EBCs. RESULTS: A total of 23 RARs were defined in the discovery set. Six were commonly detected in both stage I and II groups (> 50%), suggesting their connection with early breast tumorigenesis. There were gains on 1q21.2-q21.3, 8q24.13, 8q24.13-21, 8q24.3, and 8q24.3 and a loss on 8p23.1-p22. Among the 23 RARs, copy number gains on 16p11.2 (NUPR1) and 17q12 (ERBB2) showed a significant association with poor survival (P = 0.0186 and P = 0.0186, respectively). The patients simultaneously positive for both gains had a significantly worse prognosis (P = 0.0001). In the independent replication, the patients who were double-positive for NUPR1-ERBB2 gains also had a significantly poorer prognosis on multivariate analysis (HR = 7.31, 95% CI 2.65-20.15, P = 0.0001). CONCLUSIONS: The simultaneous gain of NUPR1 and ERBB2 can be a significant predictor of poor prognosis in EBC. Our study will help to elucidate the molecular mechanisms underlying early-stage breast cancer tumorigenesis. This study also highlights the potential for using combinations of copy number alterations as prognosis predictors for EBC.

CNVRuler: a copy number variation-based case-control association analysis tool.

SUMMARY: The method for genome-wide association study (GWAS) based on copy number variation (CNV) is not as well established as that for single nucleotide polymorphism (SNP)-GWAS. Although there are several tools for CNV association studies, most of them do not provide appropriate definitions of CNV regions (CNVRs), which are essential for CNV-association studies. Here we present a user-friendly program called CNVRuler for CNV-association studies. Outputs from the 10 most common CNV defining algorithms can be directly used as input files for determining the three different definitions of CNVRs. Once CNVRs are defined, CNVRuler supports four kinds of statistical association tests and options for population stratification. CNVRuler is based on the open-source programs R and Java from Sun Microsystems. AVAILABILITY: CNVRuler software is available with an online manual at the website, www.ircgp.com/CNVRuler/index.html.

Common variants at the promoter region of the APOM confer a risk of rheumatoid arthritis.

Although the genetic component in the etiology of rheumatoid arthritis (RA) has been consistently suggested, many novel genetic loci remain to uncover. To identify RA risk loci, we performed a genome-wide association study (GWAS) with 100 RA cases and 600 controls using Affymetrix SNP array 5.0. The candidate risk locus (APOM gene) was re-sequenced to discover novel promoter and coding variants in a group of the subjects. Replication was performed with the independent case-control set comprising of 578 RAs and 711 controls. Through GWAS, we identified a novel SNP associated with RA at the APOM gene in the MHC class III region on 6p21.33 (rs805297, odds ratio (OR) = 2.28, P = 5.20 × 10-7). Three more polymorphisms were identified at the promoter region of the APOM by the re-sequencing. For the replication, we genotyped the four SNP loci in the independent case-control set. The association of rs805297 identified by GWAS was successfully replicated (OR = 1.40, P = 6.65 × 10-5). The association became more significant in the combined analysis of discovery and replication sets (OR = 1.56, P = 2.73 × 10-10). The individuals with the rs805297 risk allele (A) at the promoter region showed a significantly lower level of APOM expression compared with those with the protective allele (C) homozygote. In the logistic regressions by the phenotype status, the homozygote risk genotype (A/A) consistently showed higher ORs than the heterozygote one (A/C) for the phenotype-positive RAs. These results indicate that APOM promoter polymorphisms are significantly associated with the susceptibility to RA.