A copy number variation in human NCF1 and its pseudogenes

Low C4A copy numbers and higher HERV gene insertion contributes to increased risk of SLE, with absence of association with disease phenotype and disease activity

Low copy numbers (CNs) of C4 genes are associated with systemic autoimmune disorders and affects autoantibody diversity and disease subgroups. The primary objective of this study was to characterize diversity of complement (C4) and C4-Human Endogenous Retrovirus (HERV) gene copy numbers in SLE. We also sought to assess the association of C4 and C4-HERV CNs with serum complement levels, autoantibodies, disease phenotypes and activity. Finally, we checked the association of C4 and HERV CNs with specific HLA alleles. Genomic DNA from 70 SLE and 90 healthy controls of south Indian Tamil origin were included. Demographic, clinical and serological data was collected in a predetermined proforma. CNs of C4A and C4B genes and the frequency of insertion of 6.4kb HERV within C4 gene (C4AL, C4BL) was determined using droplet digital polymerase chain reaction (ddPCR). A four digit high resolution HLA genotyping was done using next generation sequencing. In our cohort, the total C4 gene copies ranged from 2 to 6. Compared to controls, presence of two or less copies of C4A gene was associated with SLE risk (p = 0.005; OR = 2.79; 95% CI = 1.29-6.22). Higher frequency of HERV insertion in C4A than in C4B increases such risk (p = 0.000; OR = 12.67; 95% CI = 2.80-115.3). AL-AL-AL-BS genotype was significantly higher in controls than SLE (9%vs1%, p = 0.04; OR = 0.15, 95% CI = 0.00-0.16). Distribution of HLA alleles was not different in SLE compared to controls as well as in SLE subjects with ≤ 2 copies and > 2 copies of C4A, but HLA allele distribution was diverse in subjects with C4B ≤ 2 copies and > 2 copies. Finally, there was no correlation between the C4 and the C4-HERV diversity and complement levels, autoantibodies, disease phenotypes and activity. In conclusion, our data show that, low C4A copy number and higher insertion of HERV-K in C4A increases the risk for SLE. C4 and C4-HERV CNs did not correlate with serum complements, autoantibodies, disease phenotypes and activity in SLE. Further validation in a larger homogenous SLE cohort is needed.

Subcutaneous adipose tissue alteration in aging process associated with thyroid hormone signaling

BACKGROUND

Functional changes in subcutaneous adipose tissue (SAT) occur earlier in the aging process and play an important role in the occurrence and development of age-related metabolic diseases. The mechanism of this phenomenon is still unclear, and the change in adipose tissue with age is poorly understood.

METHODS

We used transcriptome sequencing (RNA seq) to screen differentially expressed genes at the mRNA level, and analyzed the functional characteristics of the differential genes through GO and KEGG analysis in human SAT of all ages. In order to clarify the specific mechanism of the functional change, we analyzed the chromatin accessibility in the promoter region in the same SAT used in the RNA seq by the assay for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) and obtained the functional genes in SAT changed with age. To verify these changes, we enlarged our sample content of human SAT. The primary mice adipocytes were extracted and stimulated by thyroid hormone of different concentration to construct an animal model, and the expression of the genes were determined through real-time Polymerase Chain Reaction(RT-PCR). The oxygen consumption test and immunofluorescence staining were used to determine the mitochondrial function of SAT.

RESULTS

RNA-seq showed characteristic gene expression of young and old human SAT, in which 331 genes were up-regulated and 349 genes were down-regulated. ATAC-seq, RNA-seq, combined with the mouse prediction model, determined the functional changed characteristics of seven genes. All these genes expressed differently in SAT of different ages, in which, NCF1, NLRP3, DUOX1 showed positive correlation with age; The expression of IFI30, P2RX1, P2RX6, PRODH, however, decreased with age. And all these genes showed dose dependent alternations under treatment of triiodothyroxine in mice SAT. The oxygen consumption rate revealed significant changes of the mitochondrial function and ROS accumulation in human SAT of different ages.

CONCLUSION

