Simple regression for correcting ΔCt bias in RT-qPCR low-density array data normalization.

BACKGROUND: Reverse transcription quantitative PCR (RT-qPCR) is considered the gold standard for quantifying relative gene expression. Normalization of RT-qPCR data is commonly achieved by subtracting the Ct values of the internal reference genes from the Ct values of the target genes to obtain ΔCt. ΔCt values are then used to derive ΔΔCt when compared to a control group or to conduct further statistical analysis. RESULTS: We examined two rheumatoid arthritis RT-qPCR low density array datasets and found that this normalization method introduces substantial bias due to differences in PCR amplification efficiency among genes. This bias results in undesirable correlations between target genes and reference genes, which affect the estimation of fold changes and the tests for differentially expressed genes. Similar biases were also found in multiple public mRNA and miRNA RT-qPCR array datasets we analysed. We propose to regress the Ct values of the target genes onto those of the reference genes to obtain regression coefficients, which are then used to adjust the reference gene Ct values before calculating ΔCt. CONCLUSIONS: The per-gene regression method effectively removes the ΔCt bias. This method can be applied to both low density RT-qPCR arrays and individual RT-qPCR assays.

Gene expression patterns in peripheral blood cells associated with radiographic severity in African Americans with early rheumatoid arthritis.

Gene expression profiling may be used to stratify patients by disease severity to test the hypothesis that variable disease outcome has a genetic component. In order to define unique expression signatures in African American rheumatoid arthritis (RA) patients with severe erosive disease, we undertook a gene expression study using samples of RNA from peripheral blood mononuclear cells (PBMCs). RNA from baseline PBMC samples of 96 African American RA patients with early RA (<2 years disease duration) was hybridized to cDNA probes of the Illumina Human HT-V3 expression array. Expression analyses were performed using the ca. 25,000 cDNA probes, and then expression levels were compared to the total number of erosions in radiographs of the hands and feet at baseline and 36 months. Using a false discovery rate cutoff of Q = 0.30, 1,138 genes at baseline and 680 genes at 36 months significantly correlated with total erosions. No evidence of a signal differentiating disease progression, or change in erosion scores between baseline and 36 months, was found. Further analyses demonstrated that the differential gene expression signature was localized to the patients with the most erosive disease (>10 erosions). Ingenuity Pathway Analysis demonstrated that genes with fold change greater than 1.5 implicated immune pathways such as CTLA signaling in cytotoxic T lymphocytes. These results demonstrate that CLEAR patients with early RA having the most severe erosive disease, as compared to more mild cases (<10 erosions), may be characterized by a set of differentially expressed genes that represent biological pathways with relevance to autoimmune disease.

Association of single-nucleotide polymorphisms in CCR6, TAGAP, and TNFAIP3 with rheumatoid arthritis in African Americans.

