Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

Susceptibility to osteoporotic fracture is influenced by genetic factors that can be dissected by whole-genome linkage analysis in experimental animal crosses. The aim of this study was to characterize quantitative trait loci (QTLs) for biomechanical and two-dimensional dual-energy X-ray absorptiometry (DXA) phenotypes in reciprocal F2 crosses between diabetic GK and normo-glycemic F344 rat strains and to identify possible co-localization with previously reported QTLs for bone size and structure. The biomechanical measurements of rat tibia included ultimate force, stiffness and work to failure while DXA was used to characterize tibial area, bone mineral content (BMC) and areal bone mineral density (aBMD). F2 progeny (108 males, 98 females) were genotyped with 192 genome-wide markers followed by sex- and reciprocal cross-separated whole-genome QTL analyses. Significant QTLs were identified on chromosome 8 (tibial area; logarithm of odds (LOD) = 4.7 and BMC; LOD = 4.1) in males and on chromosome 1 (stiffness; LOD = 5.5) in females. No QTLs showed significant sex-specific interactions. In contrast, significant cross-specific interactions were identified on chromosome 2 (aBMD; LOD = 4.7) and chromosome 6 (BMC; LOD = 4.8) for males carrying F344mtDNA, and on chromosome 15 (ultimate force; LOD = 3.9) for males carrying GKmtDNA, confirming the effect of reciprocal cross on osteoporosis-related phenotypes. By combining identified QTLs for biomechanical-, size- and qualitative phenotypes (pQCT and 3D CT) from the same population, overlapping regions were detected on chromosomes 1, 3, 4, 6, 8 and 10. These are strong candidate regions in the search for genetic risk factors for osteoporosis.     

Genetic control of susceptibility to infection with Candida albicans in mice

Candida albicans is an opportunistic pathogen that causes acute disseminated infections in immunocompromised hosts, representing an important cause of morbidity and mortality in these patients. To study the genetic control of susceptibility to disseminated C. albicans in mice, we phenotyped a group of 23 phylogenetically distant inbred strains for susceptibility to infection as measured by extent of fungal replication in the kidney 48 hours following infection. Susceptibility was strongly associated with the loss-of-function mutant complement component 5 (C5/Hc) allele, which is known to be inherited by approximately 40% of inbred strains. Our survey identified 2 discordant strains, AKR/J (C5-deficient, resistant) and SM/J (C5-sufficient, susceptible), suggesting that additional genetic effects may control response to systemic candidiasis in these strains. Haplotype association mapping in the 23 strains using high density SNP maps revealed several putative loci regulating the extent of C. albicans replication, amongst which the most significant were C5 (P value = 2.43×10⁻¹¹) and a novel effect on distal chromosome 11 (P value = 7.63×10⁻⁹). Compared to other C5-deficient strains, infected AKR/J strain displays a reduced fungal burden in the brain, heart and kidney, and increased survival, concomitant with uniquely high levels of serum IFNγ. C5-independent genetic effects were further investigated by linkage analysis in an [A/JxAKR/J]F2 cross (n = 158) where the mutant Hc allele is fixed. These studies identified a chromosome 11 locus (Carg4, Candida albicans resistance gene 4; LOD = 4.59), and a chromosome 8 locus (Carg3; LOD = 3.95), both initially detected by haplotype association mapping. Alleles at both loci were inherited in a co-dominant manner. Our results verify the important effect of C5-deficiency in inbred mouse strains, and further identify two novel loci, Carg3 and Carg4, which regulate resistance to C. albicans infection in a C5-independent manner.     

Common variants in a novel gene, FONG on chromosome 2q33.1 confer risk of osteoporosis in Japanese