In elderly individuals, the function, in addition to distribution, of SAT undergoes significant changes, primarily in mitochondria, which may be due to insensitivity to thyroid hormone signaling. These results identified seven novel genes regulated by thyroid hormone, exhibiting significant changes in SAT of different age, and are probably related to the dysfunction of the aged SAT due to the mitochondrial damage and ROS accumulation.

Translational enhancement by base editing of the Kozak sequence rescues haploinsufficiency

A variety of single-gene human diseases are caused by haploinsufficiency, a genetic condition by which mutational inactivation of one allele leads to reduced protein levels and functional impairment. Translational enhancement of the spare allele could exert a therapeutic effect. Here we developed BOOST, a novel gene-editing approach to rescue haploinsufficiency loci by the change of specific single nucleotides in the Kozak sequence, which controls translation by regulating start codon recognition. We evaluated for translational strength 230 Kozak sequences of annotated human haploinsufficient genes and 4621 derived variants, which can be installed by base editing, by a high-throughput reporter assay. Of these variants, 149 increased the translation of 47 Kozak sequences, demonstrating that a substantial proportion of haploinsufficient genes are controlled by suboptimal Kozak sequences. Validation of 18 variants for 8 genes produced an average enhancement in an expression window compatible with the rescue of the genetic imbalance. Base editing of the NCF1 gene, whose monoallelic loss causes chronic granulomatous disease, resulted in the desired increase of NCF1 (p47phox) protein levels in a relevant cell model. We propose BOOST as a fine-tuned approach to modulate translation, applicable to the correction of dozens of haploinsufficient monogenic disorders independently of the causing mutation.

Transcriptomic Analysis of Mycobacterium leprae-Stimulated Response in Peripheral Blood Mononuclear Cells Reveal Potential Biomarkers for Early Diagnosis of Leprosy

We aimed to identify an unique host transcriptional signature in peripheral blood mononuclear cells (PBMCs) in response to Mycobacterium leprae antigens to distinguish between patients with leprosy and non-leprosy controls for early diagnosis of the disease. Sixteen individuals were enrolled in the discovery cohort [eight patients with leprosy, comprising four multibacillary (MB) and four paucibacillary (PB); and eight non-leprosy controls, comprising four healthy house contacts (HHCs) and four endemic controls (ECs)]. The differences in the transcriptome response of PBMCs to M. leprae sonicate antigen were evaluated between leprosy patients and non-leprosy controls, and 12 differentially expressed genes (CCL2/MCP-1, IL-8, JAKM, ATP, ND1, SERP, FLJ10489, LINC00659, LOC34487, LOC101928143, MIR22, and NCF1C) were identified. The accuracy of the 12 differentially expressed genes was further validated for the diagnosis of leprosy using real-time quantitative PCR in 82 individuals (13 MB, 10 PB, 37 HHCs, and 22 ECs) in the validation cohort. We found that a 5 gene signature set IL-8, CCL2/MCP-1, SERP, LINC00659 and FLJ10489 had a suitable performance in discriminating leprosy from ECs. In addition, elevated expression of IL-8, CCL2/MCP-1, SERP and LINC00659 was associated with MB diagnosis compared with ECs, whereas increased expression of IL-8, CCL2/MCP-1, SERP and FLJ10489 was found to be useful biomarkers for PB diagnosis from ECs. Moreover, we found decreased expression of NCF1C among leprosy patients could distinguish leprosy from HHCs, whereas higher expression of CCL2 among MB than PB could distinguish different leprosy patients. In conclusion, among the 12 candidate host genes identified, a three gene signature IL-8, CCL2/MCP-1, and SERP showed the best performance in distinguishing leprosy patients from healthy controls. These findings may have implications for developing a rapid blood-based test for early diagnosis of leprosy.

Neutrophil Cytosolic Factor-1 Genotyping in Acne Vulgaris

