Spatiotemporal dynamics and genome‐wide association analysis of desiccation tolerance in Drosophila melanogaster

Water availability is a major environmental challenge to a variety of terrestrial organisms. In insects, desiccation tolerance varies predictably over spatial and temporal scales and is an important physiological determinant of fitness in natural populations. Here, we examine the dynamics of desiccation tolerance in North American populations of Drosophila melanogaster using: (a) natural populations sampled across latitudes and seasons; (b) experimental evolution in field mesocosms over seasonal time; (c) genome‐wide associations to identify SNPs/genes associated with variation for desiccation tolerance; and (d) subsequent analysis of patterns of clinal/seasonal enrichment in existing pooled sequencing data of populations sampled in both North America and Australia. A cline in desiccation tolerance was observed, for which tolerance exhibited a positive association with latitude; tolerance also varied predictably with culture temperature, demonstrating a significant degree of thermal plasticity. Desiccation tolerance evolved rapidly in field mesocosms, although only males showed differences in desiccation tolerance between spring and autumn collections from natural populations. Water loss rates did not vary significantly among latitudinal or seasonal populations; however, changes in metabolic rates during prolonged exposure to dry conditions are consistent with increased tolerance in higher latitude populations. Genome‐wide associations in a panel of inbred lines identified twenty‐five SNPs in twenty‐one loci associated with sex‐averaged desiccation tolerance, but there is no robust signal of spatially varying selection on genes associated with desiccation tolerance. Together, our results suggest that desiccation tolerance is a complex and important fitness component that evolves rapidly and predictably in natural populations.     

Multiplex SNP genotyping in pooled DNA samples by a four-colour microarray system

We selected 125 candidate single nucleotide polymorphisms (SNPs) in genes belonging to the human type 1 interferon (IFN) gene family and the genes coding for proteins in the main type 1 IFN signalling pathway by screening databases and by in silico comparison of DNA sequences. Using quantitative analysis of pooled DNA samples by solid-phase mini-sequencing, we found that only 20% of the candidate SNPs were polymorphic in the Finnish and Swedish populations. To allow more effective validation of candidate SNPs, we developed a four-colour microarray-based mini-sequencing assay for multiplex, quantitative allele frequency determination in pooled DNA samples. We used cyclic mini-sequencing reactions with primers carrying 5′-tag sequences, followed by capture of the products on microarrays by hybridisation to complementary tag oligonucleotides. Standard curves prepared from mixtures of known amounts of SNP alleles demonstrate the applicability of the system to quantitative analysis, and showed that for about half of the tested SNPs the limit of detection for the minority allele was below 5%. The microarray-based genotyping system established here is universally applicable for genotyping and quantification of any SNP, and the validated system for SNPs in type 1 IFN-related genes should find many applications in genetic studies of this important immunoregulatory pathway.     

Allele frequency matching between SNPs reveals an excess of linkage disequilibrium in genic regions of the human genome

Significant interest has emerged in mapping genetic susceptibility for complex traits through whole-genome association studies. These studies rely on the extent of association, i.e., linkage disequilibrium (LD), between single nucleotide polymorphisms (SNPs) across the human genome. LD describes the nonrandom association between SNP pairs and can be used as a metric when designing maximally informative panels of SNPs for association studies in human populations. Using data from the 1.58 million SNPs genotyped by Perlegen, we explored the allele frequency dependence of the LD statistic r(2) both empirically and theoretically. We show that average r(2) values between SNPs unmatched for allele frequency are always limited to much less than 1 (theoretical approximately 0.46 to 0.57 for this dataset). Frequency matching of SNP pairs provides a more sensitive measure for assessing the average decay of LD and generates average r(2) values across nearly the entire informative range (from 0 to 0.89 through 0.95). Additionally, we analyzed the extent of perfect LD (r(2) = 1.0) using frequency-matched SNPs and found significant differences in the extent of LD in genic regions versus intergenic regions. The SNP pairs exhibiting perfect LD showed a significant bias for derived, nonancestral alleles, providing evidence for positive natural selection in the human genome.     

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis‐regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low‐frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis‐regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic‐acid pathway to local adaptation.     

Associations of ATR and CHEK1 Single Nucleotide Polymorphisms with Breast Cancer

