Comparison of methods to account for relatedness in genome-wide association studies with family-based data.

Approaches based on linear mixed models (LMMs) have recently gained popularity for modelling population substructure and relatedness in genome-wide association studies. In the last few years, a bewildering variety of different LMM methods/software packages have been developed, but it is not always clear how (or indeed whether) any newly-proposed method differs from previously-proposed implementations. Here we compare the performance of several LMM approaches (and software implementations, including EMMAX, GenABEL, FaST-LMM, Mendel, GEMMA and MMM) via their application to a genome-wide association study of visceral leishmaniasis in 348 Brazilian families comprising 3626 individuals (1972 genotyped). The implementations differ in precise details of methodology implemented and through various user-chosen options such as the method and number of SNPs used to estimate the kinship (relatedness) matrix. We investigate sensitivity to these choices and the success (or otherwise) of the approaches in controlling the overall genome-wide error-rate for both real and simulated phenotypes. We compare the LMM results to those obtained using traditional family-based association tests (based on transmission of alleles within pedigrees) and to alternative approaches implemented in the software packages MQLS, ROADTRIPS and MASTOR. We find strong concordance between the results from different LMM approaches, and all are successful in controlling the genome-wide error rate (except for some approaches when applied naively to longitudinal data with many repeated measures). We also find high correlation between LMMs and alternative approaches (apart from transmission-based approaches when applied to SNPs with small or non-existent effects). We conclude that LMM approaches perform well in comparison to competing approaches. Given their strong concordance, in most applications, the choice of precise LMM implementation cannot be based on power/type I error considerations but must instead be based on considerations such as speed and ease-of-use.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis.

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P(combined) = 2.76 × 10(-17) and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species.

Cytokine responses to novel antigens in a peri-urban population in Brazil exposed to Leishmania infantum chagasi.

Visceral leishmaniasis (VL) is fatal if untreated, and there are no vaccines for this disease. High levels of CD4-derived interferon-γ (IFN-γ) in the presence of low levels of interleukin-10 (IL-10) predicts vaccine success. Tumor necrosis factor-α (TNF-α) is also important in this process. We characterized human immune responses in three groups exposed to Leishmania infantum chagasi in Brazil: 1) drug-cured VL patients (recovered VL); 2) asymptomatic persons with positive Leishmania-specific delayed-type hypersensitivity skin reactions (DTH+); and 3) DTH-negative household contacts. Magnitude of DTH correlated with crude Leishmania antigen-driven IFN-γ, TNF-α, and IL-5, but not IL-10. DTH+ persons showed equivalent levels of IFN-γ, but higher levels of IL-10, to tryparedoxin peroxidase and Leishmania homolog of receptor for activated C kinase compared with recovered VL patients. The IFN-γ:IL-10 and TNF-α:IL-10 ratios were higher in recovered VL patients than in DTH+ persons. Seven of 11 novel candidates (R71, L37, N52, L302.06, M18, J41, and M22) elicited cytokine responses (36-71% of responders) in recovered VL patients and DTH+ persons. This result confirmed their putative status as cross-species vaccine/immunotherapeutic candidates.

Genetic and functional evidence implicating DLL1 as the gene that influences susceptibility to visceral leishmaniasis at chromosome 6q27.

BACKGROUND: Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. METHODS: Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. RESULTS: Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 × 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 × 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 × 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. CONCLUSIONS: DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

CXCR1 and SLC11A1 polymorphisms affect susceptibility to cutaneous leishmaniasis in Brazil: a case-control and family-based study.