Acne vulgaris (AV) is a very common inflammatory dermatosis. It has a complex pathogenesis in which oxidative stress plays an important role. Neutrophil cytosolic factor (NCF)-1 gene encodes for NCF1 protein which shares in reactive oxygen species (ROS) production. Copy number variation (CNV) is a type of genetic variance in which gene copies are duplicated or deleted. The current work aimed to detect the association between NCF1 CNV and NCF-1 genotypes and AV to explore their possible role in increased disease risk or influencing its clinical presentation. Twenty-five cases with AV and 25 age- and gender-matched healthy volunteers were selected. NCF1 CNV and genotypes were determined using quantitative real-time polymerase chain reaction. NCF1 copy number was significantly increased in patients compared to the control group (p = 0.02). Higher copy number increased the risk of occurrence of AV by about 4-fold. The NCF1 genotype was more prevalent in patients (72%) compared to NCF1B (24%) and NCF1C (4%) variants, while NCF1B and NCF1C variants (68%) were more prevalent in the control group. The NCF1B genotype decreased the risk of occurrence of AV by 0.2-fold. NCF1 was significantly associated with cases more than controls (p = 0.005). It increased the risk of occurrence of acne by 5.4-fold. There was significant association between NCF1 copy number and disease duration where higher number was associated with long disease duration (p = 0.03). Higher copy number was also associated with the NCF1 genotype (p = 0.01). This study suggests that increased copy number of NCF1 gene may be a predisposing factor for AV development. However, the presence of NCF1B and NCF1C variants lowers ROS production and subsequently decreases the risk of development of AV.

Differential Transgene Silencing of Myeloid-Specific Promoters in the AAVS1 Safe Harbor Locus of Induced Pluripotent Stem Cell-Derived Myeloid Cells

Targeted integration into a genomic safe harbor, such as the AAVS1 locus on chromosome 19, promises predictable transgene expression and reduces the risk of insertional mutagenesis in the host genome. The application of gamma-retroviral long terminal repeat (LTR)-driven vectors, which semirandomly integrate into the genome, has previously caused severe adverse events in some clinical studies due to transactivation of neighboring proto-oncogenes. Consequently, the site-specific integration of a therapeutic transgene into a genomic safe harbor locus would allow stable genetic correction with a reduced risk of insertional mutagenesis. However, recent studies revealed that transgene silencing, especially in case of weaker cell type-specific promoters, can occur in the AAVS1 locus of human pluripotent stem cells (PSCs) and can impede transgene expression during differentiation. In this study, we aimed to correct p47^phox deficiency, which is the second most common cause of chronic granulomatous disease, by insertion of a therapeutic p47^phox transgene into the AAVS1 locus of human induced PSCs (iPSCs) using CRISPR-Cas9. We analyzed transgene expression and functional correction from three different myeloid-specific promoters (miR223, CatG/cFes, and myeloid-related protein 8 [MRP8]). Upon myeloid differentiation of corrected iPSC clones, we observed that the miR223 and CatG/cFes promoters achieved therapeutically relevant levels of p47^phox expression and nicotinamide adenine dinucleotide phosphate oxidase activity, whereas the MRP8 promoter was less efficient. Analysis of the different promoters revealed high CpG methylation of the MRP8 promoter in differentiated cells, which correlated with the transgene expression data. In summary, we identified the miR223 and CatG/cFes promoters as cell type-specific promoters that allow stable transgene expression in the AAVS1 locus of iPSC-derived myeloid cells. Our findings further indicate that promoter silencing can occur in the AAVS1 safe harbor locus in differentiated hematopoietic cells and that a comparison of different promoters is necessary to achieve optimal transgene expression for therapeutic application of iPSC-derived cells.

NCF1 (p47phox)-deficient chronic granulomatous disease: comprehensive genetic and flow cytometric analysis

