Mapping of PARK2 and PACRG Overlapping Regulatory Region Reveals LD Structure and Functional Variants in Association with Leprosy in Unrelated Indian Population Groups

Leprosy is a chronic infectious disease caused by Mycobacterium Leprae, where the host genetic background plays an important role toward the disease pathogenesis. Various studies have identified a number of human genes in association with leprosy or its clinical forms. However, non-replication of results has hinted at the heterogeneity among associations between different population groups, which could be due to differently evolved LD structures and differential frequencies of SNPs within the studied regions of the genome. A need for systematic and saturated mapping of the associated regions with the disease is warranted to unravel the observed heterogeneity in different populations. Mapping of the PARK2 and PACRG gene regulatory region with 96 SNPs, with a resolution of 1 SNP per 1 Kb for PARK2 gene regulatory region in a North Indian population, showed an involvement of 11 SNPs in determining the susceptibility towards leprosy. The association was replicated in a geographically distinct and unrelated population from Orissa in eastern India. In vitro reporter assays revealed that the two significantly associated SNPs, located 63.8 kb upstream of PARK2 gene and represented in a single BIN of 8 SNPs, influenced the gene expression. A comparison of BINs between Indian and Vietnamese populations revealed differences in the BIN structures, explaining the heterogeneity and also the reason for non-replication of the associated genomic region in different populations.     

Allelic variation in BTNL2 and susceptibility to tuberculosis in a South African population

Tuberculosis and sarcoidosis show phenotypic features of granulomatous disease. The bacterium Mycobacterium tuberculosis can induce the expression of the sarcoidosis susceptibility gene BTNL2 in monocyte-derived macrophages. BTNL2 was therefore investigated as a candidate gene for tuberculosis in a case-control association study in the South African Coloured population. We sequenced the coding regions of BTNL2 to detect known and novel polymorphisms and genotyped 18 SNPs in 432 pulmonary tuberculosis cases and 482 controls. We did not find a significant association between the truncating rs2076530 SNP, previously associated with sarcoidosis, and tuberculosis. No association was found between any of the other SNPs studied and disease and none of the estimated haplotypes showed any association with TB. Comparative analyses with the South African data from this study and published data on German and American populations revealed that, for a segment of BTNL2, the admixed, but not stratified, South African population resembles the African-Americans more than white populations.     

<Emphasis Type="Italic">TNF</Emphasis><Emphasis Type="Italic">α</Emphasis> and <Emphasis Type="Italic">LT</Emphasis><Emphasis Type="Italic">α</Emphasis> gene polymorphisms as additional markers of celiac disease susceptibility in a DQ2-positive population

TNFalpha and TNFbeta, or linfotoxin (LTalpha), are two molecules playing an important role in inflammation. Their genes map on Chromosome 6, between the HLA class II and class I loci. Polymorphisms in, or near, TNF genes have been associated with susceptibility to several autoimmune diseases. Studies of TNF genes in celiac disease (CD) have presented contradictory results. We have assessed the role of TNFalpha and linfotoxin alpha (TNFbeta) in CD and their relative value as CD markers in addition to the presence of DQ2. The TNFA -308 polymorphism and the polymorphism at the first intron of the LTA gene were typed in CD patients and healthy controls and the results were correlated with the presence of DQ2. Significant differences were found in genotype and allele frequencies for the TNFA and LTA genes between CD patients and controls, with an increase in the presence of the TNFA*2 and LTA*1 alleles in CD patients. These differences increase when DQ2-positive CD patients and DQ2-positive controls are compared. In DQ2-positive individuals, allele 2 (A) in position -308 of the promoter of TNFA and allele 1 (G) of the NcoI RFLP in the first intron of LTA are additional risk markers for CD.     

A second susceptibility gene for developing rheumatoid arthritis in the human MHC is localized within a 70-kb interval telomeric of the TNF genes in the HLA class III region