DNA damage and replication checkpoints mediated by the ATR-CHEK1 pathway are key to the maintenance of genome stability, and both ATR and CHEK1 have been proposed as potential breast cancer susceptibility genes. Many novel variants recently identified by the large resequencing projects have not yet been thoroughly tested in genome-wide association studies for breast cancer susceptibility. We therefore used a tagging SNP (tagSNP) approach based on recent SNP data available from the 1000 genomes projects, to investigate the roles of ATR and CHEK1 in breast cancer risk and survival. ATR and CHEK1 tagSNPs were genotyped in the Sheffield Breast Cancer Study (SBCS; 1011 cases and 1024 controls) using Illumina GoldenGate assays. Untyped SNPs were imputed using IMPUTE2, and associations between genotype and breast cancer risk and survival were evaluated using logistic and Cox proportional hazard regression models respectively on a per allele basis. Significant associations were further examined in a meta-analysis of published data or confirmed in the Utah Breast Cancer Study (UBCS). The most significant associations for breast cancer risk in SBCS came from rs6805118 in ATR (p=7.6x10^-5) and rs2155388 in CHEK1 (p=3.1x10^-6), but neither remained significant after meta-analysis with other studies. However, meta-analysis of published data revealed a weak association between the ATR SNP rs1802904 (minor allele frequency is 12%) and breast cancer risk, with a summary odds ratio (confidence interval) of 0.90 (0.83-0.98) [p=0.0185] for the minor allele. Further replication of this SNP in larger studies is warranted since it is located in the target region of 2 microRNAs. No evidence of any survival effects of ATR or CHEK1 SNPs were identified. We conclude that common alleles of ATR and CHEK1 are not implicated in breast cancer risk or survival, but we cannot exclude effects of rare alleles and of common alleles with very small effect sizes.     

Polymorphisms in recent GWA identified asthma genes CA10, SGK493, and CTNNA3 are associated with disease severity and treatment response in childhood asthma

Recent genome-wide association studies (GWAs) have identified several new genetic risk factors for asthma; however, their influence on disease behavior and treatment response is still unclear. The aim of our study was the association analysis of the most significant single nucleotide polymorphisms (SNPs) recently reported by GWAs in different phenotypes of childhood asthma and analysis of correlation between these SNPs and clinical parameters. We have genotyped 288 children with asthma and 276 healthy controls. We provided here first replication of bivariate associations between CA10 (p?=?0.001) and SGK493 (p?=?0.011) with asthma. In addition, we have identified new correlation between SNPs in CA10, SGK493, and CTNNA3 with asthma behavior and glucocorticoid treatment response. Asthma patients who carried G allele in SNP rs967676 in gene CA10 were associated with more pronounced airway obstruction, higher bronchial hyper-reactivity, and increased inflammation. Higher bronchial hyper-reactivity was also associated with C allele in SNP rs1440095 in gene SGK493 but only in nonatopic asthmatics. In addition, we found that patients who carried at least one T allele in SNP rs1786929 in CTNNA3 (p?=?0.022) and atopic patients who carried at least one G allele in SNP rs967676 in gene CA10 (p?=?0.034) had higher increase in pulmonary function after glucocorticoid therapy. Our results suggest genetic heterogeneity between atopic and nonatopic asthma. We provided further evidence that treatment response in childhood asthma is genetically predisposed, and we report here two novel SNPs in genes CA10 and CTNNA3 as potential pharmacogenetic biomarkers that could be used in personalized treatment in childhood asthma.     

Detection of quantitative trait loci in Bos indicus and Bos taurus cattle using genome-wide association studies

Background

The apparent effect of a single nucleotide polymorphism (SNP) on phenotype depends on the linkage disequilibrium (LD) between the SNP and a quantitative trait locus (QTL). However, the phase of LD between a SNP and a QTL may differ between Bos indicus and Bos taurus because they diverged at least one hundred thousand years ago. Here, we test the hypothesis that the apparent effect of a SNP on a quantitative trait depends on whether the SNP allele is inherited from a Bos taurus or Bos indicus ancestor.

Methods

Phenotype data on one or more traits and SNP genotype data for 10 181 cattle from Bos taurus, Bos indicus and composite breeds were used. All animals had genotypes for 729 068 SNPs (real or imputed). Chromosome segments were classified as originating from B. indicus or B. taurus on the basis of the haplotype of SNP alleles they contained. Consequently, SNP alleles were classified according to their sub-species origin. Three models were used for the association study: (1) conventional GWAS (genome-wide association study), fitting a single SNP effect regardless of subspecies origin, (2) interaction GWAS, fitting an interaction between SNP and subspecies-origin, and (3) best variable GWAS, fitting the most significant combination of SNP and sub-species origin.

Results

Fitting an interaction between SNP and subspecies origin resulted in more significant SNPs (i.e. more power) than a conventional GWAS. Thus, the effect of a SNP depends on the subspecies that the allele originates from. Also, most QTL segregated in only one subspecies, suggesting that many mutations that affect the traits studied occurred after divergence of the subspecies or the mutation became fixed or was lost in one of the subspecies.

