Mitochondrial haplogroups in patients with rheumatoid arthritis: No association with disease and disease manifestations

Differential mitochondrial genome in patients with Rheumatoid Arthritis

BACKGROUND

Mitochondria play an important role in cell survival, function and lineage differentiation. Changes in mitochondrial DNA (mtDNA) may control mitochondrial functions and thus may impart an alternative cellular state thereby leading to a disease condition in the body. Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease wherein immune cells become self-reactive causing joint inflammation, swelling and pain in patients. The changes in mtDNA may alter cellular functions thereby directing the immune cells towards an inflammatory phenotype in RA. Therefore, it becomes pertinent to identify changes in mtDNA sequence in immune cells of RA patients to understand the pathogenesis and progression of RA.

METHODS

mtDNA from peripheral blood mono-nuclear cells (PBMCs) of 23 RA patients and 17 healthy controls (HCs) were sequenced using next-generation sequencing (NGS). Further, single nucleotide polymorphisms (SNPs) and other variable changes in mtDNA hypervariable and coding regions, amino acid changes with a putative impact on disease, levels of heteroplasmy, copy number variations and haplogroup analysis in RA patients and HCs were analysed and compared to identify any association of mtDNA changes and RA disease.

RESULTS

A total of 382 single nucleotide mtDNA variants were observed, 91 (23.82%) were present in hypervariable region and 291 (76.18%) in coding region of patients and HC. The variant 513 GCA > ACA, with G present in HVR-III, known to control the mitochondrial translation function, was significantly present in RA patients. The CYTB gene had larger number of SNPs in HC samples while RNR2 was more variable in RA patients. A non-synonymous heteroplasmy in ND1 gene was found at a single nucleotide position 3533 in an increased number of RA patients as compared to the controls. A significant increase in mtDNA duplication and a higher frequency of the haplogroup U was also characteristic of RA. Also, the presence of SNPs in mitochondrial tRNA genes at two positions 12308 A > G and 15924 A > G were found to be pathogenic.

CONCLUSION

We herein observed an altered mtDNA sequence in immune cells of RA patients and thus a possible role of mitochondrial genome in the development of RA. The observed nucleotide changes in mtDNA control region, RNR2 gene, increased heteroplasmy and mtDNA duplication in RA patients may alter sites for transcription factor binding thereby influencing mtDNA gene expression, as well as copy numbers thereby affecting the mitochondrial proteins and their functions. These changes in mtDNA could be one of the probable reasons among many leading to the progression of RA.

Germline and somatic mtDNA mutation spectrum of rheumatoid arthritis patients in the Taizhou area, China

Objective

Reactive oxygen species are believed to be involved in the onset of RA, and the association between nuclear-encoded mitochondrial respiratory chain-related variants and RA has recently been revealed. However, little is known about the landscape of mitochondrial DNA (mtDNA) variants in RA.

Methods

Next-generation sequencing was conducted to profile mtDNA germline and somatic variants in 124 RA patients and 123 age- and sex-matched healthy controls in the Taizhou area, China. Fisher’s exact test, SKAT and SKAT-O were used for gene-burden tests to investigate RA-related variants of mitochondrial genes. Predictive tools were applied to evaluate the pathogenicity of mtDNA variants, and mtDNA haplogroups were assigned according to mtDNA mutations recorded in PhyloTree database. The frequency distribution of mtDNA haplogroups between the groups was compared using χ2 analysis.

Results

We identified 467 RA-unique and 341 healthy control-unique mtDNA variants, with 443 common variants. Only MT-ATP6 with a significant burden of variants was identified by Fisher’s exact test, SKAT and SKAT-O, even after Bonferroni adjustment, and the enrichment variants in MT-ATP6 was mainly driven by m.8830C>A, m.8833G>C and m.8843T>A variants. Besides, four frequently low-heteroplasmic variants including the three variants above and m.14135T>G of MT-ND5 were detected in RA only; except for m.8830C>A, they are considered potential pathogenicity based on functional predictions. χ2 analysis before Bonferroni adjustment revealed haplogroup F1/F1a to be negatively associated with RA (P < 0.05).

Conclusion

These results profiled the landscape of germline and somatic mtDNA variants in RA and supported the effects of mitochondrial genes on RA.

Environmental factors modulated ancient mitochondrial DNA variability and the prevalence of rheumatic diseases in the Basque Country

Among the factors that would explain the distribution of mitochondrial lineages in Europe, climate and diseases may have played an important role. A possible explanation lies in the nature of the mitochondrion, in which the energy generation process produces reactive oxygen species that may influence the development of different diseases. The present study is focused on the medieval necropolis of San Miguel de Ereñozar (13^th-16^th centuries, Basque Country), whose inhabitants presented a high prevalence of rheumatic diseases and lived during the Little Ice Age (LIA). Our results indicate a close relationship between rheumatic diseases and mitochondrial haplogroup H, and specifically between spondyloarthropathies and sub-haplogroup H2. One possible explanation may be the climate change that took place in the LIA that favoured those haplogroups that were more energy-efficient, such as haplogroup H, to endure lower temperatures and food shortage. However, it had a biological trade-off: the increased risk of developing rheumatic diseases.