BACKGROUND: L. braziliensis causes cutaneous (CL) and mucosal (ML) leishmaniasis. Wound healing neutrophil (PMN) and macrophage responses made following the bite of the vector sand fly contribute to disease progression in mice. To look at the interplay between PMN and macrophages in disease progression in humans we asked whether polymorphisms at genes that regulate their infiltration or function are associated with different clinical phenotypes. Specifically, CXCR1 (IL8RA) and CXCR2 (IL8RB) are receptors for chemokines that attract PMN to inflammatory sites. They lie 30-260 kb upstream of SLC11A1, a gene known primarily for its role in regulating macrophage activation, resistance to leishmaniasis, and wound healing responses in mice, but also known to be expressed in PMN, macrophages and dendritic cells. METHODS: Polymorphic variants at CXCR1, CXCR2 and SLC11A1 were analysed using Taqman or ABI fragment separation technologies in cases (60 CL; 60 ML), unrelated controls (n = 120), and multicase families (104 nuclear families; 88 ML, 250 CL cases) from Brazil. Logistic regression analysis, family-based association testing (FBAT) and haplotype analysis (TRANSMIT) were performed. RESULTS: Case-control analysis showed association between the common C allele (OR 2.38; 95% CI 1.23-4.57; P = 0.009) of CXCR1_rs2854386 and CL, supported by family-based (FBAT; Z score 2.002; P = 0.045) analysis (104 nuclear families; 88 ML, 250 CL cases). ML associated with the rarer G allele (Z score 1.999; P = 0.046). CL associated with a 3' insertion/deletion polymorphism at SLC11A1 (Z score 2.549; P = 0.011). CONCLUSIONS: The study supports roles for CXCR1 and SLC11A1 in the outcome of L. braziliensis infection in humans. Slc11a1 does not influence cutaneous lesion development following needle injection of Leishmania in mice, suggesting that its role here might relate to the action of PMN, macrophage and/or dendritic cells in the wound healing response to the sand fly bite. Together with the CXCR1 association, the data are consistent with hypotheses relating to the possible role of PMN in initiation of a lesion following the delivery of parasites via the sand fly bite. Association of ML with the rare derived G allele suggests that PMN also have an important positive role to play in preventing this form of the disease.

Classification and regression tree and spatial analyses reveal geographic heterogeneity in genome wide linkage study of Indian visceral leishmaniasis.

BACKGROUND: Genome wide linkage studies (GWLS) have provided evidence for loci controlling visceral leishmaniasis on Chromosomes 1p22, 6q27, 22q12 in Sudan and 6q27, 9p21, 17q11-q21 in Brazil. Genome wide studies from the major focus of disease in India have not previously been reported. METHODS AND FINDINGS: We undertook a GWLS in India in which a primary ∼10 cM (515 microsatellites) scan was carried out in 58 multicase pedigrees (74 nuclear families; 176 affected, 353 total individuals) and replication sought in 79 pedigrees (102 nuclear families; 218 affected, 473 total individuals). The primary scan provided evidence (≥2 adjacent markers allele-sharing LOD≥0.59; nominal P≤0.05) for linkage on Chromosomes 2, 5, 6, 7, 8, 10, 11, 20 and X, with peaks at 6p25.3-p24.3 and 8p23.1-p21.3 contributed to largely by 31 Hindu families and at Xq21.1-q26.1 by 27 Muslim families. Refined mapping confirmed linkage across all primary scan families at 2q12.2-q14.1 and 11q13.2-q23.3, but only 11q13.2-q23.3 replicated (combined LOD = 1.59; P = 0.0034). Linkage at 6p25.3-p24.3 and 8p23.1-p21.3, and at Xq21.1-q26.1, was confirmed by refined mapping for primary Hindu and Muslim families, respectively, but only Xq21.1-q26.1 replicated across all Muslim families (combined LOD 1.49; P = 0.0045). STRUCTURE and SMARTPCA did not identify population genetic substructure related to religious group. Classification and regression tree, and spatial interpolation, analyses confirm geographical heterogeneity for linkages at 6p25.3-p24.3, 8p23.1-p21.3 and Xq21.1-q26.1, with specific clusters of families contributing LOD scores of 2.13 (P = 0.0009), 1.75 (P = 0.002) and 1.84 (P = 0.001), respectively. CONCLUSIONS: GWLS has identified novel loci that show geographical heterogeneity in their influence on susceptibility to VL in India.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9)

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95 percent confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9).

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95% confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