Rheumatoid arthritis (RA) is a chronic inflammatory joint disease with a multifactorial genetic basis. However, pathogenic genes for RA other than the human leukocyte antigen (HLA)-DRB1 gene have yet to be identified. Here, we investigated whether there is a second susceptibility locus for RA within the human major histocompatibility complex using 18 microsatellite markers distributed from the centromeric (HSET) to the telomeric end (P5-15) of the 3.6-Mb HLA region. Statistical studies of associated alleles on each microsatellite locus showed that one pathogenic gene for RA in the HLA region is localized in the DRB1 gene, as expected. Further, a second susceptibility gene of RA was suggested to be present in the HLA class III region, narrowed to 70 kb, that is just telomeric of the TNF gene cluster (TNFA and LTA) and that is located between the microsatellites TNFa and C1-2-A. In this critical segment, four expressed genes have been thus far identified, NFKBIL1 (IkappaBL), ATP6G, BAT1, and MICB, all of which are candidate genes for determining susceptibility to RA. These results exclude the possibility of involvement of the TNFA genes (TNF-alpha) in the development of RA, which was suggested previously to be a strong candidate for RA in the class III region.     

The BTNL2 gene and sarcoidosis susceptibility in African Americans and Whites

The BTNL2 gene is a member of the B7 receptor family that probably functions as a T-cell costimulatory molecule. It resides in the class II major histocompatibility complex (MHC) region of chromosome 6p and has recently been associated with sarcoidosis susceptibility in a white German population. We sought to replicate the BTNL2 association in an African American family-based study population (n=219 nuclear families) and two case-control populations--one African American (n=295 pairs) and one white (n=366 pairs). Ten SNPs were detected within a 490-bp region spanning exon/intron 5 of BTNL2. Haplotype variation within this region was significantly associated with sarcoidosis in all three study populations but more so in whites (P=.0006) than in the African American case-control (P=.02) or family-based (P=.03) samples. The previously reported BTNL2 SNP with the strongest sarcoidosis association, rs2076530, was also the SNP with the strongest association in our white population (P<.0001). The A allele of rs2076530 results in a premature exon-splice site and increases risk for sarcoidosis (odds ratio=2.03; 95% confidence interval 1.32-3.12). Although rs2076530 was not associated with sarcoidosis in either African American sample, a three-locus haplotype that included rs2076530 was associated with sarcoidosis across all three study samples. Multivariable logistic regression analyses showed that BTNL2 effects are independent of human leukocyte antigen class II genes in whites but may interact antagonistically in African Americans. Our results underscore the complexity of genetic risk for sarcoidosis emanating from the MHC region.     

Characterization of Single-Nucleotide Polymorphisms in the Tumor Necrosis Factor α Promoter Region and in Lymphotoxin α in Squamous Intraepithelial Lesions, Precursors of Cervical Cancer

Development of cervical cancer is a long process of abnormal cancerous cell growth in the cervix and is primarily the result of infection with specific high-risk types of human papillomavirus (HPV). The cytokines tumor necrosis factor α (TNFα) and lymphotoxin α (LTA) have an important role in all stages of cervical cancer and have the ability to induce the regression or promote the development of human tumors. Biologically important single-nucleotide polymorphisms (SNPs) occur within the TNFα and LTA genes. Therefore, the purpose of this study was to investigate the SNPs in the TNFα promoter region (-163, -238, -244, -308, -376, -857, -863, and -1031) and in the first intron of LTA (+252) in women with precursor lesions of cervical cancer. Overall, we studied 396 women from Mexico City. A total of 191 patients with HPV infection and precursor cervical lesions were subdivided in two groups: those with low-grade squamous intraepithelial lesions (n = 132) and those with high-grade squamous intraepithelial lesions (n = 59). Women (n = 205) negative for HPV and without cervical lesions were also included in the study. DNA was extracted from peripheral white blood cells and from cervical samples, and detection of biallelic polymorphisms of TNFα and LTA was performed using the polymerase chain reaction-sequence-specific oligonucleotide probe and restriction fragment length polymorphism techniques, respectively. We demonstrated that risk is associated with the genotype G/A (odds ratio = 2.48) and that protection is associated with the genotype G/G of SNP TNFα -376 (odds ratio = 0.37).     

Genetics of sputum gene expression in chronic obstructive pulmonary disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.     

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.     

Selecting SNPs for association studies based on population frequencies: a novel interactive tool and its application to polygenic diseases