Osteoporosis is a common disease characterized by low bone mass, decreased bone quality and increased predisposition to fracture. Genetic factors have been implicated in its etiology; however, the specific genes related to susceptibility to osteoporosis are not entirely known. To detect susceptibility genes for osteoporosis, we conducted a genome-wide association study in Japanese using ∼270,000 SNPs in 1,747 subjects (190 cases and 1,557 controls) followed by multiple levels of replication of the association using a total of ∼5,000 subjects (2,092 cases and 3,114 controls). Through these staged association studies followed by resequencing and linkage disequilibrium mapping, we identified a single nucleotide polymorphism (SNP), rs7605378 associated with osteoporosis. (combined P = 1.51×10⁻⁸, odds ratio = 1.25). This SNP is in a previously unknown gene on chromosome 2q33.1, FONG. FONG is predicted to encode a 147 amino-acid protein with a formiminotransferase domain in its N-terminal (FTCD_N domain) and is ubiquitously expressed in various tissues including bone. Our findings would give a new insight into osteoporosis etiology and pathogenesis.     

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.     

Association of mitochondrial DNA variations with lung cancer risk in a Han Chinese population from southwestern China

Mitochondrial DNA (mtDNA) is particularly susceptible to oxidative damage and mutation due to the high rate of reactive oxygen species (ROS) production and limited DNA-repair capacity in mitochondrial. Previous studies demonstrated that the increased mtDNA copy number for compensation for damage, which was associated with cigarette smoking, has been found to be associated with lung cancer risk among heavy smokers. Given that the common and “non-pathological” mtDNA variations determine differences in oxidative phosphorylation performance and ROS production, an important determinant of lung cancer risk, we hypothesize that the mtDNA variations may play roles in lung cancer risk. To test this hypothesis, we conducted a case-control study to compare the frequencies of mtDNA haplogroups and an 822 bp mtDNA deletion between 422 lung cancer patients and 504 controls. Multivariate logistic regression analysis revealed that haplogroups D and F were related to individual lung cancer resistance (OR = 0.465, 95%CI = 0.329–0.656, p<0.001; and OR = 0.622, 95%CI = 0.425–0.909, p = 0.014, respectively), while haplogroups G and M7 might be risk factors for lung cancer (OR = 3.924, 95%CI = 1.757–6.689, p<0.001; and OR = 2.037, 95%CI = 1.253–3.312, p = 0.004, respectively). Additionally, multivariate logistic regression analysis revealed that cigarette smoking was a risk factor for the 822 bp mtDNA deletion. Furthermore, the increased frequencies of the mtDNA deletion in male cigarette smoking subjects of combined cases and controls with haplogroup D indicated that the haplogroup D might be susceptible to DNA damage from external ROS caused by heavy cigarette smoking.     

Mitochondrial DNA Copy Number and Risk of Oral Cancer: A Report from Northeast India

Background

Oral squamous cell carcinoma (OSCC) is the sixth most common cancer globally. Tobacco consumption and HPV infection, both are the major risk factor for the development of oral cancer and causes mitochondrial dysfunction. Genetic polymorphisms in xenobiotic-metabolizing enzymes modify the effect of environmental exposures, thereby playing a significant role in gene–environment interactions and hence contributing to the individual susceptibility to cancer. Here, we have investigated the association of tobacco - betel quid chewing, HPV infection, GSTM1-GSTT1 null genotypes, and tumour stages with mitochondrial DNA (mtDNA) content variation in oral cancer patients.

Methodology/Principal Findings

The study comprised of 124 cases of OSCC and 140 control subjects to PCR based detection was done for high-risk HPV using a consensus primer and multiplex PCR was done for detection of GSTM1-GSTT1 polymorphism. A comparative ΔCt method was used for determination of mtDNA content. The risk of OSCC increased with the ceased mtDNA copy number (P_trend = 0.003). The association between mtDNA copy number and OSCC risk was evident among tobacco – betel quid chewers rather than tobacco – betel quid non chewers; the interaction between mtDNA copy number and tobacco – betel quid was significant (P = 0.0005). Significant difference was observed between GSTM1 - GSTT1 null genotypes (P = 0.04, P = 0.001 respectively) and HPV infection (P<0.001) with mtDNA content variation in cases and controls. Positive correlation was found with decrease in mtDNA content with the increase in tumour stages (P<0.001). We are reporting for the first time the association of HPV infection and GSTM1-GSTT1 null genotypes with mtDNA content in OSCC.