Genetic and epigenetic factors at COL2A1 and ABCA4 influence clinical outcome in congenital toxoplasmosis.

BACKGROUND: Primary Toxoplasma gondii infection during pregnancy can be transmitted to the fetus. At birth, infected infants may have intracranial calcification, hydrocephalus, and retinochoroiditis, and new ocular lesions can occur at any age after birth. Not all children who acquire infection in utero develop these clinical signs of disease. Whilst severity of disease is influenced by trimester in which infection is acquired by the mother, other factors including genetic predisposition may contribute. METHODS AND FINDINGS: In 457 mother-child pairs from Europe, and 149 child/parent trios from North America, we show that ocular and brain disease in congenital toxoplasmosis associate with polymorphisms in ABCA4 encoding ATP-binding cassette transporter, subfamily A, member 4. Polymorphisms at COL2A1 encoding type II collagen associate only with ocular disease. Both loci showed unusual inheritance patterns for the disease allele when comparing outcomes in heterozygous affected children with outcomes in affected children of heterozygous mothers. Modeling suggested either an effect of mother's genotype, or parent-of-origin effects. Experimental studies showed that both ABCA4 and COL2A1 show isoform-specific epigenetic modifications consistent with imprinting. CONCLUSIONS: These associations between clinical outcomes of congenital toxoplasmosis and polymorphisms at ABCA4 and COL2A1 provide novel insight into the molecular pathways that can be affected by congenital infection with this parasite.

Y chromosome lineage- and village-specific genes on chromosomes 1p22 and 6q27 control visceral leishmaniasis in Sudan.

Familial clustering and ethnic differences suggest that visceral leishmaniasis caused by Leishmania donovani is under genetic control. A recent genome scan provided evidence for a major susceptibility gene on Chromosome 22q12 in the Aringa ethnic group in Sudan. We now report a genome-wide scan using 69 families with 173 affected relatives from two villages occupied by the related Masalit ethnic group. A primary ten-centimorgan scan followed by refined mapping provided evidence for major loci at 1p22 (LOD score 5.65; nominal p = 1.72 x 10(-7); empirical p < 1 x 10(-5); lambdaS = 5.1) and 6q27 (LOD score 3.74; nominal p = 1.68 x 10(-5); empirical p < 1 x 10(-4); lambdaS = 2.3) that were Y chromosome-lineage and village-specific. Neither village supported a visceral leishmaniasis susceptibility gene on 22q12. The results suggest strong lineage-specific genes due to founder effect and consanguinity in these recently immigrant populations. These chance events in ethnically uniform African populations provide a powerful resource in the search for genes and mechanisms that regulate this complex disease.

Polymorphism in IL1B: IL1B-511 association with tuberculosis and decreased lipopolysaccharide-induced IL-1beta in IFN-gamma primed ex-vivo whole blood assay.

To determine whether variation in two interleukin 1 family genes (IL1B and interleukin 1 receptor antagonist, IL1RN) is associated with pulmonary tuberculosis (TB), two published polymorphisms at nucleotide positions -511 and +3953 in IL1B and one in the IL1RN 86 bp VNTR were genotyped in 335 smear positive Gambian TB patients, and 298 ethnically matched controls. All individuals were HIV negative. Decreased risk of pulmonary TB was associated with both heterozygosity and homozygosity for the IL1B-511-C allele (OR 0.66, P = 0.027, and OR 0.58, P = 0.015, respectively). Nonetheless, the C allele was present at a frequency of 0.66 in TB cases suggesting that whilst IL-1beta contributes to disease susceptibility, it is not the major factor. There was no association between the IL1B+3953-T/C polymorphism or the 86 bp IL1RN pentallelic repeat and TB in this population. Using an ex-vivo whole blood assay, healthy Gambian individuals who are homozygous for the IL1B-511-T allele failed to exhibit a significant increase in IL-1beta production in response to LPS after IFN-gamma priming.