Genome-wide scan for loci influencing quantitative immune response traits in the Belém family study: comparison of methods and summary of results.

Here we report the results from a genome-wide linkage scan to identify genes and chromosomal regions that influence quantitative immune response traits, using multi-case leprosy and tuberculosis families from north-eastern Brazil. Total plasma IgE, antigen-specific IgG to Mycobacterium leprae soluble antigen (MLSA), M. tuberculosis soluble antigen (MTSA) and M. tuberculosis purified protein derivative (PPD), and antigen-specific lymphocyte proliferation (stimulation index or SI) and interferon-gamma (IFN-gamma) release to MLSA and PPD, were measured in 16 tuberculosis (184 individuals) and 21 leprosy (177 individuals) families. The individuals were genotyped at 382 autosomal microsatellite markers across the genome. The adjusted immune-response phenotypes were analysed using a variety of variance components and regression-based methods. These analyses highlighted a number of practical issues and problems with regard to implementation of the methods and, interestingly, differences were observed between several standard statistical and genetic analysis packages used. From this we determined that, for this set of traits in these pedigrees, significant p values for linkage using variance components analysis, supported by significance using the Visscher-Hopper modification of the Haseman-Elston method, provided the most compelling evidence for linkage. Using these criteria, linkage (5.8 x 10(-5) < p < 0.008) was seen for: total plasma IgE on chromosome 2; IgG to MLSA on chromosomes 8, 17 and 21; IgG to PPD on chromosome 12; SI to PPD on chromosome 1; IFN-gamma to MLSA on chromosomes 6, 7, 10, 12 and 14; and IFN-gamma to PPD on chromosomes 1, 16 and 19.

Evidence for a cluster of genes on chromosome 17q11-q21 controlling susceptibility to tuberculosis and leprosy in Brazilians.

The region of conserved synteny on mouse chromosome 11/human 17q11-q21 is known to carry a susceptibility gene(s) for intramacrophage pathogens. The region is rich in candidates including NOS2A, CCL2/MCP-1, CCL3/MIP-1alpha, CCL4/MIP-1beta, CCL5/RANTES, CCR7, STAT3 and STAT5A/5B. To examine the region in man, we studied 92 multicase tuberculosis (627 individuals) and 72 multicase leprosy (372 individuals) families from Brazil. Multipoint nonparametric analysis (ALLEGRO) using 16 microsatellites shows two peaks of linkage for leprosy at D17S250 (Z(lr) score 2.34; P=0.01) and D17S1795 (Z(lr) 2.67; P=0.004) and a single peak for tuberculosis at D17S250 (Z(lr) 2.04; P=0.02). Combined analysis shows significant linkage (peak Z(lr) 3.38) at D17S250, equivalent to an allele sharing LOD score 2.48 (P=0.0004). To determine whether one or multiple genes contribute, 49 informative single nucleotide polymorphisms were typed in candidate genes. Family-based allelic association testing that was robust to family clustering demonstrated significant associations with tuberculosis susceptibility at four loci separated by intervals (NOS2A-8.4 Mb-CCL18-32.3 kb-CCL4-6.04 Mb-STAT5B) up to several Mb. Stepwise conditional logistic regression analysis using a case/pseudo-control data set showed that the four genes contributed separate main effects, consistent with a cluster of susceptibility genes across 17q11.2.

Genome-wide scans for leprosy and tuberculosis susceptibility genes in Brazilians.

Genome-wide scans were conducted for tuberculosis and leprosy per se in Brazil. At stage 1, 405 markers (10 cM map) were typed in 16 (178 individuals) tuberculosis and 21 (173 individuals) leprosy families. Nonparametric multipoint analysis detected 8 and 9 chromosomal regions respectively with provisional evidence (P<0.05) for linkage. At stage 2, 58 markers from positive regions were typed in a second set of 22 (176 individuals) tuberculosis families, with 22 additional markers typed in all families; 42 positive markers in 50 (192 individuals) new leprosy families, and 30 additional markers in all families. Three regions (10q26.13, 11q12.3, 20p12.1) retained suggestive evidence (peak LOD scores 1.31, 1.85, 1.78; P=0.007, 0.0018, 0.0021) for linkage to tuberculosis, 3 regions (6p21.32, 17q22, 20p13) to leprosy (HLA-DQA, 3.23, P=5.8 x 10(-5); D17S1868, 2.38, P=0.0005; D20S889, 1.51, P=0.004). The peak at D20S889 for leprosy is 3.5 Mb distal to that reported at D20S115 for leprosy in India. (151 words).