Conclusion

Our results indicate that the mtDNA content in tumour tissues changes with tumour stage and tobacco-betel quid chewing habits while low levels of mtDNA content suggests invasive thereby serving as a biomarker in detection of OSCC.     

The genetics of reading disability in an often excluded sample: novel loci suggested for reading disability in rolandic epilepsy

Background

Reading disability (RD) is a common neurodevelopmental disorder with genetic basis established in families segregating “pure” dyslexia. RD commonly occurs in neurodevelopmental disorders including Rolandic Epilepsy (RE), a complex genetic disorder. We performed genomewide linkage analysis of RD in RE families, testing the hypotheses that RD in RE families is genetically heterogenenous to pure dyslexia, and shares genetic influences with other sub-phenotypes of RE.

Methods

We initially performed genome-wide linkage analysis using 1000 STR markers in 38 US families ascertained through a RE proband; most of these families were multiplex for RD. We analyzed the data by two-point and multipoint parametric LOD score methods. We then confirmed the linkage evidence in a second US dataset of 20 RE families. We also resequenced the SEMA3C gene at the 7q21 linkage locus in members of one multiplex RE/RD pedigree and the DISC1 gene in affected pedigrees at the 1q42 locus.

Results

In the discovery dataset there was suggestive evidence of linkage for RD to chromosome 7q21 (two-point LOD score 3.05, multipoint LOD 3.08) and at 1q42 (two-point LOD 2.87, multipoint LOD 3.03). Much of the linkage evidence at 7q21 derived from families of French-Canadian origin, whereas the linkage evidence at 1q42 was well distributed across all the families. There was little evidence for linkage at known dyslexia loci. Combining the discovery and confirmation datasets increased the evidence at 1q42 (two-point LOD = 3.49, multipoint HLOD = 4.70), but decreased evidence at 7q21 (two-point LOD = 2.28, multipoint HLOD  = 1.81), possibly because the replication sample did not have French Canadian representation.

Discussion

Reading disability in rolandic epilepsy has a genetic basis and may be influenced by loci at 1q42 and, in some populations, at 7q21; there is little evidence of a role for known DYX loci discovered in “pure” dyslexia pedigrees. 1q42 and 7q21 are candidate novel dyslexia loci.     

mtDNA nt13708A variant increases the risk of multiple sclerosis

Background

Mitochondrial DNA (mtDNA) polymorphism is a possible factor contributing to the maternal parent-of-origin effect in multiple sclerosis (MS) susceptibility.

Methods and Findings

In order to investigate the role of mtDNA variations in MS, we investigated six European MS case-control cohorts comprising >5,000 individuals. Three well matched cohorts were genotyped with seven common, potentially functional mtDNA single nucleotide polymorphisms (SNPs). A SNP, nt13708 G/A, was significantly associated with MS susceptibility in all three cohorts. The nt13708A allele was associated with an increased risk of MS (OR = 1.71, 95% CI 1.28–2.26, P = 0.0002). Subsequent sequencing of the mtDNA of 50 individuals revealed that the nt13708 itself, rather than SNPs linked to it, was responsible for the association. However, the association of nt13708 G/A with MS was not significant in MS cohorts which were not well case-control matched, indicating that the significance of association was affected by the population structure of controls.

Conclusions

Taken together, our finding identified the nt13708A variant as a susceptibility allele to MS, which could contribute to defining the role of the mitochondrial genome in MS pathogenesis.     

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.     

Quantitative Trait Loci (QTL) Associated with Resistance to a Monogenean Parasite (Benedenia seriolae) in Yellowtail (Seriola quinqueradiata) through Genome Wide Analysis