Conclusions

The results imply that GWAS and genomic selection could gain power by distinguishing SNP alleles based on their subspecies origin, and that only few QTL segregate in both B. indicus and B. taurus cattle. Thus, the QTL that segregate in current populations likely resulted from mutations that occurred in one of the subspecies and can have both positive and negative effects on the traits. There was no evidence that selection has increased the frequency of alleles that increase body weight.     

Analysis of polymorphisms in the insulin-like growth factor 1 receptor (IGF1R) gene from Japanese quail selected for body weight

Insulin-like growth factor 1 receptor ( IGF1R ) is essential for the signalling of growth. In this study, we performed single nucleotide polymorphism (SNP) detection in the Japanese quail IGF1R coding region and an association study between SNPs and body weight in two lines (SS and LL) selected for large and small body weight. Of 21 SNPs obtained, a SNP at position AB292766:c.2293G>A led to the replacement of a valine with an isoleucine (V765I). The two lines were fixed for alternate alleles, with allele encoding valine fixed in the LL line. A significant effect of the SNP genotype was found on 10-week body weight ( P  < 0.01) and on 4- to 10-week and 6- to 10-week average daily gain ( P  < 0.05) in the F₂ family obtained from lines LL and SS. In six populations maintained in Japan or France, the frequency of allele encoding valine was higher than the allele encoding isoleucine.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.     

Replication Study Confirms Link between TSPAN18 Mutation and Schizophrenia in Han Chinese

Schizophrenia (SCZ) is a severe psychiatric disorder associated with many different risk factors, both genetic and environmental. A recent genome-wide association study (GWAS) of Han Chinese identified three single-nucleotide polymorphisms (SNPs rs11038167, rs11038172, and rs835784) in the tetraspanins gene TSPAN18 as possible susceptibility loci for schizophrenia. Hoping to validate these findings, we conducted a case-control study of Han Chinese with 1093 schizophrenia cases and 1022 healthy controls. Using the LDR-PCR method to genotype polymorphisms in TSPAN18, we found no significant differences (P>0.05) between patients and controls in either the allele or genotype frequency of the SNPs rs11038167 and rs11038172. We did find, however, that the frequency of the ‘A’ allele of SNP rs835784 is significantly higher in patients than in controls. We further observed a significant association (OR  = 1.197, 95%CI  = 1.047–1.369) between risk for SCZ and this ‘A’ allele. These results confirm the significant association, in Han Chinese populations, of increased SCZ risk and the variant of the TSPAN18 gene containing the ‘A’ allele of SNP rs835784.     

SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries

High-density single-nucleotide polymorphism (SNP) arrays have revolutionized the ability of genome-wide association studies to detect genomic regions harboring sequence variants that affect complex traits. Extensive numbers of validated SNPs with known allele frequencies are essential to construct genotyping assays with broad utility. We describe an economical, efficient, single-step method for SNP discovery, validation and characterization that uses deep sequencing of reduced representation libraries (RRLs) from specified target populations. Using nearly 50 million sequences generated on an Illumina Genome Analyzer from DNA of 66 cattle representing three populations, we identified 62,042 putative SNPs and predicted their allele frequencies. Genotype data for these 66 individuals validated 92% of 23,357 selected genome-wide SNPs, with a genotypic and sequence allele frequency correlation of r = 0.67. This approach for simultaneous de novo discovery of high-quality SNPs and population characterization of allele frequencies may be applied to any species with at least a partially sequenced genome.     

Comparison and confirmation of SNP-bud burst associations in European beech populations in Germany

European beech (Fagus sylvatica L.) is one of the most important forest tree species in Central Europe. Climate change scenarios predict an increase in annual mean temperatures that may cause earlier bud burst in spring, potentially leading to an increased late frost-risk. Despite the ecologic and economic importance of beech, knowledge about the molecular basis of bud burst is still scarce in this species. Here, an association analysis was used to detect SNPs that are significantly associated with beech bud burst. A translocation experiment was established with progenies of six different beech populations from three widely separated regions in Northeast, Central, and Southwest Germany. In total, 600 individuals of the translocation experiment were genotyped using a set of 46 SNPs located in bud burst candidate genes. The association analysis revealed seven SNPs significantly associated with bud burst, each SNP explaining only a few percent of the observed phenotypic variation. Since the same SNP set was used in a previous association analysis with European beech, we were able to compare and confirm significant associations between SNPs and bud burst in distinct beech populations growing in different environments.     

