Toll-like receptor 1 N248S single-nucleotide polymorphism is associated with leprosy risk and regulates immune activation during mycobacterial infection

Conflicting findings about the association between leprosy and TLR1 variants N248S and I602S have been reported. Here, we performed case-control and family based studies, followed by replication in 2 case-control populations from Brazil, involving 3162 individuals. Results indicated an association between TLR1 248S and leprosy in the case-control study (SS genotype odds ratio [OR], 1.81; P = .004) and the family based study (z = 2.02; P = .05). This association was consistently replicated in other populations (combined OR, 1.51; P < .001), corroborating the finding that 248S is a susceptibility factor for leprosy. Additionally, we demonstrated that peripheral blood mononuclear cells (PBMCs) carrying 248S produce a lower tumor necrosis factor/interleukin-10 ratio when stimulated with Mycobacterium leprae but not with lipopolysaccharide or PAM3cysK4. The same effect was observed after infection of PBMCs with the Moreau strain of bacillus Calmette-Guerin but not after infection with other strains. Finally, molecular dynamics simulations indicated that the Toll-like receptor 1 structure containing 248S amino acid is different from the structure containing 248N. Our results suggest that TLR1 248S is associated with an increased risk for leprosy, consistent with its hypoimmune regulatory function.

The common G-allele of interleukin-18 single-nucleotide polymorphism is a genetic risk factor for atopic asthma. The SAPALDIA Cohort Study

BACKGROUND: IL-18 is a pleiotrophic cytokine involved in both, T-helper type 1 (Th1) and Th2 differentiation. Recently genetic variants in the IL-18 gene have been associated with increased risk of atopy and asthma. OBJECTIVE: To examine the relationship of a genetic, haplotype-tagging promotor variant -137G/C in the IL-18 gene with atopic asthma in a large, well-characterized and population-based study of adults. METHODS: Prospective cohort study design was used to collect interview and biological measurement data at two examination time-points 11 years apart. Multivariate logistic regression analysis was used to assess the association of genotype with asthma and atopy. RESULTS: The G-allele of the IL-18 promotor variant (-137G/C) was associated with a markedly increased risk for the prevalence of physician-diagnosed asthma with concomitant skin reactivity to common allergens. Stratification of the asthma cases by skin reactivity to common allergens revealed an exclusive association of IL-18 -137 G-allele with an increased prevalence of atopic asthma (adjusted odds ratio (OR): 3.63; 95% confidence interval: (1.64-8.02) for GC or GG carriers vs. CC carriers), and no according association with asthma and concomitant negative skin reactivity (adjusted OR: 1.13; 0.66-1.94). The interaction between IL-18 -137G/C genotype and positive skin prick test was statistically significant (P=0.029). None of 74 incident asthma cases with atopy at baseline exhibited the CC genotype. CONCLUSION: Our results strongly suggest that this variant of the IL-18 gene is an important genetic determinant involved in the development of atopic asthma.

The HCP5 single-nucleotide polymorphism: a simple screening tool for prediction of hypersensitivity reaction to abacavir

The HLA-B 5701 allele is predictive of hypersensitivity reaction to abacavir, a response herein termed "ABC-HSR." This study of 1,103 individuals infected with human immunodeficiency virus assessed the usefulness of genotyping a HCP5 single-nucleotide polymorphism (SNP), rs2395029, in relation to ABC-HSR. In populations with European ancestry, rs2395029 is in linkage disequilibrium with HLA-B 5701. The HCP5 SNP was present in all 98 HLA-B 5701-positive individuals and was absent in 999 of 1005 HLA-B 5701-negative individuals. rs2395029 was overrepresented in 25 individuals with clinically likely ABC-HSR, compared with its frequency in 175 ABC-tolerant individuals (80% vs. 2%, respectively; P < .0001). Therefore, HCP5 genotyping could serve as a simple screening tool for ABC-HSR, particularly in settings where sequence-based HLA typing is not available.

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

Alkaline-mediated differential interaction (AMDI): A simple automatable single-nucleotide polymorphism assay