Common complex polygenic diseases as autoimmune diseases have not been completely understood on a molecular level. While many genes are known to be involved in the pathways responsible for the phenotype, explicit causes for the susceptibility of the disease remain to be elucidated. The susceptibility to disease is thought to be the result of genetic epistatic interactions between common polymorphic genes. This polymorphism is mostly caused by single nucleotide polymorphisms (SNPs). Human subpopulations are known to differ in the susceptibility to the diseases and generally in the distribution of single nucleotide polymorphisms. The here presented approach retrieves SNPs with the most divergent frequencies for selected human subpopulations to help defining properties for the experimental verification of SNPs within defined regions. A web-accessible program implementing this approach was evaluated for multiple sclerosis (MS), a common human polygenic disease. A link to a summary of data from "The SNP Consortium" (TSC) with sex-dependencies of SNPs is available. Associations of SNPs to genes, genetic markers and chromosomal loci are retrieved from the Ensembl project. This tool is recommended to be used in conjunction with microarray analyses or marker association studies that link genes or chromosomal loci to particular diseases.     

Linkage disequilibrium pattern and age-at-diagnosis are critical for replicating genetic associations across ethnic groups in leprosy

One of the persistent challenges of genetic association studies is the replication of genetic marker-disease associations across ethnic groups. Here, we conducted high-density association mapping of PARK2/PACRG SNPs with leprosy and identified 69 SNPs significantly associated with leprosy in 198 single-case Vietnamese leprosy families. A total of 56 associated SNPs localized to the overlapping promoter regions of PARK2/PACRG. For this region, multivariate analysis identified four SNPs belonging to two major SNP bins (rs1333955, rs7744433) and two single SNP bins (rs2023004, rs6936895) that capture the combined statistical evidence (P = 1.1 × 10^?5) for association among Vietnamese patients. Next, we enrolled a case–control sample of 364 leprosy cases and 370 controls from Northern India. We genotyped all subjects for 149 SNPs that capture >80 % of the genetic variation in the Vietnamese sample and found 24 SNPs significantly associated with leprosy. Multivariate analysis identified three SNPs (rs1333955, rs9356058 and rs2023004) that capture the association with leprosy (P < 10^?8). Hence, two SNPs (rs1333955 and rs2023004) were replicated by multivariate analysis between both ethnic groups. Marked differences in the linkage disequilibrium pattern explained some of the differences in univariate analysis between the two ethnic groups. In addition, the strength of association for two promoter region SNP bins was significantly stronger among young leprosy patients in the Vietnamese sample. The same trend was observed in the Indian sample, but due to the higher age-at-diagnosis of the patients the age effect was less pronounced.     

Discovery and genetic mapping of single nucleotide polymorphisms in candidate genes for pathogen defence response in perennial ryegrass (<Emphasis Type="Italic">Lolium perenne</Emphasis> L.)

Susceptibility to foliar pathogens commonly causes significant reductions in productivity of the important temperate forage perennial ryegrass. Breeding for durable disease resistance involves not only the deployment of major genes but also the additive effects of minor genes. An approach based on in vitro single nucleotide polymorphism (SNP) discovery in candidate defence response (DR) genes has been used to develop potential diagnostic genetic markers. SNPs were predicted, validated and mapped for representatives of the pathogenesis-related (PR) protein-encoding and reactive oxygen species (ROS)-generating gene classes. The F(1)(NA(6) x AU(6)) two-way pseudo-test cross population was used for SNP genetic mapping and detection of quantitative trait loci (QTLs) in response to a crown rust field infection. Novel resistance QTLs were coincident with mapped DR gene SNPs. QTLs on LG3 and LG7 also coincided with both herbage quality QTLs and candidate genes for lignin biosynthesis. Multiple DR gene SNP loci additionally co-located with QTLs for grey leaf spot, bacterial wilt and crown rust resistance from other published studies. Further functional validation of DR gene SNP loci using methods such as fine-mapping and association genetics will improve the efficiency of parental selection based on superior allele content.     

Deep resequencing identifies candidate functional genes in leprosy GWAS loci

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.     

An investigation of genome-wide studies reported susceptibility loci for ulcerative colitis shows limited replication in north Indians