Evidence for a cluster of genes on chromosome 17q11-q21 controlling susceptibility to tuberculosis and leprosy in Brazilians

The region of conserved synteny on mouse chromosome 11/human 17q11-q21 is known to carry a susceptibility gene(s) for intramacrophage pathogens. The region is rich in candidates including NOS2A, CCL2/MCP-1, CCL3/MIP-1 alpha, CCL4/MIP-1 beta, CCL5/RANTES, CCR7, STAT5A/5B. To examine the region in man, we studied 92 multicase tuberculosis (627 individuals) and 72 multicase leprosy (372 indiciduals) families from Brazil. Multipoint nonparametric analysis (ALLEGRO) using 16 microsatellites shows two peaks of linkage for leprosy at D17S250 (Zir score 2.34; P=0.01) and D17S1795 (Zir 2.67; P=O.004) and a single peack for tuberculosis at D17S250 (Zir 2.04; P=0.02). Combined analysis shows significant linkage (peak Zir 3.38) at D17S250, equivalent to an allele sharing LOD score 2.48 (P=0.0004). To determine whether one or multiple genes contribute, 49 informative single nucleotide polymorphisms were typed in candidate genes. Family-based allelic association testing that was robust to family clustering demonstrated significant associations with tuberculosis susceptibility at four loci separated by intervals (NOS2A-8.4 Mb-CCL 18-32.3 kb-CCL4-6.04 Mb-STAT5B) up to several Mb. Stepwise conditional logistic regression analysis using a case/pseudo-control data set showed that the four genes contributed separate main effects, consistent with a cluster of susceptibilitty genes acros 17q11.2

Genome-wide scans for leprosy and tuberculosis susceptibility genes in Brazilians

Genome-wide scans were conducted for tuberculosis and leprosy per se in Brazil. At stage 1,405 markers (10 cM map) were typed in 16 (178 individuals) tuberculosis and 21 (173 individuals) leprosy families. Nonparametric multipoint analysis detected 8 and 9 chromosomal regions respectively with provisional evidence (P<0.05) for linkage. A stage 2, 58 markers from positive regions were typed in a second set of 22 (176 individuals) tuberculosis families, with 22 additional markers types in all families; 42 positive markers in 50 (192 individuals) new leprosy families, and 30 additional markers in all families. Three regions (10q26.13, 11q12.3, 20p12.1) retained suggestive evidence (peak LOD scores 1.31, 1.78, 1.78; P=0.007, 0.0018, 0.0021) for linkage to tuberculosis, 3 regions (6p21.32, 17q22, 20p13) to leprosy (HLA-DQA, 3.23, P=5.8 x 10-5; D17S1868.2.38, P=0.0005; D20S889, 1.51, P=0.004). The peak at D20S889 for leprosy is 3.5 Mb distal to that reported at D20S115 for leprosy in India

Association and linkage of leprosy phenotypes with HLA class II and tumour necrosis factor genes.

Previous analyses indicate major gene control of susceptibility to leprosy per se and the HLA class II region has been implicated in determining susceptibility and control of clinical phenotype. Segregation analysis using data from 76 Brazilian leprosy multi-case pedigrees (1166 individuals) supported a two locus model as the best fit: a recessive major gene and a recessive modifier gene(s) (single locus vs two locus model, P = 0.0007). Combined segregation and linkage analysis to the major locus, showed strong linkage to HLA class II (HLA-DQB1 P = 0.000002, HLA-DQA1 P = 0.000002, HLA-DRB1 P = 0.0000003) and tumour necrosis factor genes (TNF P = 0.00002, LTA P = 0.003). Extended transmission disequilibrium testing, using multiple affected family members, demonstrated that the common allele TNF*1 of the -308 promoter region polymorphism showed linkage and/or association with disease per se, at a high level of significance (P < 0.0001). Two locus transmission disequilibrium testing suggested susceptibility (TNF*1/LTA*2) and protective (TNF*2/LTA*2) haplotypes in the class iii region. Taken together the segregation and HLA analyses suggest the possibility of more than one susceptibility locus in the MHC.

Immunogenetics of leishmanial and mycobacterial infections: the Belem Family Study.

