Definition of arthritis candidate risk genes by combining rat linkage-mapping results with human case-control association data

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

Identification of candidate risk gene variations by whole-genome sequence analysis of four rat strains commonly used in inflammation research

Background

The DA rat strain is particularly susceptible to the induction of a number of chronic inflammatory diseases, such as models for rheumatoid arthritis and multiple sclerosis. Here we sequenced the genomes of two DA sub-strains and two disease resistant strains, E3 and PVG, previously used together with DA strains in genetically segregating crosses.

Results

The data uncovers genomic variations, such as single nucleotide variations (SNVs) and copy number variations that underlie phenotypic differences between the strains. Comparisons of regional differences between the two DA sub-strains identified 8 genomic regions that discriminate between the strains that together cover 38 Mbp and harbor 302 genes. We analyzed 10 fine-mapped quantitative trait loci and our data implicate strong candidates for genetic variations that mediate their effects. For example we could identify a single SNV candidate in a regulatory region of the gene Il21r, which has been associated to differential expression in both rats and human MS patients. In the APLEC complex we identified two SNVs in a highly conserved region, which could affect the regulation of all APLEC encoded genes and explain the polygenic differential expression seen in the complex. Furthermore, the non-synonymous SNV modifying aa153 of the Ncf1 protein was confirmed as the sole causative factor.

Conclusion

This complete map of genetic differences between the most commonly used rat strains in inflammation research constitutes an important reference in understanding how genetic variations contribute to the traits of importance for inflammatory diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-391) contains supplementary material, which is available to authorized users.

Genetic variability in the rat Aplec C-type lectin gene cluster regulates lymphocyte trafficking and motor neuron survival after traumatic nerve root injury

BACKGROUND

C-type lectin (CLEC) receptors are important for initiating and shaping immune responses; however, their role in inflammatory reactions in the central nervous system after traumatic injuries is not known. The antigen-presenting lectin-like receptor gene complex (Aplec) contains a few CLEC genes, which differ genetically among inbred rat strains. It was originally thought to be a region that regulates susceptibility to autoimmune arthritis, autoimmune neuroinflammation and infection.

METHODS

The inbred rat strains DA and PVG differ substantially in degree of spinal cord motor neuron death following ventral root avulsion (VRA), which is a reproducible model of localized nerve root injury. A large F2 (DAxPVG) intercross was bred and genotyped after which global expressional profiling was performed on spinal cords from F2 rats subjected to VRA. A congenic strain, Aplec, created by transferring a small PVG segment containing only seven genes, all C-type lectins, ontoDA background, was used for further experiments together with the parental strains.

RESULTS

Global expressional profiling of F2 (DAxPVG) spinal cords after VRA and genome-wide eQTL mapping identified a strong cis-regulated difference in the expression of Clec4a3 (Dcir3), a C-type lectin gene that is a part of the Aplec cluster. Second, we demonstrate significantly improved motor neuron survival and also increased T-cell infiltration into the spinal cord of congenic rats carrying Aplec from PVG on DA background compared to the parental DA strain. In vitro studies demonstrate that the Aplec genes are expressed on microglia and upregulated upon inflammatory stimuli. However, there were no differences in expression of general microglial activation markers between Aplec and parental DA rats, suggesting that the Aplec genes are involved in the signaling events rather than the primary activation of microglia occurring upon nerve root injury.

CONCLUSIONS

In summary, we demonstrate that a genetic variation in Aplec occurring among inbred strains regulates both survival of axotomized motor neurons and the degree of lymphocyte infiltration. These results demonstrate a hitherto unknown role for CLECs for intercellular communication that occurs after damage to the nervous system, which is relevant for neuronal survival.

Superoxide production by macrophages and T cells is critical for the induction of autoreactivity and type 1 diabetes

OBJECTIVE

The role of reactive oxygen species (ROS) and their dissipation in type 1 diabetes pathogenesis have garnered considerable controversy. Our recent work has demonstrated the importance of NADPH oxidase (NOX) activity for type 1 diabetes development and modulating T-cell autoreactivity. We previously linked decreased monocyte ROS with diabetes resistance in the alloxan-resistant mouse, and NOD-Ncf1(m1J) mice with a genetic ablation of NOX activity had reduced and delayed type 1 diabetes compared with NOD mice.

RESEARCH DESIGN AND METHODS

To determine the required cellular sources of ROS that are necessary for type 1 diabetes initiation, we used antibody depletion and adoptive transfer experiments into NOD and NOD-Scid females, respectively. After receiving treatment, female mice were monitored for hyperglycemia and overt diabetes.

RESULTS

Depletion of macrophages and neutrophils fully protected NOD mice from type 1 diabetes. However, elimination of neutrophils alone showed no significant reduction or delay. Type 1 diabetes induction in NOD-Scid mice by adoptive transfer with NOD-Ncf1(m1J) splenocytes was significantly delayed compared with NOD splenocytes, suggesting macrophage ROS and modulation of effector responses are critical for diabetes. The adaptive immune response was also altered by the absence of NOX activity, as purified T cells from NOD-Ncf1(m1J) mice exhibited delayed transfer kinetics. Cotransfer experiments demonstrated the defect was intrinsic to NOX-deficient CD8(+) T cells. After stimulation, cytotoxic T cells exhibited decreased effector function in the absence of superoxide production.