Mutations in NCF1 (p47phox) cause autosomal recessive chronic granulomatous disease (CGD) with abnormal dihydrorhodamine (DHR) assay and absent p47phox protein. Genetic identification of NCF1 mutations is complicated by adjacent highly conserved (>98%) pseudogenes (NCF1B and NCF1C). NCF1 has GTGT at the start of exon 2, whereas the pseudogenes each delete 1 GT (ΔGT). In p47phox CGD, the most common mutation is ΔGT in NCF1 (c.75_76delGT; p.Tyr26fsX26). Sequence homology between NCF1 and its pseudogenes precludes reliable use of standard Sanger sequencing for NCF1 mutations and for confirming carrier status. We first established by flow cytometry that neutrophils from p47phox CGD patients had negligible p47phox expression, whereas those from p47phox CGD carriers had ∼60% of normal p47phox expression, independent of the specific mutation in NCF1 We developed a droplet digital polymerase chain reaction (ddPCR) with 2 distinct probes, recognizing either the wild-type GTGT sequence or the ΔGT sequence. A second ddPCR established copy number by comparison with the single-copy telomerase reverse transcriptase gene, TERT We showed that 84% of p47phox CGD patients were homozygous for ΔGT NCF1 The ddPCR assay also enabled determination of carrier status of relatives. Furthermore, only 79.2% of normal volunteers had 2 copies of GTGT per 6 total (NCF1/NCF1B/NCF1C) copies, designated 2/6; 14.7% had 3/6, and 1.6% had 4/6 GTGT copies. In summary, flow cytometry for p47phox expression quickly identifies patients and carriers of p47phox CGD, and genomic ddPCR identifies patients and carriers of ΔGT NCF1, the most common mutation in p47phox CGD.

Gene expression patterns of chicken neuregulin 3 in association with copy number variation and frameshift deletion

Background

Neuregulin 3 (NRG3) plays a key role in central nervous system development and is a strong candidate for human mental disorders. Thus, genetic variation in NRG3 may have some impact on a variety of phenotypes in non-mammalian vertebrates. Recently, genome-wide screening for short insertions and deletions in chicken (Gallus gallus) genomes has provided useful information about structural variation in functionally important genes. NRG3 is one such gene that has a putative frameshift deletion in exon 2, resulting in premature termination of translation. Our aims were to characterize the structure of chicken NRG3 and to compare expression patterns between NRG3 isoforms.

Results

Depending on the presence or absence of the 2-bp deletion in chicken NRG3, 3 breeds (red junglefowl [RJF], Boris Brown [BB], and Hinai-jidori [HJ]) were genotyped using flanking primers. In the commercial breeds (BB and HJ), approximately 45% of individuals had at least one exon 2 allele with the 2-bp deletion, whereas there was no deletion allele in RJF. The lack of a homozygous mutant indicated the existence of duplicated NRG3 segments in the chicken genome. Indeed, highly conserved elements consisting of exon 1, intron 1, exon 2, and part of intron 2 were found in the reference RJF genome, and quantitative PCR detected copy number variation (CNV) between breeds as well as between individuals. The copy number of conserved elements was significantly higher in chicks harboring the 2-bp deletion in exon 2. We identified 7 novel transcript variants using total mRNA isolated from the amygdala. Novel isoforms were found to lack the exon 2 cassette, which probably harbored the premature termination codon. The relative transcription levels of the newly identified isoforms were almost the same between chick groups with and without the 2-bp deletion, while chicks with the deletion showed significant suppression of the expression of previously reported isoforms.

Conclusions

A putative frameshift deletion and CNV in chicken NRG3 are structural mutations that occurred before the establishment of commercial chicken lines. Our results further suggest that the putative frameshift deletion in exon 2 may potentially affect the expression level of particular isoforms of chicken NRG3.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0537-z) contains supplementary material, which is available to authorized users.

A single nucleotide polymorphism in theNCF1gene leading to reduced oxidative burst is associated with systemic lupus erythematosus

Objectives

Ncf1polymorphisms leading to low production of reactive oxygen species (ROS) are strongly associated with autoimmune diseases in animal models. The humanNCF1gene is very complex with both functional and non-functional gene copies and genotyping requires assays specific for functionalNCF1genes. We aimed at investigating association and function of the missense single nucleotide polymorphism (SNP), rs201802880 (here denoted NCF1-339) inNCF1with systemic lupus erythematosus (SLE).

Methods

We genotyped the NCF1-339 SNP in 973 Swedish patients with SLE and 1301 controls, using nested PCR and pyrosequencing. ROS production and gene expression of type 1 interferon-regulated genes were measured in isolated cells from subjects with different NCF1-339 genotypes.

Results