In the 1970s and 1980s, analysis of recombinant inbred, congenic and recombinant haplotype mouse strains permitted us to effectively 'scan' the murine genome for genes controlling resistance and susceptibility to leishmanial infections. Five major regions of the genome were implicated in the control of infections caused by different Leishmania species which, because they show conserved synteny with regions of the human genome, immediately provides candidate gene regions for human disease susceptibility genes. A common intramacrophage niche for leishmanial and mycobacterial pathogens, and a similar spectrum of immune response and disease phenotypes, also led to the prediction that the same genes/candidate gene regions might be responsible for genetic susceptibility to mycobacterial infections such as leprosy and tuberculosis. Indeed, one of the murine genes (Nramp1) was identified for its role in controlling a range of intramacrophage pathogens including leishmania, salmonella and mycobacterium infections. In recent studies, multicase family data on visceral leishmaniasis and the mycobacterial diseases, tuberculosis and leprosy, have been collected from north-eastern Brazil and analysed to determine the role of these candidate genes/regions in determining disease susceptibility. Complex segregation analysis provides evidence for one or two major genes controlling susceptibility to tuberculosis in this population. Family-based linkage analyses (combined segregation and linkage analysis; sib-pair analysis), which have the power to detect linkage between marker loci in candidate gene regions and the putative disease susceptibility genes over 10-20 centimorgans, and transmission disequilibrium testing, which detects allelic associations over 1 centimorgan (ca. 1 megabase), have been used to examine the role of four regions in determining disease susceptibility and/or immune response phenotype. Our results demonstrate: (i) the major histocompatibility complex (MHC: H-2 in mouse, HLA in man: mouse chromosome 17/human 6p; candidates class II and class III including TNF alpha/beta genes) shows both linkage to, and allelic association with, leprosy per se, but is only weakly associated with visceral leishmaniasis and shows neither linkage to nor allelic association with tuberculosis; (ii) no evidence for linkage between NRAMP1, the positionally cloned candidate for the murine macrophage resistance gene Ity/Lsh/Bcg (mouse chromosome 1/human 2q35), and susceptibility to tuberculosis or visceral leishmaniasis could be demonstrated in this Brazilian population; (iii) the region of human chromosome 17q (candidates NOS2A, SCYA2-5) homologous with distal mouse chromosome 11, originally identified as carrying the Scl1 gene controlling healing versus nonhealing responses to Leishmania major, is linked to tuberculosis susceptibility; and (iv) the 'T helper 2' cytokine gene cluster (proximal murine chromosome 11/human 5q; candidates IL4, IL5, IL9, IRF1, CD14) controlling later phases of murine L. major infection, is not linked to human disease susceptibility for any of the three infections, but shows linkage to and highly significant allelic association with ability to mount an immune response to mycobacterial antigens. These studies demonstrate that the 'mouse-to-man' strategy, refined by our knowledge of the human immune response to infection, can lead to the identification of important candidate gene regions in man.

Genetic and physical mapping of 2q35 in the region of the NRAMP and IL8R genes: identification of a polymorphic repeat in exon 2 of NRAMP.

Recent interest has focused on the region of conserved synteny between mouse chromosome 1 and human 2q33-q37, particularly over the region encoding the murine macrophage resistance gene Ity/Lsh/Bcg (candidate Nramp) and members of the Il8r interleukin-8 (IL8) receptor gene cluster. In this paper, identification of a restriction fragment length polymorphism in the IL8RB gene in 35 pedigrees previously typed for markers in the 2q33-q37 interval provided evidence (lod scores > 3) for linkage between IL8RB and the 2q34-q35 markers FN1, TNP1, VIL1, and DES. Physical mapping, using yeast artificial chromosomes isolated with VIL1, confirmed that IL8RA, IL8RB, and the IL8RB pseudogene map within the NRAMP-VIL1 interval, with the physical distance (155 kb) from 5' LSH to 3' VIL1 representing approximately 3-fold that observed in the mouse. Partial sequencing of NRAMP confirmed the presence of the N-terminal proline/serine-rich putative SH3 binding domain in exon 2 of the human gene. Further analysis of Brazilian leprosy and visceral leishmaniasis pedigrees identified a rare second allele varying in a 9-nucleotide repeat motif of the exon 2 sequence but segregating independently of the disease phenotype.