Genome-wide association analysis reveals cryptic alleles as an important factor in heterosis for fatness in chicken F2 population

Genome-wide association studies have become possible in the chicken because of the recent availability of the complete genome sequence, a polymorphism map and high-density single nucleotide polymorphism (SNP) genotyping platforms. We used these tools to study the genetic basis of a very high level of heterosis that was previously observed for fatness in two F₂ populations established by crossing one outbred broiler (meat-type) sire with dams from two unrelated, highly inbred, light-bodied lines (Fayoumi and Leghorn). In each F₂ population, selective genotyping was carried out using phenotypically extreme males for abdominal fat percentage (AF) and about 3000 SNPs. Single-point association analysis of about 500 informative SNPs per cross showed significant association ( P  < 0.01) of 15 and 24 markers with AF in the Broiler × Fayoumi and Broiler × Leghorn crosses respectively. These SNPs were on 10 chromosomes (GGA1, 2, 3, 4, 7, 8, 10, 12, 15 and 27). Interestingly, of the 39 SNPs that were significantly associated with AF, there were about twice as many homozygous genotypes associated with higher AF that traced back to the inbred lines alleles, although the broiler line had on average higher AF. These SNPs are considered to be associated with QTL with cryptic alleles. This study reveals cryptic alleles as an important factor in heterosis for fatness observed in two chicken F₂ populations, and suggests epistasis as the common underlying mechanism for heterosis and cryptic allele expression. The results of this study also demonstrate the power of high marker-density SNP association studies in discovering QTL that were not detected by previous microsatellite-based genotyping studies.     

Association analysis of the RGS4 gene in Han Chinese and Scottish populations with schizophrenia

We investigated the RGS4 as a susceptibility gene for schizophrenia in Chinese Han (184 trios and 138 sibling pairs, a total of 322 families) and Scottish (580 cases and 620 controls) populations using both a family trio and case-control design. Both the samples had statistical power greater than 70% to detect a heterozygote genotype relative risk of >1.2 for frequent RGS4-risk alleles. We genotyped four single nucleotide polymorphisms (SNPs) which have previously been associated with schizophrenia as either individually or part of haplotypes. Allele frequencies and linkage disequilibrium between the SNPs was similar in the two populations. In the Chinese sample, no individual SNPs or any of their haplotypes were associated with schizophrenia. In the Scottish population, one SNP (SNP7) was significantly over-represented in the cases compared with the controls (0.44 vs. 0.38; A allele; chi(2) 7.08, P = 0.011 after correction for correlation between markers by permutation testing). One two-marker haplotype, composed of alleles T and A of SNP4 and SNP7, respectively, showed individual significance after correction by permutation testing (chi(2) 6.8; P = 0.04). None of the full four-marker haplotypes showed association, including the G-G-G-G haplotype previously associated with schizophrenia in more than one sample and the A-T-A-A haplotype. Thus, our data do not directly replicate previous associations of RGS4, but association with SNP 7 in the Scottish population provides some support for a role in schizophrenia susceptibility. We cannot conclusively exclude RGS4, as associated haplotypes are likely to be surrogates for unknown causative alleles, whose relationship with overlying haplotypes may differ between the population groups. Differences in the association seen across the two populations could result from methodological factors such as diagnostic differences but most likely result from ethnic differences in haplotype structures within RGS4.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.     

Population structure of chum salmon and selection on the markers collected for stock identification

Genetic stock identification (GSI) is a major management tool of Pacific salmon (Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single‐nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon (Oncorhynchus keta), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region, and NADH‐3 combined) in samples from 495 locations in the same distribution range (n = 61,813). TreeMix analysis of the microsatellite loci identified the greatest convergence toward Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The largest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was largest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. A regression of SNP population structures on those of microsatellites indicated the selection of the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites and performed regression analyses of SNP allele frequencies on a 2‐dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F _ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis showed differentiation of Japanese/Korean populations. The top five outlier SNPs included mtDNA3, U502241 (unknown), GnRH373, ras1362, and TCP178, which were identified by principal component analysis. We summarized the functions of 53 nuclear genes surrounding SNPs and the mtDNA3 locus by referring to a gene database system and propose how they may influence the fitness of chum salmon.     

Association of adiponectin promoter variants with traits and clusters of metabolic syndrome in Arabs: Family-based study