CONCLUSIONS

These data demonstrate that the impaired autoreactive response of NOX-deficient NOD-Ncf1(m1J) immune system results from an alteration in the antigen-presenting cell-T-cell axis rather than failure of neutrophils to act as effector cells and that ROS signaling is important for the initiation of β-cell-directed autoimmunity by T cells.

NADPH oxidase deficiency regulates Th lineage commitment and modulates autoimmunity

Reactive oxygen species are used by the immune system to eliminate infections; however, they may also serve as signaling intermediates to coordinate the efforts of the innate and adaptive immune systems. In this study, we show that by eliminating macrophage and T cell superoxide production through the NADPH oxidase (NOX), T cell polarization was altered. After stimulation with immobilized anti-CD3 and anti-CD28 or priming recall, T cells from NOX-deficient mice exhibited a skewed Th17 phenotype, whereas NOX-intact cells produced cytokines indicative of a Th1 response. These findings were corroborated in vivo by studying two different autoimmune diseases mediated by Th17 or Th1 pathogenic T cell responses. NOX-deficient NOD mice were Th17 prone with a concomitant susceptibility to experimental allergic encephalomyelitis and significant protection against type 1 diabetes. These data validate the role of superoxide in shaping Th responses and as a signaling intermediate to modulate Th17 and Th1 T cell responses.

Rare variation at the TNFAIP3 locus and susceptibility to rheumatoid arthritis

Genome-wide association studies (GWAS) conducted using commercial single nucleotide polymorphisms (SNP) arrays have proven to be a powerful tool for the detection of common disease susceptibility variants. However, their utility for the detection of lower frequency variants is yet to be practically investigated. Here we describe the application of a rare variant collapsing method to a large genome-wide SNP dataset, the Wellcome Trust Case Control Consortium rheumatoid arthritis (RA) GWAS. We partitioned the data into gene-centric bins and collapsed genotypes of low frequency variants (defined here as MAF ≤0.05) into a single count coupled with univariate analysis. We then prioritised gene regions for further investigation in an independent cohort of 3,355 cases and 2,427 controls based on rare variant signal p value and prior evidence to support involvement in RA. A total of 14,536 gene bins were investigated in the primary analysis and signals mapping to the TNFAIP3 and chr17q24 loci were selected for further investigation. We detected replicating association to low frequency variants in the TNFAIP3 gene (combined p = 6.6 × 10⁻⁶). Even though rare variants are not well-represented and can be difficult to genotype in GWAS, our study supports the application of low frequency variant collapsing methods to genome-wide SNP datasets as a means of exploiting data that are routinely ignored.

Genome-wide association studies and the problem of relatedness among advanced intercross lines and other highly recombinant populations

Model organisms offer many advantages for the genetic analysis of complex traits. However, identification of specific genes is often hampered by a lack of recombination between the genomes of inbred progenitors. Recently, genome-wide association studies (GWAS) in humans have offered gene-level mapping resolution that is possible because of the large number of accumulated recombinations among unrelated human subjects. To obtain analogous improvements in mapping resolution in mice, we used a 34th generation advanced intercross line (AIL) derived from two inbred strains (SM/J and LG/J). We used simulations to show that familial relationships among subjects must be accounted for when analyzing these data; we then used a mixed model that included polygenic effects to address this problem in our own analysis. Using a combination of F(2) and AIL mice derived from the same inbred progenitors, we identified genome-wide significant, subcentimorgan loci that were associated with methamphetamine sensitivity, (e.g., chromosome 18; LOD = 10.5) and non-drug-induced locomotor activity (e.g., chromosome 8; LOD = 18.9). The 2-LOD support interval for the former locus contains no known genes while the latter contains only one gene (Csmd1). This approach is broadly applicable in terms of phenotypes and model organisms and allows GWAS to be performed in multigenerational crosses between and among inbred strains where familial relatedness is often unavoidable.

Finemapping of the arthritis QTL Pia7 reveals co-localization with Oia2 and the APLEC locus

In this study, we sought to determine the effect of the quantitative trait locus Pia7 on arthritis severity. The regulatory locus derived from the arthritis-resistant E3 rat strain was introgressed into the arthritis-susceptibility DA strain through continuous backcrossing. Congenic rats were studied for their susceptibility to experimental arthritis using pristane and adjuvant oil. In addition, cell number and function of various leukocyte populations were analyzed either under naive or stimulated conditions. We found that the minimal congenic fragment of DA.E3-Pia7 rats overlapped with the minimal fragment in DA.PVG-Oia2 congenic rats, which has been positionally cloned to the antigen-presenting lectin-like receptor complex (APLEC) genes. DA.E3-Pia7 congenic rats were protected from both PIA and OIA, but the protection was more pronounced in OIA. In adoptive transfer experiments we observed that the Pia7 locus controlled the priming of arthritogenic T cells and not the effector phase. In addition, Pia7 congenic rats had a significant higher frequency of B cells and granulocytes as well as TNFalpha production after stimulation, indicating a higher activation state of cells of the innate immune system. In conclusion, this study shows that the APLEC locus is a major locus regulating the severity of experimentally induced arthritis in rats.