We found an increased frequency of the NCF1-339 T allele in patients with SLE, 11% compared with 4% in controls, OR 3.0, 95% CI 2.4 to 3.9, p=7.0×10−20. The NCF1-339 T allele reduced extracellular ROS production in neutrophils (p=0.004) and led to an increase expression of type 1 interferon-regulated genes. In addition, the NCF1-339 T allele was associated with a younger age at diagnosis of SLE; mean age 30.3 compared with 35.9, p=2.0×1−6.

Conclusions

These results clearly demonstrate that a genetically controlled reduced production of ROS increases the risk of developing SLE and confirm the hypothesis that ROS regulate chronic autoimmune inflammatory diseases.

Gene-edited pseudogene resurrection corrects p47phox-deficient chronic granulomatous disease

Pseudogenes are duplicated genes with mutations rendering them nonfunctional. For single-gene disorders with homologous pseudogenes, the pseudogene might be a target for genetic correction. Autosomal-recessive p47phox-deficient chronic granulomatous disease (p47-CGD) is a life-threatening immune deficiency caused by mutations in NCF1, a gene with 2 pseudogenes, NCF1B and NCF1C. The most common NCF1 mutation, a GT deletion (ΔGT) at the start of exon 2 (>90% of alleles), is constitutive to NCF1B and NCF1C. NCF1 ΔGT results in premature termination, undetectable protein expression, and defective production of antimicrobial superoxide in neutrophils. We examined strategies for p47-CGD gene correction using engineered zinc-finger nucleases targeting the exon 2 ΔGT in induced pluripotent stem cells or CD34+ hematopoietic stem cells derived from p47-CGD patients. Correction of ΔGT in NCF1 pseudogenes restores oxidase function in p47-CGD, providing the first demonstration that targeted restoration of pseudogene function can correct a monogenic disorder.

Genetic disorders coupled to ROS deficiency

Maintaining the redox balance between generation and elimination of reactive oxygen species (ROS) is critical for health. Disturbances such as continuously elevated ROS levels will result in oxidative stress and development of disease, but likewise, insufficient ROS production will be detrimental to health. Reduced or even complete loss of ROS generation originates mainly from inactivating variants in genes encoding for NADPH oxidase complexes. In particular, deficiency in phagocyte Nox2 oxidase function due to genetic variants (CYBB, CYBA, NCF1, NCF2, NCF4) has been recognized as a direct cause of chronic granulomatous disease (CGD), an inherited immune disorder. More recently, additional diseases have been linked to functionally altered variants in genes encoding for other NADPH oxidases, such as for DUOX2/DUOXA2 in congenital hypothyroidism, or for the Nox2 complex, NOX1 and DUOX2 as risk factors for inflammatory bowel disease. A comprehensive overview of novel developments in terms of Nox/Duox-deficiency disorders is presented, combined with insights gained from structure-function studies that will aid in predicting functional defects of clinical variants.

Functional pseudogenes inhibit the superoxide production

We recently discovered a dynamic copy number variation on the NCF1 (neutrophil cytosolic factor 1) pseudogenes in human populations. In this study, we investigated whether these pseudogenes are functional or junk as described before. We sequenced the RNAs transcribed from the genome of this locus, and discovered over 10 splicing isoforms from the NCF1 pseudogenes. We cloned 4 splicing isoforms into expression vectors and introduced them into human vascular endothelial cells by transient transfection. We then used two chemical approaches to measure the superoxide production in the cells with and without these pseudogene overexpression. Our data showed that three pseudogene splicing products remarkably reduced the superoxide production after the GFP (Green fluorescent protein) normalization. We used an anti-HA (Hemagglutinin A) tag antibody to stain the cells and confirmed that the proteins transcribed from the NCF1 pseudogene were exclusively localized in the cytoplasm in the perinuclear area in the transient transfection assays. We further examined the tissue distribution of these splicing isoforms of NCF1 pseudogenes in human tissues by quantitative real-time PCR analysis. Our data showed that although these splicing variants are ubiquitously expressed in non-immune tissues in human, they seem to be under a tight control of transcription regulation and show a non-random distribution pattern across tissues. This study challenges the concept that pseudogenes in human genome are only junks without biological functions. Moreover, it suggests that those pseudogenes in human genome may serve as a natural resource for novel drug discovery.