Benedenia infections caused by the monogenean fluke ectoparasite Benedenia seriolae seriously impact marine finfish aquaculture. Genetic variation has been inferred to play a significant role in determining the susceptibility to this parasitic disease. To evaluate the genetic basis of Benedenia disease resistance in yellowtail (Seriola quinqueradiata), a genome-wide and chromosome-wide linkage analyses were initiated using F₁ yellowtail families (n = 90 per family) based on a high-density linkage map with 860 microsatellite and 142 single nucleotide polymorphism (SNP) markers. Two major quantitative trait loci (QTL) regions on linkage groups Squ2 (BDR-1) and Squ20 (BDR-2) were identified. These QTL regions explained 32.9–35.5% of the phenotypic variance. On the other hand, we investigated the relationship between QTL for susceptibility to B. seriolae and QTL for fish body size. The QTL related to growth was found on another linkage group (Squ7). As a result, this is the first genetic evidence that contributes to detailing phenotypic resistance to Benedenia disease, and the results will help resolve the mechanism of resistance to this important parasitic infection of yellowtail.     

Mapping of a Chromosome 12 Region Associated with Airway Hyperresponsiveness in a Recombinant Congenic Mouse Strain and Selection of Potential Candidate Genes by Expression and Sequence Variation Analyses

In a previous study we determined that BcA86 mice, a strain belonging to a panel of AcB/BcA recombinant congenic strains, have an airway responsiveness phenotype resembling mice from the airway hyperresponsive A/J strain. The majority of the BcA86 genome is however from the hyporesponsive C57BL/6J strain. The aim of this study was to identify candidate regions and genes associated with airway hyperresponsiveness (AHR) by quantitative trait locus (QTL) analysis using the BcA86 strain. Airway responsiveness of 205 F2 mice generated from backcrossing BcA86 strain to C57BL/6J strain was measured and used for QTL analysis to identify genomic regions in linkage with AHR. Consomic mice for the QTL containing chromosomes were phenotyped to study the contribution of each chromosome to lung responsiveness. Candidate genes within the QTL were selected based on expression differences in mRNA from whole lungs, and the presence of coding non-synonymous mutations that were predicted to have a functional effect by amino acid substitution prediction tools. One QTL for AHR was identified on Chromosome 12 with its 95% confidence interval ranging from 54.6 to 82.6 Mbp and a maximum LOD score of 5.11 (p = 3.68×10⁻³). We confirmed that the genotype of mouse Chromosome 12 is an important determinant of lung responsiveness using a Chromosome 12 substitution strain. Mice with an A/J Chromosome 12 on a C57BL/6J background have an AHR phenotype similar to hyperresponsive strains A/J and BcA86. Within the QTL, genes with deleterious coding variants, such as Foxa1, and genes with expression differences, such as Mettl21d and Snapc1, were selected as possible candidates for the AHR phenotype. Overall, through QTL analysis of a recombinant congenic strain, microarray analysis and coding variant analysis we identified Chromosome 12 and three potential candidate genes to be in linkage with airway responsiveness.     

A case-control study of peripheral blood mitochondrial DNA copy number and risk of renal cell carcinoma

Background

Low mitochondrial DNA (mtDNA) copy number is a common feature of renal cell carcinoma (RCC), and may influence tumor development. Results from a recent case-control study suggest that low mtDNA copy number in peripheral blood may be a marker for increased RCC risk. In an attempt to replicate that finding, we measured mtDNA copy number in peripheral blood DNA from a U.S. population-based case-control study of RCC.

Methodology/Principal Findings

Relative mtDNA copy number was measured in triplicate by a quantitative real-time PCR assay using DNA extracted from peripheral whole blood. Cases (n = 603) had significantly lower mtDNA copy number than controls (n = 603; medians 0.85, 0.91 respectively; P = 0.0001). In multiple logistic regression analyses, the lowest quartile of mtDNA copy number was associated with a 60% increase in RCC risk relative to the highest quartile (OR = 1.6, 95% CI = 1.1–2.2; P _trend = 0.009). This association remained in analyses restricted to cases treated by surgery alone (OR _Q1 = 1.4, 95% CI = 1.0–2.1) and to localized tumors (2.0, 1.3–2.8).

Conclusions/Significance

Our findings from this investigation, to our knowledge the largest of its kind, offer important confirmatory evidence that low mtDNA copy number is associated with increased RCC risk. Additional research is needed to assess whether the association is replicable in prospective studies.     

