Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)

Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype.

BACKGROUND: To date, no genome-wide association study (GWAS) has considered the combined phenotype of asthma with hay fever. Previous analyses of family data from the Tasmanian Longitudinal Health Study provide evidence that this phenotype has a stronger genetic cause than asthma without hay fever. OBJECTIVE: We sought to perform a GWAS of asthma with hay fever to identify variants associated with having both diseases. METHODS: We performed a meta-analysis of GWASs comparing persons with both physician-diagnosed asthma and hay fever (n = 6,685) with persons with neither disease (n = 14,091). RESULTS: At genome-wide significance, we identified 11 independent variants associated with the risk of having asthma with hay fever, including 2 associations reaching this level of significance with allergic disease for the first time: ZBTB10 (rs7009110; odds ratio [OR], 1.14; P = 4 × 10(-9)) and CLEC16A (rs62026376; OR, 1.17; P = 1 × 10(-8)). The rs62026376:C allele associated with increased asthma with hay fever risk has been found to be associated also with decreased expression of the nearby DEXI gene in monocytes. The 11 variants were associated with the risk of asthma and hay fever separately, but the estimated associations with the individual phenotypes were weaker than with the combined asthma with hay fever phenotype. A variant near LRRC32 was a stronger risk factor for hay fever than for asthma, whereas the reverse was observed for variants in/near GSDMA and TSLP. Single nucleotide polymorphisms with suggestive evidence for association with asthma with hay fever risk included rs41295115 near IL2RA (OR, 1.28; P = 5 × 10(-7)) and rs76043829 in TNS1 (OR, 1.23; P = 2 × 10(-6)). CONCLUSION: By focusing on the combined phenotype of asthma with hay fever, variants associated with the risk of allergic disease can be identified with greater efficiency.

Brief report: high-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population.

OBJECTIVE: Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. METHODS: A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. RESULTS: We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly(149) Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly(149) Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). CONCLUSION: This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk.

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies.

A novel serogenetic approach determines the community prevalence of celiac disease and informs improved diagnostic pathways.

BACKGROUND: Changing perspectives on the natural history of celiac disease (CD), new serology and genetic tests, and amended histological criteria for diagnosis cast doubt on past prevalence estimates for CD. We set out to establish a more accurate prevalence estimate for CD using a novel serogenetic approach. METHODS: The human leukocyte antigen (HLA)-DQ genotype was determined in 356 patients with 'biopsy-confirmed' CD, and in two age-stratified, randomly selected community cohorts of 1,390 women and 1,158 men. Sera were screened for CD-specific serology. RESULTS: Only five 'biopsy-confirmed' patients with CD did not possess the susceptibility alleles HLA-DQ2.5, DQ8, or DQ2.2, and four of these were misdiagnoses. HLA-DQ2.5, DQ8, or DQ2.2 was present in 56% of all women and men in the community cohorts. Transglutaminase (TG)-2 IgA and composite TG2/deamidated gliadin peptide (DGP) IgA/IgG were abnormal in 4.6% and 5.6%, respectively, of the community women and 6.9% and 6.9%, respectively, of the community men, but in the screen-positive group, only 71% and 75%, respectively, of women and 65% and 63%, respectively, of men possessed HLA-DQ2.5, DQ8, or DQ2.2. Medical review was possible for 41% of seropositive women and 50% of seropositive men, and led to biopsy-confirmed CD in 10 women (0.7%) and 6 men (0.5%), but based on relative risk for HLA-DQ2.5, DQ8, or DQ2.2 in all TG2 IgA or TG2/DGP IgA/IgG screen-positive subjects, CD affected 1.3% or 1.9%, respectively, of females and 1.3% or 1.2%, respectively, of men. Serogenetic data from these community cohorts indicated that testing screen positives for HLA-DQ, or carrying out HLA-DQ and further serology, could have reduced unnecessary gastroscopies due to false-positive serology by at least 40% and by over 70%, respectively. CONCLUSIONS: Screening with TG2 IgA serology and requiring biopsy confirmation caused the community prevalence of CD to be substantially underestimated. Testing for HLA-DQ genes and confirmatory serology could reduce the numbers of unnecessary gastroscopies.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease.

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6 × 10(-10), odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6 × 10(-4). OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6 × 10(-5), OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2 × 10(-5), OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2 × 10(-5), OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1 × 10(-4), OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9 × 10(-5), OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9 × 10(-4), OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis.

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10(-12)) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P < 10(-14)). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Identification of IL6R and chromosome 11q13.5 as risk loci for asthma.