Genetic mutations potentially cause two novel NCF1 splice variants up-regulated in the mammary gland, blood and neutrophil of cows infected by Escherichia coli

Neutrophil cytosolic factor 1 (NCF1) plays a crucial role in host defense against microbial pathogens. In this study, we investigated the potential alternative splicing patterns, expression and splice-relevant single nucleotide polymorphisms (SNPs) of the bovine NCF1 gene to increase insights into its potential role against bovine mastitis caused by Escherichia coli infection. Using RT-PCR and clone sequencing methods, we found two novel splice variants designed as NCF1-TV1 (retained intron 6) and NCF1-TV2 (retained part of intron 8), respectively, encoding two putative truncated proteins (239AA and 283AA). Two splice variants were drastically up-regulated in the mastitis-infected cows' mammary tissues, blood and neutrophils compared with these of healthy cows using real-time RT-PCR. Genomic sequencing analysis identified four novel SNPs g.10112 G>A, g.10766 T>C, SNPs g.12085 G>A and g.12430 T>C at the ends of intron 6 and intron 8 of NCF1. ESE motif predicted that the SNP (g.10766 T>C) might affect the binding with splicing-related factors and subsequently caused the production of aberrant splice variant NCF1-TV1, which is one of the potential reasons that the functional SNP was associated with increased milk somatic cell score in cow. Our results would help in better understanding the NCF1 gene function in the process against pathogen infection, and the effect of splicing-related SNP on the production of aberrant splice variant and careful functional characterization in dairy cattle.

Different unequal cross-over events between NCF1 and its pseudogenes in autosomal p47(phox)-deficient chronic granulomatous disease

Chronic granulomatous disease (CGD) is a rare congenital disorder in which phagocytes cannot generate superoxide (O2(-)) and other microbicidal oxidants due to mutations in one of the five components of the O2(-)-generating NADPH oxidase complex. The most common autosomal subtype of CGD is caused by mutations in NCF1, encoding the NADPH subunit p47(phox). Usually, these mutations are the result of unequal exchange of chromatid between NCF1 and one of its two pseudogenes. We have now investigated in detail the breakpoints within or between these (pseudo) NCF1 genes in 43 families with p47(phox)-deficient CGD by means of multiplex ligase-dependent probe amplification (MLPA). In 24 families the patients totally lacked NCF1 sequences, indicating that in these families the cross-over points are located between NCF1 and its pseudogenes. Six other families were compound heterozygous for a total NCF1 deletion and another mutation in NCF1 on the other allele. In 8 families, the patients lacked NCF1 exons 1-4 but had retained NCF1 exons 6-10, indicating that a cross-over point is located within NCF1 between exons 4 and 6. Similarly, in 4 families a cross-over point was located within NCF1 between exons 2 and 4. Similar cross-overs, in heterozygous form, were observed in family members of the patients. Several patients were compound heterozygous for total and partial NCF1 deletions. Thus, at least three different cross-over points exist within the NCF1 gene cluster, indicating that autosomal p47(phox)-deficient CGD is genetically heterogeneous but can be dissected in detail by MLPA.

Diabetes, renal and cardiovascular disease in p47 phox-/- chronic granulomatous disease

Chronic granulomatous disease is a rare immunodeficiency due to defects in the phagocyte NADPH oxidase. The X-linked form (gp91 (phox) deficiency) accounts for about 70 % of cases; autosomal recessive p47 (phox) deficiency accounts for about 25 % of cases. We identified a 10 % incidence of diabetes in p47 (phox) deficient CGD, but none in X-linked CGD. Renal and cardiovascular diseases were also higher in p47 (phox) deficiency. p47 (phox) deficient CGD has non-infectious morbidities distinct from those in X-linked CGD.

Copy number variations of the ATP-binding cassette transporter ABCC6 gene and its pseudogenes

BACKGROUND

The ATP-binding cassette transporter ABCC6 gene is located on chromosome 16 between its two pseudogenes (ABCC6P1 and ABCC6P2). Previously, we have shown that ABCC6P1 is transcribed and affects ABCC6 at the transcriptional level. In this study we aimed to determine copy number variations of ABCC6, ABCC6P1 and ABCC6P2 in different populations. Moreover, we sought to study the transcription pattern of ABCC6 and ABCC6 pseudogenes in 39 different human tissues.