Plasma levels of adiponectin are decreased in type 2 diabetes, obesity and hypertension. Our aim was to use a family-based analysis to identify the genetic variants of the adiponectin (ADIPOQ) gene that are associated with obesity, insulin resistance, dyslipidemia and hypertension, among Arabs. We screened 328 Arabs in one large extended family for single nucleotide polymorphisms (SNPs) in the promoter region of the ADIPOQ gene. Two common SNPs were detected: rs17300539 and rs266729. Evidences of association between traits related to the metabolic syndrome and the SNPs were studied by implementing quantitative genetic association analysis. Results showed that SNP rs266729 was significantly associated with body weight (p-value = 0.001), waist circumference (p-value = 0.037), BMI (p-value = 0.015) and percentage of total body fat (p-value = 0.003). Up to 4.1% of heritability of obesity traits was explained by the rs266729 locus. Further cross-sectional analysis showed that carriers of the G allele had significantly higher values of waist circumference, BMI and percentage of total body fat (p-values 0.014, 0.004 and 0.032, respectively). No association was detected between SNP rs266729 and other clusters of metabolic syndrome or their traits except for HOMA-IR and fasting plasma insulin levels, p-values 0.035 and 0.004, respectively. In contrast, both measured genotype and cross-sectional analysis failed to detect an association between the SNP rs17300539 with traits and clusters of metabolic syndrome. In conclusion, we showed family-based evidence of association of SNP rs266729 at ADIPOQ gene with traits defining obesity in Arab population. This is important for future prediction and prevention of obesity in population where obesity is in an increasing trend.     

TLR7 single-nucleotide polymorphisms in the 3' untranslated region and intron 2 independently contribute to systemic lupus erythematosus in Japanese women: a case-control association study

Introduction

The Toll-like receptor 7 (TLR7) gene, encoded on human chromosome Xp22.3, is crucial for type I interferon production. A recent multicenter study in East Asian populations, comprising Chinese, Korean and Japanese participants, identified an association of a TLR7 single-nucleotide polymorphism (SNP) located in the 3' untranslated region (3' UTR), rs3853839, with systemic lupus erythematosus (SLE), especially in males, although some difference was observed among the tested populations. To test whether additional polymorphisms contribute to SLE in Japanese, we systematically analyzed the association of TLR7 with SLE in a Japanese female population.

Methods

A case-control association study was conducted on eight tag SNPs in the TLR7 region, including rs3853839, in 344 Japanese females with SLE and 274 healthy female controls.

Results

In addition to rs3853839, two SNPs in intron 2, rs179019 and rs179010, which were in moderate linkage disequilibrium with each other (r ² = 0.53), showed an association with SLE (rs179019: P = 0.016, odds ratio (OR) 2.02, 95% confidence interval (95% CI) 1.15 to 3.54; rs179010: P = 0.018, OR 1.75, 95% CI 1.10 to 2.80 (both under the recessive model)). Conditional logistic regression analysis revealed that the association of the intronic SNPs and the 3' UTR SNP remained significant after we adjusted them for each other. When only the patients and controls carrying the risk genotypes at the 3' UTR SNPpositionwere analyzed, the risk of SLE was significantly increased when the individuals also carried the risk genotypes at both of the intronic SNPs (P = 0.0043, OR 2.45, 95% CI 1.31 to 4.60). Furthermore, the haplotype containing the intronic risk alleles in addition to the 3' UTR risk allele was associated with SLE under the recessive model (P = 0.016, OR 2.37, 95% CI 1.17 to 4.80), but other haplotypes were not associated with SLE.

Conclusions

The TLR7 intronic SNPs rs179019 and rs179010 are associated with SLE independently of the 3' UTR SNP rs3853839 in Japanese women. Our findings support a role of TLR7 in predisposition for SLE in Asian populations.     

Association of polymorphisms in the Chr18q11.2 locus with tuberculosis in Chinese population

A GWAS study has reported that two single nucleotide polymorphisms (SNPs) were associated with predisposition to tuberculosis (TB) in African populations. These two loci represented the long-waited GWAS hits for TB susceptibility. To determine whether these two SNPs are associated with TB in Chinese population, we attempted an replication in a cohort of over one thousand Chinese TB patients and 1,280 healthy controls using melting temperature shift allele-specific genotyping analysis. We found that only SNP rs4331426 was significantly associated with TB in Chinese population (p = 0.011). However, the effect was opposite. The G allele of the SNP in Chinese population is a protective allele (OR = 0.62, 95 % CI 0.44–0.87), while it was the risk allele for African population (OR = 1.19, 95 % CI 1.12–1.26). No significance was found for SNP rs2335704. The results provided an independent support for a role in susceptibility to TB for SNP rs4331426. However, it also indicated that direct predisposition element to TB and the association effects may vary across ethnic groups.