The key requirements for high-throughput single-nucleotide polymorphism (SNP) typing of DNA samples in large-scale disease case-control studies are automatability, simplicity, and robustness, coupled with minimal cost. In this paper we describe a fluorescence technique for the detection of SNPs that have been amplified by using the amplification refractory mutation system (ARMS)-PCR procedure. Its performance was evaluated using 32 sequence-specific primer mixes to assign the HLA-DRB alleles to 80 lymphoblastoid cell line DNAs chosen from our database for their diversity. All had been typed previously by alternative methods, either direct sequencing or gel electrophoresis. We believe the detection system that we call AMDI (alkaline-mediated differential interaction) satisfies the above criteria and is suitable for general high-throughput SNP typing.

A single nucleotide polymorphism in the RAD51 gene modifies cancer risk in BRCA2 but not BRCA1 carriers

BRCA1 and BRCA2 carriers are at increased risk for both breast and ovarian cancer, but estimates of lifetime risk vary widely, suggesting their penetrance is modified by other genetic and/or environmental factors. The BRCA1 and BRCA2 proteins function in DNA repair in conjunction with RAD51. A preliminary report suggested that a single nucleotide polymorphism in the 5′ untranslated region of RAD51 (135C/G) increases breast cancer risk in BRCA1 and BRCA2 carriers. To investigate this effect we studied 257 female Ashkenazi Jewish carriers of one of the common BRCA1 (185delAG, 5382insC) or BRCA2 (6174delT) mutations. Of this group, 164 were affected with breast and/or ovarian cancer and 93 were unaffected. RAD51 genotyping was performed on all subjects. Among BRCA1 carriers, RAD51-135C frequency was similar in healthy and affected women [6.1% (3 of 49) and 9.9% (12 of 121), respectively], and RAD-135C did not influence age of cancer diagnosis [Hazard ratio (HR) = 1.18 for disease in RAD51-135C heterozygotes, not significant]. However, in BRCA2 carriers, RAD51-135C heterozygote frequency in affected women was 17.4% (8 of 46) compared with 4.9% (2 of 41) in unaffected women (P = 0.07). Survival analysis in BRCA2 carriers showed RAD51-135C increased risk of breast and/or ovarian cancer with an HR of 4.0 [95% confidence interval 1.6–9.8, P = 0.003]. This effect was largely due to increased breast cancer risk with an HR of 3.46 (95% confidence interval 1.3–9.2, P = 0.01) for breast cancer in BRCA2 carriers who were RAD51-135C heterozygotes. RAD51 status did not affect ovarian cancer risk. These results show RAD51-135C is a clinically significant modifier of BRCA2 penetrance, specifically in raising breast cancer risk at younger ages.

Development of single nucleotide polymorphism (SNP) markers from the mango (Mangifera indica) transcriptome for mapping and estimation of genetic diversity

The development of molecular markers for genomic studies in Mangifera indica (mango) will allow marker-assisted selection and identification of genetically diverse germplasm, greatly aiding mango breeding programs. We report here our identification of thousands of unambiguous molecular markers that can be easily assayed across genotypes of the species. With origin centered in Southeast Asia, mangos are grown throughout the tropics and subtropics as a nutritious fruit that exhibits remarkable intraspecific phenotypic diversity. With the goal of building a high density genetic map, we have undertaken discovery of sequence variation in expressed genes across a broad range of mango cultivars. A transcriptome sequence reference was built de novo from extensive sequencing and assembly of RNA from cultivar 'Tommy Atkins'. Single nucleotide polymorphisms (SNPs) in protein coding transcripts were determined from alignment of RNA reads from 24 mango cultivars of diverse origins: 'Amin Abrahimpur' (India), 'Aroemanis' (Indonesia), 'Burma' (Burma), 'CAC' (Hawaii), 'Duncan' (Florida), 'Edward' (Florida), 'Everbearing' (Florida), 'Gary' (Florida), 'Hodson' (Florida), 'Itamaraca' (Brazil), 'Jakarata' (Florida), 'Long' (Jamaica), 'M. Casturi Purple' (Borneo), 'Malindi' (Kenya), 'Mulgoba' (India), 'Neelum' (India), 'Peach' (unknown), 'Prieto' (Cuba), 'Sandersha' (India), 'Tete Nene' (Puerto Rico), 'Thai Everbearing' (Thailand), 'Toledo' (Cuba), 'Tommy Atkins' (Florida) and 'Turpentine' (West Indies). SNPs in a selected subset of protein coding transcripts are currently being converted into Fluidigm assays for genotyping of mapping populations and germplasm collections. Using an alternate approach, SNPs (144) discovered by sequencing of candidate genes in 'Kensington Pride' have already been converted and used for genotyping.