Mitochondrial DNA haplogroup background affects LHON, but not suspected LHON, in Chinese patients

Recent studies have shown that mtDNA background could affect the clinical expression of Leber hereditary optic neuropathy (LHON). We analyzed the mitochondrial DNA (mtDNA) variation of 304 Chinese patients with m.11778G>A (sample #1) and of 843 suspected LHON patients who lack the three primary mutations (sample #2) to discern mtDNA haplogroup effect on disease onset. Haplogroup frequencies in the patient group was compared to frequencies in the general Han Chinese population (n = 1,689; sample #3). The overall matrilineal composition of the suspected LHON population resembles that of the general Han Chinese population, suggesting no association with mtDNA haplogroup. In contrast, analysis of these LHON patients confirms mtDNA haplogroup effect on LHON. Specifically, the LHON sample significantly differs from the general Han Chinese and suspected LHON populations by harboring an extremely lower frequency of haplogroup R9, in particular of its main sub-haplogroup F (#1 vs. #3, P-value = 1.46×10⁻¹⁷, OR = 0.051, 95% CI: 0.016–0.162; #1 vs. #2, P-value = 4.44×10⁻¹⁷, OR = 0.049, 95% CI: 0.015–0.154; in both cases, adjusted P-value <10⁻⁵) and higher frequencies of M7b (#1 vs. #3, adjusted P-value = 0.001 and #1 vs. #2, adjusted P-value = 0.004). Our result shows that mtDNA background affects LHON in Chinese patients with m.11778G>A but not suspected LHON. Haplogroup F has a protective effect against LHON, while M7b is a risk factor.     

Sex-specific genetic architecture of human fatness in Chinese: the SAPPHIRe Study

To dissect the genetic architecture of sexual dimorphism in obesity-related traits, we evaluated the sex–genotype interaction, sex-specific heritability and genome-wide linkages for seven measurements related to obesity. A total of 1,365 non-diabetic Chinese subjects from the family study of the Stanford Asia–Pacific Program of Hypertension and Insulin Resistance were used to search for quantitative trait loci (QTLs) responsible for the obesity-related traits. Pleiotropy and co-incidence effects from the QTLs were also examined using the bivariate linkage approach. We found that sex-specific differences in heritability and the genotype–sex interaction effects were substantially significant for most of these traits. Several QTLs with strong linkage evidence were identified after incorporating genotype by sex (G × S) interactions into the linkage mapping, including one QTL for hip circumference [maximum LOD score (MLS) = 4.22, empirical p = 0.000033] and two QTLs: for BMI on chromosome 12q with MLS 3.37 (empirical p = 0.0043) and 3.10 (empirical p = 0.0054). Sex-specific analyses demonstrated that these linkage signals all resulted from females rather than males. Most of these QTLs for obesity-related traits replicated the findings in other ethnic groups. Bivariate linkage analyses showed several obesity traits were influenced by a common set of QTLs. All regions with linkage signals were observed in one gender, but not in the whole sample, suggesting the genetic architecture of obesity-related traits does differ by gender. These findings are useful for further identification of the liability genes for these phenotypes through candidate genes or genome-wide association analysis.     

T Cell Activation Markers and African Mitochondrial DNA Haplogroups among Non-Hispanic Black Participants in AIDS Clinical Trials Group Study 384

Introduction

Mitochondrial function influences T cell dynamics and is affected by mitochondrial DNA (mtDNA) variation. We previously reported an association between African mtDNA haplogroup L2 and less robust CD4 cell recovery on antiretroviral therapy (ART) in non-Hispanic black ACTG 384 subjects. We explored whether additional T cell parameters in this cohort differed by mtDNA haplogroup.

Methods

ACTG 384 randomized ART-naïve subjects to two different nucleoside regimens with efavirenz, nelfinavir, or both. CD4 and CD8 memory and activation markers were available at baseline and week 48 on most subjects. mtDNA sequencing was performed on whole blood DNA, and haplogroups were determined. We studied non-Hispanic black subjects with HIV RNA <400 copies/mL at week 48. Analyses included Wilcoxon ranksum test and linear regression.

Results

Data from 104 subjects were included. Major African mtDNA haplogroups included L1 (N = 25), L2 (N = 31), and L3 (N = 32). Baseline age, HIV RNA, and CD4 cells did not differ between L2 and non-L2 haplogroups. Compared to non-L2 haplogroups, L2 subjects had lower baseline activated CD4 cells (median 12% vs. 17%; p = 0.03) and tended toward lower activated CD8 cells (41% vs. 47%; p = 0.06). At 48 weeks of ART, L2 subjects had smaller decreases in activated CD4 cells (−4% vs. −11%; p = 0.01), and smaller CD4 cell increases (+95 vs. +178; p = 0.002). In models adjusting for baseline age, CD4 cells, HIV RNA, and naïve-to-memory CD4 cell ratio, haplogroup L2 was associated with lower baseline (p = 0.04) and 48-week change in (p = 0.01) activated CD4 cells.

Conclusions

Among ART-naïve non-Hispanic blacks, mtDNA haplogroup L2 was associated with baseline and 48-week change in T cell activation, and poorer CD4 cell recovery. These data suggest mtDNA variation may influence CD4 T cell dynamics by modulating T cell activation. Further study is needed to replicate these associations and identify mechanisms.     

Mitochondrial DNA Haplogroup Confers Genetic Susceptibility to Nasopharyngeal Carcinoma in Chaoshanese from Guangdong,China

Recent studies have shown association of mtDNA background with cancer development. We analyzed mitochondrial DNA (mtDNA) control region variation of 201 patients with nasopharyngeal carcinoma (NPC) and of 201 normal controls from Chaoshan Han Chinese to discern mtDNA haplogroup effect on the disease onset. Binary logistic regression analysis with adjustment for gender and age revealed that the haplogroup R9 (P = 0.011, OR = 1.91, 95% CI = 1.16–3.16), particularly its sub-haplogroup F1 (P = 0.015, OR = 2.43, 95% CI = 1.18–5.00), were associated significantly with increased NPC risk. These haplogroups were further confirmed to confer high NPC risk in males and/or individuals ≥40 years of age, but not in females or in subjects <40 years old. Our results indicated that mtDNA background confers genetic susceptibility to NPC in Chaoshan Han Chinese, and R9, particularly its sub-haplogroup F1, is a risk factor for NPC.     

Use of Allele-Specific FAIRE to Determine Functional Regulatory Polymorphism Using Large-Scale Genotyping Arrays

Following the widespread use of genome-wide association studies (GWAS), focus is turning towards identification of causal variants rather than simply genetic markers of diseases and traits. As a step towards a high-throughput method to identify genome-wide, non-coding, functional regulatory variants, we describe the technique of allele-specific FAIRE, utilising large-scale genotyping technology (FAIRE-gen) to determine allelic effects on chromatin accessibility and regulatory potential. FAIRE-gen was explored using lymphoblastoid cells and the 50,000 SNP Illumina CVD BeadChip. The technique identified an allele-specific regulatory polymorphism within NR1H3 (coding for LXR-α), rs7120118, coinciding with a previously GWAS-identified SNP for HDL-C levels. This finding was confirmed using FAIRE-gen with the 200,000 SNP Illumina Metabochip and verified with the established method of TaqMan allelic discrimination. Examination of this SNP in two prospective Caucasian cohorts comprising 15,000 individuals confirmed the association with HDL-C levels (combined beta = 0.016; p = 0.0006), and analysis of gene expression identified an allelic association with LXR-α expression in heart tissue. Using increasingly comprehensive genotyping chips and distinct tissues for examination, FAIRE-gen has the potential to aid the identification of many causal SNPs associated with disease from GWAS.     

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD_C) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD_C, as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD_C associations that had p<1×10⁻⁵ in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD_C in all cohorts (overall p = 2×10⁻¹⁴, n = 5739). Each minor allele was associated with a decrease in BMD_C of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males −6.77mg/cm³ per C allele, p = 2×10⁻⁶; females −2.79 mg/cm³ per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.     

Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.