Genome-Wide Association studies (GWAS) of both Crohn's Disease (CD) and Ulcerative Colitis (UC) have unearthed over 40 risk conferring variants. Recently, a meta-analysis on UC revealed several loci, most of which were either previously associated with UC or CD susceptibility in populations of European origin. In this study, we attempted to replicate these findings in an ethnically distinct north Indian UC cohort. 648 UC cases and 850 controls were genotyped using Infinium Human 660W-quad. Out of 59 meta-analysis index SNPs, six were not in the SNP array used in the study. Of the remaining 53 SNPs, four were found monomorphic. Association (p<0.05) at 25 SNPs was observed, of which 15 were CD specific. Only five SNPs namely rs2395185 (HLA-DRA), rs3024505 (IL10), rs6426833 (RNF186), rs3763313 (BTNL2) and rs2066843 (NOD2) retained significance after Bonferroni correction. These results (i) reveal limited replication of Caucasian based meta-analysis results; (ii) reiterate overlapping molecular mechanism(s) in UC and CD; (iii) indicate differences in genetic architecture between populations; and (iv) suggest that resources such as HapMap need to be extended to cover diverse ethnic populations. They also suggest a systematic GWAS in this terrain may be insightful for identifying population specific IBD risk conferring loci and thus enable cross-ethnicity fine mapping of disease loci.     

<Emphasis Type="Italic">LIG1</Emphasis> polymorphisms: the Indian scenario

Elucidation of the genetic diversity and relatedness of the subpopulations of India may provide a unique resource for future analysis of genetic association of several critical community-specific complex diseases. We performed a comprehensive exploration of single nucleotide polymorphisms (SNPs) within the gene DNA ligase 1 (LIG1) among a multiethnic panel of Indian subpopulations representative of the ethnic, linguistic and geographical diversity of India using a two-stage design involving DNA resequencing-based SNP discovery followed by SNP validation using sequenom-based genotyping. Thirty SNPs were identified in LIG1 gene using DNA resequencing including three promoter SNPs and one coding SNP. Following SNP validation, the SNPs rs20580/C19008A and rs3730862/C8804T were found to have the most widespread prevalence with noticeable variations in minor allele frequencies both between the Indian subpopulation groups and also from those reported on other major world populations. Subsequently, SNPs found in Indian subpopulations were analysed using bioinformatics-based approaches and compared with SNP data available on major world populations. Further, we also performed genotype–phenotype association analysis of LIG1 SNPs with publicly available data on LIG1 mRNA expression in HapMap samples. Results showed polymorphisms in LIG1 affect its expression and may therefore change its function. Our results stress upon the uniqueness of the Indian population with respect to the worldwide scenario and suggest that any epidemiological study undertaken on the global population should take this distinctiveness in consideration and avoid making generalized conclusions.     

Regulatory polymorphisms in the cyclophilin A gene, PPIA, accelerate progression to AIDS

Human cyclophilin A, or CypA, encoded by the gene peptidyl prolyl isomerase A (PPIA), is incorporated into the HIV type 1 (HIV-1) virion and promotes HIV-1 infectivity by facilitating virus uncoating. We examined the effect of single nucleotide polymorphisms (SNPs) and haplotypes within the PPIA gene on HIV-1 infection and disease progression in five HIV-1 longitudinal history cohorts. Kaplan-Meier survival statistics and Cox proportional hazards model were used to assess time to AIDS outcomes. Among eight SNPs tested, two promoter SNPs (SNP3 and SNP4) in perfect linkage disequilibrium were associated with more rapid CD4(+) T-cell loss (relative hazard = 3.7, p = 0.003) in African Americans. Among European Americans, these alleles were also associated with a significant trend to more rapid progression to AIDS in a multi-point categorical analysis (p = 0.005). Both SNPs showed differential nuclear protein-binding efficiencies in a gel shift assay. In addition, one SNP (SNP5) located in the 5' UTR previously shown to be associated with higher ex vivo HIV-1 replication was found to be more frequent in HIV-1-positive individuals than in those highly exposed uninfected individuals. These results implicate regulatory PPIA polymorphisms as a component of genetic susceptibility to HIV-1 infection or disease progression, affirming the important role of PPIA in HIV-1 pathogenesis.