FINDINGS

Genomic DNA from healthy individuals from five populations, Chinese (n = 24), Middle East (n = 20), Mexicans (n = 24), Caucasians (n = 50) and Africans (n = 24), were examined for copy number variations of ABCC6 and its pseudogenes by pyrosequencing and quantitative PCR. Copy number variation of ABCC6 was very rare (2/142; 1.4%). However, one or three copies of ABCC6P1 were relatively common (3% and 8%, respectively). Only one person had a single copy of ABCC6P2 while none had three copies. In Chinese, deletions or duplications of ABCC6P1 were more frequent than in any other population (9/24; 37.5%). The transcription pattern of ABCC6P2 was highly similar to ABCC6 and ABCC6P1, with highest transcription in liver and kidney. Interestingly, the total transcription level of pseudogenes, ABCC6P1 + ABCC6P2, was higher than ABCC6 in most tissues, including liver and kidney.

CONCLUSIONS

Copy number variations of the ABCC6 pseudogenes are quite common, especially in populations of Chinese ancestry. The expression pattern of ABCC6P2 in 39 human tissues was highly similar to that of ABCC6 and ABCC6P1 suggesting similar regulatory mechanisms for ABCC6 and its pseudogenes.

Copy number variation of the gene NCF1 is associated with rheumatoid arthritis

AIMS

The aim of this study was to investigate genetic variants in the gene neutrophil cytosolic factor 1 (NCF1) for association with rheumatoid arthritis (RA). In rodent models, a single-nucleotide polymorphism (SNP) in Ncf1 has been shown to be a major locus regulating severity of arthritis. Ncf1 encodes one of five subunits of the NADPH oxidase complex. In humans the genomic structure of NCF1 is complex, excluding it from genome-wide association screens and complicating genetic analysis. In addition to copy number variation of NCF1, there are also two nonfunctional pseudogenes, nearly identical in sequence to NCF1. We have characterized copy number variation and SNPs in NCF1, and investigated these variants for association with RA.

RESULTS

We find that RA patients are less likely to have an increased copy number of NCF1, 7.6%, compared with 11.6% in controls; p=0.037. We also show that the T-allele of NCF1-339 (rs13447) is expressed in NCF1 and significantly reduces reactive oxygen species production.

INNOVATION

This is the first finding of genetic association of NCF1 with RA. The detailed characterization of genetic variants in NCF1 also helps elucidate the complexity of the NCF1 gene.

CONCLUSION

These data suggest that an increased copy number of NCF1 can be protective against developing RA and add support to previous findings of a role of NCF1 and the phagocyte NADPH oxidase complex in RA pathogenesis.

Positioning of a polymorphic quantitative trait nucleotide in the Ncf1 gene controlling oxidative burst response and arthritis severity in rats

The Ncf1 gene, encoding the P47(PHOX) protein that regulates production of reactive oxygen species (ROS) by the phagocyte NADPH oxidase (NOX2) complex, is associated with autoimmunity and arthritis severity in rats. We have now identified that the single-nucleotide polymorphism (SNP) resulting in an M153T amino acid substitution mediates arthritis resistance and thus explains the molecular polymorphism underlying the earlier identified Ncf1 gene effect. We identified the SNP in position 153 to regulate ROS production using COS(PHOX) cells transfected with mutated Ncf1. To determine the role of this SNP for control of arthritis, we used the Wistar strain, identified to carry only the postulated arthritis resistant SNP in position 153. When this Ncf1 allele was backcrossed to the arthritis susceptible DA strain, both granulocyte ROS production and arthritis resistance were restored. Position 153 is located in the hinge region between the PX and SH3 domains of P47(PHOX). Mutational analysis of this position revealed a need for an -OH group in the side chain but we found no evidence for phosphorylation. The polymorphism did not affect assembly of the P47(PHOX)/P67(PHOX) complex in the cytosol or membrane localization, but is likely to operate downstream of assembly, affecting activity of the membrane NOX2 complex.