The Maximum-Likelihood-Binomial method revisited: a robust approach for model-free linkage analysis of quantitative traits in large sibships

An eQTL variant of ZXDC is associated with IFN-γ production following Mycobacterium tuberculosis antigen-specific stimulation

There is a large inter-individual variability in the response to Mycobacterium tuberculosis infection. In previous linkage analyses, we identified a major locus on chromosome region 8q controlling IFN-γ production after stimulation with live BCG (Bacillus Calmette-Guérin), and a second locus on chromosome region 3q affecting IFN-γ production triggered by the 6-kDa early secretory antigen target (ESAT-6), taking into account the IFN-γ production induced by BCG (IFNγ-ESAT6_BCG). High-density genotyping and imputation identified ~100,000 variants within each linkage region, which we tested for association with the corresponding IFN-γ phenotype in families from a tuberculosis household contact study in France. Significant associations were replicated in a South African familial sample. The most convincing association observed was that between the IFNγ-ESAT6_BCG phenotype and rs9828868 on chromosome 3q (p = 9.8 × 10⁻⁶ in the French sample). This variant made a significant contribution to the linkage signal (p < 0.001), and a trend towards the same association was observed in the South African sample. This variant was reported to be an eQTL of the ZXDC gene, biologically linked to monocyte IL-12 production through CCL2/MCP1. The identification of rs9828868 as a genetic driver of IFNγ production in response to mycobacterial antigens provides new insights into human anti-tuberculosis immunity.

Major Loci on Chromosomes 8q and 3q Control Interferon γ Production Triggered by Bacillus Calmette-Guerin and 6-kDa Early Secretory Antigen Target, Respectively, in Various Populations

Background. Interferon γ (IFN-γ) release assays (IGRAs) provide an in vitro measurement of antimycobacterial immunity that is widely used as a test for Mycobacterium tuberculosis infection. IGRA outcomes are highly heritable in various populations, but the nature of the involved genetic factors remains unknown.

Methods. We conducted a genome-wide linkage analysis of IGRA phenotypes in families from a tuberculosis household contact study in France and a replication study in families from South Africa to confirm the loci identified.

Results. We identified a major locus on chromosome 8q controlling IFN-γ production in response to stimulation with live bacillus Calmette-Guerin (BCG; LOD score, 3.81; P = 1.40 × 10⁻⁵). We also detected a second locus, on chromosome 3q, that controlled IFN-γ levels in response to stimulation with 6-kDa early secretory antigen target, when accounting for the IFN-γ production shared with that induced by BCG (LOD score, 3.72; P = 1.8 × 10⁻⁵). Both loci were replicated in South African families, where tuberculosis is hyperendemic. These loci differ from those previously identified as controlling the response to the tuberculin skin test (TST1 and TST2) and the production of TNF-α (TNF1).

Conclusions. The identification of 2 new linkage signals in populations of various ethnic origins living in different M. tuberculosis exposure settings provides new clues about the genetic control of human antimycobacterial immunity.

Tuberculin skin test negativity is under tight genetic control of chromosomal region 11p14-15 in settings with different tuberculosis endemicities

A substantial proportion of subjects exposed to a contagious tuberculosis case display lack of tuberculin skin test (TST) reactivity. We previously mapped a major locus (TST1) controlling lack of TST reactivity in families from an area in South Africa where tuberculosis is hyperendemic. Here, we conducted a household tuberculosis contact study in a French area where the endemicity of tuberculosis is low. A genome-wide analysis of TST negativity identified a significant linkage signal (P < 3 × 10⁻⁵) in close vicinity of TST1. Combined analysis of the 2 samples increased evidence of linkage (P = 2.4 × 10⁻⁶), further implicating genetic factors located on 11p14-15. This region overlaps the TNF1 locus controlling mycobacteria-driven tumor necrosis factor α production.

Combined linkage and association studies show that HLA class II variants control levels of antibodies against Epstein-Barr virus antigens

Over 95% of the adult population worldwide is infected with Epstein-Barr virus (EBV). EBV infection is associated with the development of several cancers, including Hodgkin lymphoma (HL). Elevated levels of anti-EBV antibodies have been associated with increased risk of HL. There is growing evidence that genetic factors control the levels of antibodies against EBV antigens. Here, we conducted linkage and association studies to search for genetic factors influencing either anti-viral capsid antigen (VCA) or anti-Epstein Barr nuclear antigen-1 (EBNA-1) IgG levels in a unique cohort of 424 individuals of European origin from 119 French families recruited through a Hodgkin lymphoma (HL) patient. No major locus controlling anti-VCA antibody levels was identified. However, we found that the HLA region influenced anti-EBNA-1 IgG titers. Refined association studies in this region identified a cluster of HLA class II variants associated with anti-EBNA-1 IgG titers (e.g. p = 5×10–5 for rs9268403). The major allele of rs9268403 conferring a predisposition to high anti-EBNA-1 antibody levels was also associated with an increased risk of HL (p = 0.02). In summary, this study shows that HLA class II variants influenced anti-EBNA-1 IgG titers in a European population. It further shows the role of the same variants in the risk of HL.

Identification of a major locus, TNF1, that controls BCG-triggered tumor necrosis factor production by leukocytes in an area hyperendemic for tuberculosis

BACKGROUND

Tumor necrosis factor (TNF) is a key immune regulator of tuberculosis resistance, as exemplified by the highly increased risk of tuberculosis disease among individuals receiving TNF-blocker therapy.

METHODS

We determined the extent of TNF production after stimulation with BCG or BCG plus interferon gamma (IFN-γ) using a whole blood assay in 392 children belonging to 135 nuclear families from an area hyperendemic for tuberculosis in South Africa. We conducted classical univariate and bivariate genome-wide linkage analysis of TNF production using the data from both stimulation protocols by means of an extension of the maximum-likelihood-binomial method for quantitative trait loci to multivariate analysis.

RESULTS

Stimulation of whole blood by either BCG or BCG plus IFN-γ resulted in a range of TNF release across subjects. Extent of TNF production following both stimulation protocols was highly correlated (r = 0.81). We failed to identify genetic linkage of TNF release when considering each stimulus separately. However, using a multivariate approach, we detected a major pleiotropic locus (P < 10(-5)) on chromosome region 11p15, termed TNF locus 1 (TNF1), that controlled TNF production after stimulation by both BCG alone and BCG plus IFN-γ.

CONCLUSIONS

The TNF1 locus was mapped in the vicinity of the TST1 locus, previously identified in the same family sample, that controls tuberculin skin test (TST) negativity per se, that is, T-cell-independent resistance to Mycobacterium tuberculosis infection. This suggested that there is a connection between TST negativity per se and TNF production.