A molecular platform for the diagnosis of multidrug-resistant and pre-extensively drug-resistant tuberculosis based on single nucleotide polymorphism mutations present in Colombian isolates of Mycobacterium tuberculosis

Developing a fast, inexpensive, and specific test that reflects the mutations present in Mycobacterium tuberculosis isolates according to geographic region is the main challenge for drug-resistant tuberculosis (TB) control. The objective of this study was to develop a molecular platform to make a rapid diagnosis of multidrug-resistant (MDR) and extensively drug-resistant TB based on single nucleotide polymorphism (SNP) mutations present in the rpoB, katG, inhA, ahpC, and gyrA genes from Colombian M. tuberculosis isolates. The amplification and sequencing of each target gene was performed. Capture oligonucleotides, which were tested before being used with isolates to assess the performance, were designed for wild type and mutated codons, and the platform was standardised based on the reverse hybridisation principle. This method was tested on DNA samples extracted from clinical isolates from 160 Colombian patients who were previously phenotypically and genotypically characterised as having susceptible or MDR M. tuberculosis. For our method, the kappa index of the sequencing results was 0,966, 0,825, 0,766, 0,740, and 0,625 for rpoB, katG, inhA, ahpC, and gyrA, respectively. Sensitivity and specificity were ranked between 90-100% compared with those of phenotypic drug susceptibility testing. Our assay helps to pave the way for implementation locally and for specifically adapted methods that can simultaneously detect drug resistance mutations to first and second-line drugs within a few hours.

Single-nucleotide polymorphism in two representative multidrug-resistant Mycobacterium bovis isolates collected from patients in a Spanish hospital harboring a human infection outbreak

Mycobacterium bovis is the etiological agent of tuberculosis in domestic and wild animals. Its involvement as a human pathogen has been highlighted again with the recent descriptions of transmission through dairy products (18), reactivation or primary infection in human immunodeficiency virus-infected patients (5), and association with meat industry workers, animal keepers, or hunters (3). Strains resistant to antituberculous drugs (M. bovis is naturally resistant to pyrazinamide) pose an additional risk (2). Several studies have demonstrated that mutations in target genes are associated with resistance to antituberculous drugs (4, 7, 10, 11, 16). However, most of them have been developed in Mycobacterium tuberculosis strains and limited data are available regarding M. bovis isolates. The aim of this study was to characterize by sequencing the main genes involved in antibiotic resistance in two multidrug-resistant (MDR) M. bovis isolates in a human outbreak detected in a hospital in Madrid that subsequently spread to several countries (5, 6, 15). The isolates were resistant to 11 drugs, but only their rpoB and katG genes have been analyzed so far (1, 14). We studied the first (93/R1) and last (95/R4) M. bovis isolates of this nosocomial outbreak, characterized by spoligotyping as SB0426 (hexacode 63-5F-5E-7F-FF-60 in the database at www.mbovis.org) (1, 13). Several genes involved in resistance to isoniazid (katG, ahpC, inhA, and the oxyR-ahpC intergenic region), rifampin (rpoB), streptomycin (rrs, rpsL), ethambutol (embB), and quinolones (gyrA) were studied. These genes, or fragments of genes, were amplified and sequenced as previously described (12). The sequence analysis revealed polymorphisms in five (ahpC, rpoB, rpsL, embB, and gyrA) out of nine analyzed genes (Table 1). Nucleotide substitutions in four genes cause a change in the encoded amino acid. Two additional synonymous mutations in ahpC and rpsL differentiated the first and last isolates from the outbreak.

A single nucleotide polymorphism in NEUROD1 is associated with production traits in nelore beef cattle.

