RESUMEN
Previous genotyping-based assays have identified non-coding variants of several interleukins (ILs) being associated with genetic susceptibility to leprosy. However, understanding of the involvement of coding variants within all IL family genes in leprosy was still limited. To obtain the full mutation spectrum of all ILs in leprosy, we performed a targeted deep sequencing of coding regions of 58 ILs genes in 798 leprosy patients (age 56.2 ± 14.4; female 31.5%) and 990 healthy controls (age 38.1 ± 14.0; female 44.3%) from Yunnan, Southwest China. mRNA expression alterations of ILs in leprosy skin lesions or in response to M. leprae treatment were estimated by using publicly available expression datasets. Two coding variants in IL27 (rs17855750, p.S59A, p = 4.02 × 10-8 , odds ratio [OR] = 1.748) and IL1RN (rs45507693, p.A106T, p = 1.45 × 10-5 , OR = 3.629) were significantly associated with leprosy risk. mRNA levels of IL27 and IL1RN were upregulated in whole blood cells after M. leprae stimulation. These data showed that IL27 and IL1RN are leprosy risk genes. Further functional study is required for characterizing the exact role of ILs in leprosy.
Asunto(s)
Predisposición Genética a la Enfermedad/genética , Interleucinas/genética , Lepra/genética , Polimorfismo de Nucleótido Simple/genética , Pueblo Asiatico/genética , Estudios de Casos y Controles , China , Femenino , Humanos , Masculino , Persona de Mediana Edad , ARN Mensajero/genéticaRESUMEN
Genome-wide association studies (GWASs) and genome-wide linkage studies (GWLSs) have identified numerous risk genes affecting the susceptibility to leprosy. However, most of the reported GWAS hits are noncoding variants and account for only part of the estimated heritability for this disease. In order to identify additional risk genes and map the potentially functional variants within the GWAS loci, we performed a three-stage study combining whole-exome sequencing (WES; discovery stage), targeted next-generation sequencing (NGS; screening stage), and refined validation of risk missense variants in 1,433 individuals with leprosy and 1,625 healthy control individuals from Yunnan Province, Southwest China. We identified and validated a rare damaging variant, rs142179458 (c.1045G>A [p.Asp349Asn]) in HIF1A, as contributing to leprosy risk (p = 4.95 × 10-9, odds ratio [OR] = 2.266). We were able to show that affected individuals harboring the risk allele presented with multibacillary leprosy at an earlier age (p = 0.025). We also confirmed the association between missense variant rs3764147 (c.760A>G [p.Ile254Val]) in the GWAS hit LACC1 (formerly C13orf31) and leprosy (p = 6.11 × 10-18, OR = 1.605). By using the population attributable fraction, we have shown that HIF1A and LACC1 are the major genes with missense variants contributing to leprosy risk in our study groups. Consistently, mRNA expression levels of both HIF1A and LACC1 were upregulated in the skin lesions of individuals with leprosy and in Mycobacterium leprae-stimulated cells, indicating an active role of HIF1A and LACC1 in leprosy pathogenesis.
Asunto(s)
Pueblo Asiatico/genética , Etnicidad/genética , Predisposición Genética a la Enfermedad , Subunidad alfa del Factor 1 Inducible por Hipoxia/genética , Lepra/genética , Mutación Missense/genética , Proteínas/genética , Adolescente , Adulto , Anciano , Anciano de 80 o más Años , Alelos , Niño , Estudios de Cohortes , Femenino , Secuenciación de Nucleótidos de Alto Rendimiento , Humanos , Péptidos y Proteínas de Señalización Intracelular , Masculino , Persona de Mediana Edad , Proteínas Nucleares/genética , Polimorfismo de Nucleótido Simple/genética , ARN Mensajero/genética , ARN Mensajero/metabolismo , Reproducibilidad de los Resultados , Factores de Riesgo , Transactivadores/genética , Regulación hacia Arriba/genética , Secuenciación del Exoma , Adulto JovenRESUMEN
BACKGROUND: The pathogen Mycobacterium leprae of leprosy is heavily dependent on the host energy metabolites and nutritional products for survival. Previously we and others have identified associations of several mitochondrion-related genes and mitochondrial DNA (mtDNA) copy number alterations with leprosy and/or its subtype. We hypothesized that genetic variants of mtDNA replication-related genes would affect leprosy. OBJECTIVE: We aimed to identify genetic associations between the mtDNA replication-related genes TFAM, POLG and leprosy. METHODS: Genetic association study was performed in 2898 individuals from two independent sample sets in Yunnan Province, China. We first screened 7 tag SNPs of TFAM and POLG in 527 leprosy cases and 583 controls (Sample I). Expression quantitative trait loci (eQTL) analysis and differential mRNA expression were analyzed to discern potential effect of risk variants. The entire exon region of TFAM and POLG were further analyzed in 798 leprosy cases and 990 controls (Sample II; 4327 East Asians from the ExAC dataset was included as a reference control) by using targeted gene sequencing for fine mapping potentially causal variants. RESULTS: Two tag SNPs of TFAM (rs1049432, P=0.007) and POLG (rs3176238, P=0.006) were associated with multibacillary leprosy (MB) in Sample I and the significance survived correction for multiple comparisons. SNPs rs1937 of TFAM (which was linked with rs1049432) and rs61756401 of POLG were associated with leprosy, whereas no potentially causative coding variants were identified in Sample II. The eQTL analysis showed that rs1049432 was a significant cis eQTL for TFAM in nerve tissue (P=1.20×10-12), and rs3176238 was a significant cis eQTL for POLG in nerve (P=3.90×10-13) and skin tissues (P=2.50×10-11). Consistently, mRNA level of POLG was differentially expressed in leprotic skin lesions. CONCLUSIONS: Genetic variants of TFAM and POLG were associated with leprosy in Han Chinese, presumably by affecting gene expression.