BACKGROUND: We aimed to identify novel genetic variants affecting asthma risk, since these might provide novel insights into molecular mechanisms underlying the disease. METHODS: We did a genome-wide association study (GWAS) in 2669 physician-diagnosed asthmatics and 4528 controls from Australia. Seven loci were prioritised for replication after combining our results with those from the GABRIEL consortium (n=26,475), and these were tested in an additional 25,358 independent samples from four in-silico cohorts. Quantitative multi-marker scores of genetic load were constructed on the basis of results from the GABRIEL study and tested for association with asthma in our Australian GWAS dataset. FINDINGS: Two loci were confirmed to associate with asthma risk in the replication cohorts and reached genome-wide significance in the combined analysis of all available studies (n=57,800): rs4129267 (OR 1·09, combined p=2·4×10(-8)) in the interleukin-6 receptor (IL6R) gene and rs7130588 (OR 1·09, p=1·8×10(-8)) on chromosome 11q13.5 near the leucine-rich repeat containing 32 gene (LRRC32, also known as GARP). The 11q13.5 locus was significantly associated with atopic status among asthmatics (OR 1·33, p=7×10(-4)), suggesting that it is a risk factor for allergic but not non-allergic asthma. Multi-marker association results are consistent with a highly polygenic contribution to asthma risk, including loci with weak effects that might be shared with other immune-related diseases, such as NDFIP1, HLA-B, LPP, and BACH2. INTERPRETATION: The IL6R association further supports the hypothesis that cytokine signalling dysregulation affects asthma risk, and raises the possibility that an IL6R antagonist (tocilizumab) may be effective to treat the disease, perhaps in a genotype-dependent manner. Results for the 11q13.5 locus suggest that it directly increases the risk of allergic sensitisation which, in turn, increases the risk of subsequent development of asthma. Larger or more functionally focused studies are needed to characterise the many loci with modest effects that remain to be identified for asthma. FUNDING: National Health and Medical Research Council of Australia. A full list of funding sources is provided in the webappendix.

Glaucoma risk alleles at CDKN2B-AS1 are associated with lower intraocular pressure, normal-tension glaucoma, and advanced glaucoma.

PURPOSE: Genetic variation at the 9p21 locus encompassing the CDKN2B-AS1, CDKN2A, and CDKN2B genes has been associated with primary open-angle glaucoma (POAG) in several independent studies. This study aimed to dissect the association further and to determine genotype-phenotype correlations between genetic variation at this locus and a range of glaucoma-related traits in a large cohort of POAG patients. DESIGN: Comparative case series and case-control study. PARTICIPANTS: One thousand four hundred thirty-two POAG patients and 595 unaffected controls recruited from 2 population-based and 2 cross-sectional studies. METHODS: Each patient was genotyped at 9 single nucleotide polymorphisms (SNPs) previously associated with POAG at the 9p21 locus. Each SNP was assessed for association with each outcome measure using linear regression under an additive genetic model. Associated traits were explored further including adjustment for relevant covariates. Highest recorded intraocular pressure (IOP) also was analyzed both with and without correction for central corneal thickness (CCT) and was dichotomized into high-tension glaucoma and normal-tension glaucoma (NTG). MAIN OUTCOME MEASURES: Intraocular pressure and vertical cup-to-disc ratio (VCDR). RESULTS: Glaucoma risk alleles at 9p21, particularly, rs7049105 and rs10120688, were associated with the presence of both NTG and advanced POAG. The SNP rs10120688 was associated with greater VCDR after adjustment for covariates (P = 0.003; ß = 0.016; standard error, 0.006). In addition, multiple SNPs in the region were associated with reduced IOP, before and after adjustment for CCT. The SNP most significantly associated with IOP was also rs10120688 (P = 0.001; ß = -2.135; standard error, 0.634) after adjustment for covariates under an additive model. In a comparison of high-tension versus low-tension glaucoma, this SNP was also the most significantly associated, particularly when IOP was corrected for CCT before classification of the type of glaucoma (P = 0.0009; odds ratio, 0.63; 95% confidence interval, 0.48-0.83). CONCLUSIONS: Patients with POAG carrying the glaucoma risk alleles at the 9p21 locus have larger VCDR and lower IOP than POAG patients without these alleles. Carriers of these alleles seem to be predisposed to POAG developing at lower IOP levels and exhibit stronger associations with NTG and advanced glaucoma phenotypes. This may be of relevance when setting target pressures in patients carrying these risk alleles.