Feed efficiency and carcass characteristics are late-measured traits. The detection of molecular markers associated with them can help breeding programs to select animals early in life, and to predict breeding values with high accuracy. The objective of this study was to identify polymorphisms in the functional and positional candidate gene NEUROD1 (neurogenic differentiation 1), and investigate their associations with production traits in reference families of Nelore cattle. A total of 585 steers were used, from 34 sires chosen to represent the variability of this breed. By sequencing 14 animals with extreme residual feed intake (RFI) values, seven single nucleotide polymorphisms (SNPs) in NEUROD1 were identified. The investigation of marker effects on the target traits RFI, backfat thickness (BFT), ribeye area (REA), average body weight (ABW), and metabolic body weight (MBW) was performed with a mixed model using the restricted maximum likelihood method. SNP1062, which changes cytosine for guanine, had no significant association with RFI or REA. However, we found an additive effect on ABW (P ≤ 0.05) and MBW (P ≤ 0.05), with an estimated allele substitution effect of -1.59 and -0.93 kg0.75, respectively. A dominant effect of this SNP for BFT was also found (P ≤ 0.010). Our results are the first that identify NEUROD1 as a candidate that affects BFT, ABW, and MBW. Once confirmed, the inclusion of this SNP in dense panels may improve the accuracy of genomic selection for these traits in Nelore beef cattle as this SNP is not currently represented on SNP chips.

The association of single nucleotide polymorphisms with endometrial cancer risk and prognosis