OBJECTIVE: We previously reported an analysis of single-nucleotide polymorphisms (SNPs) in 3 validated European rheumatoid arthritis (RA) susceptibility loci, TAGAP, TNFAIP3, and CCR6, in African American patients with RA. Unexpectedly, the disease-associated alleles were different in African Americans from those in Europeans. In an effort to better define their contribution, we performed additional SNP genotyping in these genes. METHODS: Seven SNPs were genotyped in 446 African American patients with RA and in 733 African American control subjects. Differences in minor allele frequency between the RA cases and controls were analyzed after controlling for the global proportion of European admixture, and pairwise linkage disequilibrium (LD) was estimated among the SNPs. RESULTS: Three SNPs were significantly associated with RA: the TNFAIP3 rs719149 A allele (OR 1.22 [95% confidence interval (95% CI) 1.03-1.44], P = 0.02), the TAGAP rs1738074 G allele (OR 0.75 [95% CI 0.63-0.89, P = 0.0012), and the TAGAP rs4709267 G allele (OR 0.74 [95% CI 0.60-0.91], P = 0.004). Pairwise LD between the TAGAP SNPs was low (r(2) = 0.034). The haplotype containing minor alleles for both TAGAP SNPs was uncommon (4.5%). After conditional analysis of each TAGAP SNP, its counterpart remained significantly associated with RA (rs1738074 for rs4709267 P = 0.00001 and rs4709267 for rs1738074 P = 0.00005), suggesting independent effects. CONCLUSION: SNPs in regulatory regions of TAGAP and an intronic SNP (TNFAIP3) are potential susceptibility loci in African Americans. Pairwise LD, haplotype analysis, and SNP conditioning analysis suggest that these 2 SNPs in TAGAP are independent susceptibility alleles. Additional fine-mapping of this gene and functional genomic studies of these SNPs should provide further insight into the role of these genes in RA.

Radiographic severity of rheumatoid arthritis in African Americans: results from a multicenter observational study.

OBJECTIVE: To describe radiographic changes in African Americans with rheumatoid arthritis (RA) from the Consortium for the Longitudinal Evaluations of African Americans with Early Rheumatoid Arthritis (CLEAR) Registry, a multicenter observational study. METHODS: Self-declared African American patients were enrolled in CLEAR I, a longitudinal cohort of early RA (disease duration of <2 years) from 2000 to 2005, or in CLEAR II, a cross-sectional cohort (any disease duration) from 2006 to the present. Demographic and clinical data were obtained, and sets of hand/wrist and foot radiographs were scored using the modified Sharp/van der Heijde scoring system. RESULTS: A total of 357 and 418 patients were enrolled in CLEAR I and CLEAR II, respectively. We report here an interim analysis of radiographic severity in these patients. For the CLEAR I cohort, 294 patients had a mean radiographic score of 2.89 at the baseline visit; 32.0% showed either erosions (25.9%) or joint space narrowing (JSN; 19.4%). At the 36-month visit, the mean score was 5.65; 44.2% had erosions, 41.5% had JSN, and 54.4% had either. Among those patients without radiographic damage at baseline, 18.9% had progressed at the 36-month visit, compared with 57.1% of those with baseline damage (P < 0.0001). For the CLEAR II cohort, of 167 patients with RA of any duration, 65.3% exhibited joint erosions, 65.3% exhibited JSN, and 74.8% exhibited either. The mean radiographic score was 33.42. CONCLUSION: To our knowledge, this is the largest radiographic study of African American RA patients. Damage occurs early in the disease and is associated with radiographic progression at 3 years of disease duration. The CLEAR Registry will provide a valuable resource for future analyses of genetic, clinical, and environmental factors associated with radiographic severity of RA in African Americans.

Association of IL4R single-nucleotide polymorphisms with rheumatoid nodules in African Americans with rheumatoid arthritis.

INTRODUCTION: To determine whether IL4R single-nucleotide polymorphisms (SNPs) rs1805010 (I50V) and rs1801275 (Q551R), which have been associated with disease severity in rheumatoid arthritis (RA) patients of European ancestry, relate to the presence of rheumatoid nodules and radiographic erosions in African Americans. METHODS: Two IL4R SNPs, rs1805010 and rs1801275, were genotyped in 749 patients from the Consortium for Longitudinal Evaluation of African-Americans with Early Rheumatoid Arthritis (CLEAR) registries. End points were rheumatoid nodules defined as present either by physical examination or by chest radiography and radiographic erosions (radiographs of hands/wrists and feet were scored using the modified Sharp/van der Heijde system). Statistical analyses were performed by using logistic regression modeling adjusted for confounding factors. RESULTS: Of the 749 patients with RA, 156 (20.8%) had rheumatoid nodules, with a mean age of 47.0 years, 84.6% female gender, and median disease duration of 1.9 years. Of the 461 patients with available radiographic data, 185 (40.1%) had erosions (score>0); their mean age was 46.7 years; 83.3% were women; and median disease duration was 1.5 years. Patients positive for HLA-DRB1 shared epitope (SE) and autoantibodies (rheumatoid factor (RF) or anti-cyclic citrullinated peptide (CCP)) had a higher risk of developing rheumatoid nodules in the presence of the AA and AG alleles of rs1801275 (odds ratio (OR)adj=8.08 (95% confidence interval (CI): 1.60-40.89), P=0.01 and ORadj=2.97 (95% CI, 1.08 to 8.17), P=0.04, respectively). Likewise, patients positive for the HLA-DRB1 SE and RF alone had a higher risk of developing rheumatoid nodules in presence of the AA and AG alleles of rs1801275 (ORadj=8.45 (95% CI, 1.57 to 45.44), P=0.01, and ORadj=3.57 (95% CI, 1.18 to 10.76), P=0.02, respectively) and in the presence of AA allele of rs1805010 (ORadj=4.52 (95% CI, 1.20 to 17.03), P=0.03). No significant association was found between IL4R and radiographic erosions or disease susceptibility, although our statistical power was limited by relatively small numbers of cases and controls. CONCLUSIONS: We found that IL4R SNPs, rs1801275 and rs1805010, are associated with rheumatoid nodules in autoantibody-positive African-American RA patients with at least one HLA-DRB1 allele encoding the SE. These findings highlight the need for analysis of genetic factors associated with clinical RA phenotypes in different racial/ethnic populations.