Genome-wide study of familial juvenile hyperuricaemic (gouty) nephropathy (FJHN) indicates a new locus, FJHN3, linked to chromosome 2p22.1-p21.

Familial juvenile hyperuricaemic (gouty) nephropathy (FJHN), is an autosomal dominant disease associated with a reduced fractional excretion of urate, and progressive renal failure. FJHN is genetically heterogeneous and due to mutations of three genes: uromodulin (UMOD), renin (REN) and hepatocyte nuclear factor-1beta (HNF-1ß) on chromosomes 16p12, 1q32.1, and 17q12, respectively. However, UMOD, REN or HNF-1ß mutations are found in only approximately 45% of FJHN probands, indicating the involvement of other genetic loci in approximately 55% of probands. To identify other FJHN loci, we performed a single nucleotide polymorphism (SNP)-based genome-wide linkage analysis, in six FJHN families in whom UMOD, HNF-1ß and REN mutations had been excluded. Parametric linkage analysis using a 'rare dominant' model established linkage in five of the six FJHN families, with a LOD score >+3, at 0% recombination, between FJHN and SNPs at chromosome 2p22.1-p21. Analysis of individual recombinants in two unrelated affected individuals defined a approximately 5.5 Mbp interval, flanked telomerically by SNP RS372139 and centromerically by RS896986 that contained the locus, designated FJHN3. The interval contains 28 genes, and DNA sequence analysis of the most likely candidate, solute carrier family 8 member 1 (SLC8A1), did not identify any abnormalities in the FJHN3 probands. FJHN3 is likely located within a approximately 5.5 Mbp interval on chromosome 2p22.1-p21, and identifying the genetic abnormality will help to further elucidate mechanisms predisposing to gout and renal failure.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population.

BACKGROUND: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. OBJECTIVE: To investigate these associations in the Han Chinese population. METHODS: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r(2)=0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. RESULTS: Robust association was detected for MMEL1 and CTLA4, and modest association was identified for another six loci: PADI4, STAT4, PRDM1, CDK6, TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6 × 10(-5) unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4, significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3 × 10(-5) unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. CONCLUSION: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese.