Endometrial cancer is one of the most common female diseases in developed nations and is the most commonly diagnosed gynaecological cancer in Australia. The disease is commonly classified by histology: endometrioid or non-endometrioid endometrial cancer. While non-endometrioid endometrial cancers are accepted to be high-grade, aggressive cancers, endometrioid cancers (comprising 80% of all endometrial cancers diagnosed) generally carry a favourable patient prognosis. However, endometrioid endometrial cancer patients endure significant morbidity due to surgery and radiotherapy used for disease treatment, and patients with recurrent disease have a 5-year survival rate of less than 50%. Genetic analysis of women with endometrial cancer could uncover novel markers associated with disease risk and/or prognosis, which could then be used to identify women at high risk and for the use of specialised treatments. Proteases are widely accepted to play an important role in the development and progression of cancer. This PhD project hypothesised that SNPs from two protease gene families, the matrix metalloproteases (MMPs, including their tissue inhibitors, TIMPs) and the tissue kallikrein-related peptidases (KLKs) would be associated with endometrial cancer susceptibility and/or prognosis. In the first part of this study, optimisation of the genotyping techniques was performed. Results from previously published endometrial cancer genetic association studies were attempted to be validated in a large, multicentre replication set (maximum cases n = 2,888, controls n = 4,483, 3 studies). The rs11224561 progesterone receptor SNP (PGR, A/G) was observed to be associated with increased endometrial cancer risk (per A allele OR 1.31, 95% CI 1.12-1.53; p-trend = 0.001), a result which was initially reported among a Chinese sample set. Previously reported associations for the remaining 8 SNPs investigated for this section of the PhD study were not confirmed, thereby reinforcing the importance of validation of genetic association studies. To examine the effect of SNPs from the MMP and KLK families on endometrial cancer risk, we selected the most significantly associated MMP and KLK SNPs from genome-wide association study analysis (GWAS) to be genotyped in the GWAS replication set (cases n = 4,725, controls n = 9,803, 13 studies). The significance of the MMP24 rs932562 SNP was unchanged after incorporation of the stage 2 samples (Stage 1 per allele OR 1.18, p = 0.002; Combined Stage 1 and 2 OR 1.09, p = 0.002). The rs10426 SNP, located 3' to KLK10 was predicted by bioinformatic analysis to effect miRNA binding. This SNP was observed in the GWAS stage 1 result to exhibit a recessive effect on endometrial cancer risk, a result which was not validated in the stage 2 sample set (Stage 1 OR 1.44, p = 0.007; Combined Stage 1 and 2 OR 1.14, p = 0.08). Investigation of the regions imputed surrounding the MMP, TIMP and KLK genes did not reveal any significant targets for further analysis. Analysis of the case data from the endometrial cancer GWAS to identify genetic variation associated with cancer grade did not reveal SNPs from the MMP, TIMP or KLK genes to be statistically significant. However, the representation of SNPs from the MMP, TIMP and KLK families by the GWAS genotyping platform used in this PhD project was examined and observed to be very low, with the genetic variation of four genes (MMP23A, MMP23B, MMP28 and TIMP1) not captured at all by this technique. This suggests that comprehensive candidate gene association studies will be required to assess the role of SNPs from these genes with endometrial cancer risk and prognosis. Meta-analysis of gene expression microarray datasets curated as part of this PhD study identified a number of MMP, TIMP and KLK genes to display differential expression by endometrial cancer status (MMP2, MMP10, MMP11, MMP13, MMP19, MMP25 and KLK1) and histology (MMP2, MMP11, MMP12, MMP26, MMP28, TIMP2, TIMP3, KLK6, KLK7, KLK11 and KLK12). In light of these findings these genes should be prioritised for future targeted genetic association studies. Two SNPs located 43.5 Mb apart on chromosome 15 were observed from the GWAS analysis to be associated with increased endometrial cancer grade, results that were validated in silico in two independent datasets. One of these SNPs, rs8035725 is located in the 5' untranslated region of a MYC promoter binding protein DENND4A (Stage 1 OR 1.15, p = 9.85 x 10P -5 P, combined Stage 1 and in silico validation OR 1.13, p = 5.24 x 10P -6 P). This SNP has previously been reported to alter the expression of PTPLAD1, a gene involved in the synthesis of very long fatty acid chains and in the Rac1 signaling pathway. Meta-analysis of gene expression microarray data found PTPLAD1 to display increased expression in the aggressive non-endometrioid histology compared with endometrioid endometrial cancer, suggesting that the causal SNP underlying the observed genetic association may influence expression of this gene. Neither rs8035725 nor significant SNPs identified by imputation were predicted bioinformatically to affect transcription factor binding sites, indicating that further studies are required to assess their potential effect on other regulatory elements. The other grade- associated SNP, rs6606792, is located upstream of an inferred pseudogene, ELMO2P1 (Stage 1 OR 1.12, p = 5 x 10P -5 P; combined Stage 1 and in silico validation OR 1.09, p = 3.56 x 10P -5 P). Imputation of the ±1 Mb region surrounding this SNP revealed a cluster of significantly associated variants which are predicted to abolish various transcription factor binding sites, and would be expected to decrease gene expression. ELMO2P1 was not included on the microarray platforms collected for this PhD, and so its expression could not be investigated. However, the high sequence homology of ELMO2P1 with ELMO2, a gene important to cell motility, indicates that ELMO2 could be the parent gene for ELMO2P1 and as such, ELMO2P1 could function to regulate the expression of ELMO2. Increased expression of ELMO2 was seen to be associated with increasing endometrial cancer grade, as well as with aggressive endometrial cancer histological subtypes by microarray meta-analysis. Thus, it is hypothesised that SNPs in linkage disequilibrium with rs6606792 decrease the transcription of ELMO2P1, reducing the regulatory effect of ELMO2P1 on ELMO2 expression. Consequently, ELMO2 expression is increased, cell motility is enhanced leading to an aggressive endometrial cancer phenotype. In summary, these findings have identified several areas of research for further study. The results presented in this thesis provide evidence that a SNP in PGR is associated with risk of developing endometrial cancer. This PhD study also reports two independent loci on chromosome 15 to be associated with increased endometrial cancer grade, and furthermore, genes associated with these SNPs to be differentially expressed according in aggressive subtypes and/or by grade. The studies reported in this thesis support the need for comprehensive SNP association studies on prioritised MMP, TIMP and KLK genes in large sample sets. Until these studies are performed, the role of MMP, TIMP and KLK genetic variation remains unclear. Overall, this PhD study has contributed to the understanding of genetic variation involvement in endometrial cancer susceptibility and prognosis. Importantly, the genetic regions highlighted in this study could lead to the identification of novel gene targets to better understand the biology of endometrial cancer and also aid in the development of therapeutics directed at treating this disease.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.