OBJECTIVES: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. METHODS: A case-control study was performed in a Han Chinese population of patients with AS (n = 775) and controls (n = 1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran-Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. RESULTS: SNPs in TNFRSF1A (rs4149577, p = 8.2 × 10⁻4), STAT3 (rs2293152, p = 0.0015; rs1053005, p = 0.017) and ERAP1 (rs27038, p = 0.0091; rs27037, p = 0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p = 0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p > 0.1). CONCLUSIONS: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

Variants in DNA double-strand break repair and DNA damage-response genes and susceptibility to lung and head and neck cancers.

Cigarette smoking is the major risk factor for lung cancer, and together with alcohol for head and neck (H--N) cancer. These genotoxics produced DNA damage and particularly double-strand breaks (DSB) that are removed by various repair pathways. To understand the initiation of these cancers, we performed a genotype analysis to correlate some variants in specific genes in a case-control study of lung and H-N cancers. In a discovery phase, we sequenced DNA samples of 32 healthy Caucasians to describe genetic variants in 30 genes involved in the repair of DSB and in DNA damage response. 625 variants were detected on 29 out of the 30 genes successfully screened by sequencing exons, parts of introns and flanking regions. These included 470 non-exonic variants, from which 33 insertions/deletions, and 155 exonic alterations, corresponding to 59 non synonymous polymorphisms. 223 of these variants were not previously described. In total, 379 variants were successfully genotyped in a case-control study restricted to smokers including 151 lung cases, 251 H-N cases, and 172 controls. To account for multiple testing, we associated to each p-value a proportion of false positives (q-value). Haplotype-analysis suggested potential associations (p < 0.05) between lung cancer and 2 genes (RECQL4 and RAD52), which came with q-value of 8%, and between H-N cancer and 1 gene (DNA-PK) but with q-value of 56%. The 3 genes are key players for regulating the efficiency of DSB repair. Large-scale studies are needed to show if any of these 3 variants are truly associated with an increased risk of cancer.

SNPs in the CpG island of NAP1L2: a possible link between DNA methylation and neural tube defects?

Deletion of the murine X-linked Nap1l2 gene causes lethality from midgestation onwards. The affected embryos exhibit neural tube defects (NTDs) closely resembling spina bifida and anencephaly in humans. X-linked familial and spontaneous cases of NTD were analyzed for sequence alterations in the human NAP1L2. No differences were found in the familial cases. However, a number of single nucleotide polymorphisms (SNPs) within the 5' region of NAP1L2 were identified both in cases of spontaneous NTD and in normal controls. Most of these SNPs lead to the replacement of guanidines or cytosines within a CpG island that is conserved between the human and the mouse promoter regions. Demethylation in vitro activates Nap1l2 transcriptional activity, suggesting the importance of the CpG island in regulating the activity of the Nap1l2/NAP1L2 genes, and the potential importance of the polymorphisms in modifying their transcriptional activity. NAP1L2/Nap1l2 expression may therefore depend on the genetic-environmental factors that are frequently associated with NTDs.

Novel risk loci for rheumatoid arthritis in Han Chinese and congruence with risk variants in Europeans.

OBJECTIVE: To investigate differences in genetic risk factors for rheumatoid arthritis (RA) in Han Chinese as compared with Europeans. METHODS: A genome-wide association study was conducted in China with 952 patients and 943 controls, and 32 variants were followed up in 2,132 patients and 2,553 controls. A transpopulation meta-analysis with results from a large European RA study was also performed to compare the genetic architecture across the 2 ethnic remote populations. RESULTS: Three non-major histocompatibility complex (non-MHC) loci were identified at the genome-wide significance level, the effect sizes of which were larger in anti-citrullinated protein antibody (ACPA)-positive patients than in ACPA-negative patients. These included 2 novel variants, rs12617656, located in an intron of DPP4 (odds ratio [OR] 1.56, P = 1.6 × 10(-21) ), and rs12379034, located in the coding region of CDK5RAP2 (OR 1.49, P = 1.1 × 10(-16) ), as well as a variant at the known CCR6 locus, rs1854853 (OR 0.71, P = 6.5 × 10(-15) ). The analysis of ACPA-positive patients versus ACPA-negative patients revealed that rs12617656 at the DPP4 locus showed a strong interaction effect with ACPAs (P = 5.3 × 10(-18) ), and such an interaction was also observed for rs7748270 at the MHC locus (P = 5.9 × 10(-8) ). The transpopulation meta-analysis showed genome-wide overlap and enrichment in association signals across the 2 populations, as confirmed by prediction analysis. CONCLUSION: This study has expanded the list of alleles that confer risk of RA, provided new insight into the pathogenesis of RA, and added empirical evidence to the emerging polygenic nature of complex trait variation driven by common genetic variants.

Comparison of the performance of two commercial genome-wide association study genotyping platforms in Han Chinese samples.

Most genome-wide association studies to date have been performed in populations of European descent, but there is increasing interest in expanding these studies to other populations. The performance of genotyping chips in Asian populations is not well established. Therefore, we sought to test the performance of widely used fixed-marker, genome-wide association studies chips in the Han Chinese population. Non-HapMap Chinese samples (n = 396) were genotyped using the Illumina OmniExpress and Affymetrix 6.0 platforms, whereas a subset also were genotyped using the Immunochip. Genotyped markers from the Affymetrix 6.0 and Illumina OmniExpress were used for full genome imputation based on the HapMap 2 JPT+CHB (Japanese from Tokyo, Japan and Chinese from Beijing, China) reference panel. The concordance between markers genotypes for the three platforms was very high whether directly genotyped or genotyped and imputed single nucleotide polymorphisms (SNPs; >99.8% for directly genotyped and >99.5% for genotyped and imputed SNPs, respectively) were compared. The OmniExpress chip data enabled more SNPs to be imputed, particularly SNPs with minor allele frequency >5%. The OmniExpress chip achieved better coverage of HapMap SNPs than the Affymetrix 6.0 chip (73.6% vs. 65.9%, respectively, for minor allele frequency >5%). The Affymetrix 6.0 and Illumina OmniExpress chip have similar genotyping accuracy and provide similar accuracy of imputed SNPs. The OmniExpress chip however provides better coverage of Asian HapMap SNPs, although its coverage of HapMap SNPs is moderate.

Association of genetic variants in the TMCO1 gene with clinical parameters related to glaucoma and characterization of the protein in the eye.

PURPOSE: Glaucoma is the leading cause of irreversible blindness worldwide. Primary open angle glaucoma (POAG) is the most common subtype. We recently reported association of genetic variants at chromosomal loci, 1q24 and 9p21, with POAG. In this study, we determined association of the most significantly associated single nucleotide polymorphism (SNP) rs4656461, at 1q24 near the TMCO1 gene, with the clinical parameters related to glaucoma risk and diagnosis, and determined ocular expression and subcellular localization of the human TMCO1 protein to understand the mechanism of its involvement in POAG. METHODS: Association of SNP rs4656461 with five clinical parameters was assessed in 1420 POAG cases using linear regression. The TMCO1 gene was screened for mutations in 95 cases with a strong family history and advanced disease. Ocular expression and subcellular localization of the TMCO1 protein were determined by immunolabeling and as GFP-fusion. RESULTS: The data suggest that individuals homozygous for the rs4656461 risk allele (GG) are 4 to 5 years younger at diagnosis than noncarriers of this allele. Our data demonstrate expression of the TMCO1 protein in most tissues in the human eye, including the trabecular meshwork and retina. However, the subcellular localization differs from that reported in other studies. We demonstrate that the endogenous protein localizes to the cytoplasm and nucleus in vivo and ex vivo. In the nucleus, the protein localizes to the nucleoli. CONCLUSIONS: This study shows a relationship between genetic variation in and around TMCO1 with age at diagnosis of POAG and provides clues to the potential cellular function/s of this gene.

Common variants at 12q14 and 12q24 are associated with hippocampal volume.

Aging is associated with reductions in hippocampal volume that are accelerated by Alzheimer's disease and vascular risk factors. Our genome-wide association study (GWAS) of dementia-free persons (n = 9,232) identified 46 SNPs at four loci with P values of <4.0 × 10(-7). In two additional samples (n = 2,318), associations were replicated at 12q14 within MSRB3-WIF1 (discovery and replication; rs17178006; P = 5.3 × 10(-11)) and at 12q24 near HRK-FBXW8 (rs7294919; P = 2.9 × 10(-11)). Remaining associations included one SNP at 2q24 within DPP4 (rs6741949; P = 2.9 × 10(-7)) and nine SNPs at 9p33 within ASTN2 (rs7852872; P = 1.0 × 10(-7)); along with the chromosome 12 associations, these loci were also associated with hippocampal volume (P < 0.05) in a third younger, more heterogeneous sample (n = 7,794). The SNP in ASTN2 also showed suggestive association with decline in cognition in a largely independent sample (n = 1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8) and neuronal migration (ASTN2), as well as enzymes targeted by new diabetes medications (DPP4), indicating new genetic influences on hippocampal size and possibly the risk of cognitive decline and dementia.

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1.

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10(-10)) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10(-9)). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10(-14), OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10(-14), OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma.