Genome sequencing differentiates a paracentric inversion from a balanced insertion enabling more accurate preimplantation genetic testing.

INTRODUCTION: Distinguishing paracentric inversions (PAIs) from chromosomal insertions has traditionally relied on fluorescent in situ hybridization (FISH) techniques, but recent advancements in high-throughput sequencing have enabled the use of genome sequencing for such differentiation. In this study, we present a 38-year-old male carrier of a paracentric inversion on chromosome 2q, inv (2)(q31.2q34), whose partner experienced recurrent miscarriages. MATERIAL AND METHODS: FISH analysis confirmed the inversion, and genome sequencing was employed for detailed characterization. RESULTS: Preimplantation genetic testing (PGT) revealed that all assessed embryos were balanced, consistent with the low risk of unbalanced offspring associated with PAIs. While PAI carriers traditionally exhibit low risk of producing unbalanced offspring, exceptions exist due to crossover events within the inversion loop. Although the sample size was limited, the findings align with existing sperm study data, supporting the rare occurrence of unbalanced progeny in PAI carriers. CONCLUSIONS: This study highlights the possibility of characterizing PAIs using genome sequencing to enable correct reproductive counseling and PGT decisions. Detailed characterization of a PAI is crucial for understanding potential outcomes and guiding PGT strategies, as accurate knowledge of the inversion size is essential for appropriate method selection in PGT. Given the very low risk of unbalanced offspring in PAI carriers, routine PGT may not be warranted but should be considered in specific cases with a history of unbalanced progeny or recurrent miscarriages. This study contributes to our understanding of PAI segregation and its implications for reproductive outcomes.

Genome sequencing with comprehensive variant calling identifies structural variants and repeat expansions in a large fraction of individuals with ataxia and/or neuromuscular disorders.

Introduction: Neuromuscular disorders (NMDs) have a heterogeneous etiology. A genetic diagnosis is key to personalized healthcare and access to targeted treatment for the affected individuals. Methods: In this study, 861 patients with NMDs were analyzed with genome sequencing and comprehensive variant calling including single nucleotide variants, small insertions/deletions (SNVs/INDELs), and structural variants (SVs) in a panel of 895 NMD genes, as well as short tandem repeat expansions (STRs) at 28 loci. In addition, for unsolved cases with an unspecific clinical presentation, the analysis of a panel with OMIM disease genes was added. Results: In the cohort, 27% (232/861) of the patients harbored pathogenic variants, of which STRs and SVs accounted for one-third of the patients (71/232). The variants were found in 107 different NMD genes. Furthermore, 18 pediatric patients harbored pathogenic variants in non-NMD genes. Discussion: Our results highlight that for children with unspecific hypotonia, a genome-wide analysis rather than a disease-based gene panel should be considered as a diagnostic approach. More importantly, our results clearly show that it is crucial to include STR- and SV-analyses in the diagnostics of patients with neuromuscular disorders.

Genome sequencing is a sensitive first-line test to diagnose individuals with intellectual disability.

PURPOSE: Individuals with intellectual disability (ID) and/or neurodevelopment disorders (NDDs) are currently investigated with several different approaches in clinical genetic diagnostics. METHODS: We compared the results from 3 diagnostic pipelines in patients with ID/NDD: genome sequencing (GS) first (N = 100), GS as a secondary test (N = 129), or chromosomal microarray (CMA) with or without FMR1 analysis (N = 421). RESULTS: The diagnostic yield was 35% (GS-first), 26% (GS as a secondary test), and 11% (CMA/FMR1). Notably, the age of diagnosis was delayed by 1 year when GS was performed as a secondary test and the cost per diagnosed individual was 36% lower with GS first than with CMA/FMR1. Furthermore, 91% of those with a negative result after CMA/FMR1 analysis (338 individuals) have not yet been referred for additional genetic testing and remain undiagnosed. CONCLUSION: Our findings strongly suggest that genome analysis outperforms other testing strategies and should replace traditional CMA and FMR1 analysis as a first-line genetic test in individuals with ID/NDD. GS is a sensitive, time- and cost-effective method that results in a confirmed molecular diagnosis in 35% of all referred patients.

Unravelling genetic variants of a swedish family with high risk of prostate cancer.

BACKGROUND: Prostate cancer is the most prevalent cancer in men worldwide. It is a polygenic disease with a substantial proportion of heritability. Identification of novel candidate biomarkers is crucial for clinical cancer prevention and the development of therapeutic strategies. Here, we describe the analysis of rare and common genetic variants that can predispose to the development of prostate cancer. METHODS: Whole-genome sequencing was performed on germline DNA of five Swedish siblings which were diagnosed with prostate cancer. The high-risk variants were identified setting the minor allele frequency < 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.5, while the low-risk variants were identified minor allele frequency > 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.1. RESULTS: We identified 38 candidate high-risk gene variants and 332 candidate low-risk gene variants, where 2 and 14 variants were in coding regions, respectively, that were shared by the brothers with prostate cancer. CONCLUSIONS: This study expanded the knowledge of potential risk factor candidates involved in hereditary and familial prostate cancer. Our findings can be beneficial when applying targeted screening in families with a high risk of developing the disease.

Whole-Exome Sequencing of Germline Variants in Non-BRCA Families with Hereditary Breast Cancer.

Breast cancer is the most prevalent malignancy among women worldwide and hereditary breast cancer (HBC) accounts for about 5−10% of the cases. Today, the most recurrent genes known are BRCA1 and BRCA2, accounting for around 25% of familial cases. Although thousands of loss-of-function variants in more than twenty predisposing genes have been found, the majority of familial cases of HBC remain unexplained. The aim of this study was to identify new predisposing genes for HBC in three non-BRCA families with autosomal dominant inheritance pattern using whole-exome sequencing and functional prediction tools. No pathogenic variants in known hereditary cancer-related genes could explain the breast cancer susceptibility in these families. Among 2122 exonic variants with maximum minor allele frequency (MMAF) < 0.1%, between 17−35 variants with combined annotation-dependent depletion (CADD) > 20 segregated with disease in the three analyzed families. Selected candidate genes, i.e., UBASH3A, MYH13, UTP11L, and PAX7, were further evaluated using protein expression analysis but no alterations of cancer-related pathways were observed. In conclusion, identification of new high-risk cancer genes using whole-exome sequencing has been more challenging than initially anticipated, in spite of selected families with pronounced family history of breast cancer. A combination of low- and intermediate-genetic-risk variants may instead contribute the breast cancer susceptibility in these families.

A Swedish Genome-Wide Haplotype Association Analysis Identifies a Novel Breast Cancer Susceptibility Locus in 8p21.2 and Characterizes Three Loci on Chromosomes 10, 11 and 16.

(1) Background: The heritability of breast cancer is partly explained but much of the genetic contribution remains to be identified. Haplotypes are often used as markers of ethnicity as they are preserved through generations. We have previously demonstrated that haplotype analysis, in addition to standard SNP association studies, could give novel and more detailed information on genetic cancer susceptibility. (2) Methods: In order to examine the association of a SNP or a haplotype to breast cancer risk, we performed a genome wide haplotype association study, using sliding window analysis of window sizes 1−25 and 50 SNPs, in 3200 Swedish breast cancer cases and 5021 controls. (3) Results: We identified a novel breast cancer susceptibility locus in 8p21.1 (OR 2.08; p 3.92 × 10−8), confirmed three known loci in 10q26.13, 11q13.3, 16q12.1-2 and further identified novel subloci within these three loci. Altogether 76 risk SNPs, 3302 risk haplotypes of window size 2−25 and 113 risk haplotypes of window size 50 at p < 5 × 10−8 on chromosomes 8, 10, 11 and 16 were identified. In the known loci haplotype analysis reached an OR of 1.48 in overall breast cancer and in familial cases OR 1.68. (4) Conclusions: Analyzing haplotypes, rather than single variants, could detect novel susceptibility loci even in small study populations but the method requires a fairly homogenous study population.

Lower Intake of Saturated Fatty Acids Is Associated with Improved Lipid Profile in a 6-Year-Old Nationally Representative Population.

To strengthen the organization of new national dietary surveys and interventions in childhood, our aim was to study macronutrient intake and blood lipid profile at 6 years of age by comparing results from two earlier population-based cohorts. Subjects were n = 131 and n = 162 in the years 2001−2002 and 2011−2012, respectively. Three-day weighed food records were used to estimate diet and calculate nutrient intake. Total cholesterol, HDL-cholesterol and triacylglycerol were measured in serum and LDL-cholesterol was calculated. The average intake of saturated fatty acids (SFA) and trans FA was lower in 2011−2012 than 2001−2002 (13.3E% vs. 14.7E%, p < 0.001, and 0.8E% vs. 1.4E%, p < 0.001, respectively), replaced by a higher intake of unsaturated fatty acids. Total cholesterol and LDL-cholesterol were significantly lower in 2011−2012 than 2001−2002 (4.6 vs. 4.4 mmol/L, p = 0.003 and 2.8 vs. 2.5 mmol/L, p < 0.001, respectively). In a multiple linear regression model, one E% increase in SFA intake was related to a 0.03 mmol/L increase in LDL cholesterol (p = 0.04). A lower intake of saturated and trans fatty acids, replaced by unsaturated fatty acids, may have contributed to an improved lipid profile in a healthy 6-year-old population. Biological data for analysis of blood lipids are important in national dietary surveys in healthy children to monitor important health outcomes of interventions.

Colorectal cancer risk susceptibility loci in a Swedish population.

To search for colorectal cancer (CRC) risk loci, Swedish samples were used for a genome-wide haplotype analysis. A logistic regression model was employed in 2663 CRC cases and 1642 controls in the discovery analysis. Three analyses were done, on all, familial-, and nonfamilial CRC samples and only results with odds ratio (OR) > 1 were analyzed. single nucleotide polymorphism (SNP) analysis did not generate any statistically significant results. Haplotype analysis suggested novel loci, on chromosome 2q36.1 (OR = 1.71, p value = 5.6924 × 10-8 ) in all CRC samples, chromosome 1q43 (OR = 4.04 p value = 3.24 × 10-8 ) in familial CRC samples, and two hits in nonfamilial CRC samples, chromosomes 2q36.1 (OR = 1.71 p value = 5.69 × 10-8 ) and 3p24.3 (OR = 1.62 p value = 6.21 × 10-9 ). Moreover, one locus on chromosome 20q13.33 was suggested in analyses of all samples, and five more novel loci were suggested on chromosomes 10q25.3, 15q,22.31, 17p11.2, 1p34.2, and 3q24. The haplotypes from the analysis of all samples were replicated in a second study of CRC cases and controls from the same part of Sweden. In summary, using haplotype analysis in Swedish CRC samples, the best hits were novel loci and the locus on chromosomes 2q36.1 and 20q13.33 suggested in the analysis of all samples were confirmed in a second cohort. The ORs were often higher than ORs from published genome-wide association study (GWAS). The study suggested it was possible that a risk locus could involve more than one gene, and that haplotypes could give information on the gene or genes possibly involved in the risk at specific locus.

Sequencing for germline mutations in Swedish breast cancer families reveals novel breast cancer risk genes.

Identifying genetic cancer risk factors will lead to improved genetic counseling, cancer prevention and cancer care. Analyzing families with a strong history of breast cancer (BC) has been a successful method to identify genes that contribute to the disease. This has led to discoveries of high-risk genes like the BRCA-genes. Nevertheless, many BC incidences are of unknown causes. In this study, exome sequencing on 59 BC patients from 24 Swedish families with a strong history of BC was performed to identify variants in known and novel BC predisposing genes. First, we screened known BC genes and identified two pathogenic variants in the BRIP1 and PALB2 genes. Secondly, to identify novel BC genes, rare and high impact variants and segregating in families were analyzed to identify 544 variants in novel BC candidate genes. Of those, 22 variants were defined as high-risk variants. Several interesting genes, either previously linked with BC or in pathways that when flawed could contribute to BC, were among the detected genes. The strongest candidates identified are the FANCM gene, involved in DNA double-strand break repair, and the RAD54L gene, involved in DNA recombination. Our study shows identifying pathogenic variants is challenging despite a strong family history of BC. Several interesting candidates were observed here that need to be further studied.

Identification of known and novel familial cancer genes in Swedish colorectal cancer families.

Identifying new candidate colorectal cancer (CRC) genes and mutations are important for clinical cancer prevention as well as in cancer care. Genetic counseling is already implemented for known high-risk variants; however, the majority of CRC are of unknown causes. In our study, 110 CRC patients in 55 Swedish families with a strong history of CRC but unknown genetic causes were analyzed with the aim of identifying novel candidate CRC predisposing genes. Exome sequencing was used to identify rare and high-impact variants enriched in the families. No clear pathogenic variants were found in known CRC predisposing genes; however, potential pathogenic variants in novel CRC predisposing genes were identified. Over 3000 variants with minor allele frequency (MAF) <0.01 and Combined Annotation Dependent Depletion (CADD) > 20 were seen aggregating in the CRC families. Of those, 27 variants with MAF < 0.001 and CADD>25 were considered high-risk mutations. Interestingly, more than half of the high-risk variants were detected in three families, suggesting cumulating contribution of several variants to CRC. In summary, our study shows that despite a strong history of CRC within families, identifying pathogenic variants is challenging. In a small number of families, few rare mutations were shared by affected family members. This could indicate that in the absence of known CRC predisposing genes, a cumulating contribution of mutations leads to CRC observed in these families.

Whole genome DNA sequencing provides an atlas of somatic mutagenesis in healthy human cells and identifies a tumor-prone cell type.

BACKGROUND: The lifelong accumulation of somatic mutations underlies age-related phenotypes and cancer. Mutagenic forces are thought to shape the genome of aging cells in a tissue-specific way. Whole genome analyses of somatic mutation patterns, based on both types and genomic distribution of variants, can shed light on specific processes active in different human tissues and their effect on the transition to cancer. RESULTS: To analyze somatic mutation patterns, we compile a comprehensive genetic atlas of somatic mutations in healthy human cells. High-confidence variants are obtained from newly generated and publicly available whole genome DNA sequencing data from single non-cancer cells, clonally expanded in vitro. To enable a well-controlled comparison of different cell types, we obtain single genome data (92% mean coverage) from multi-organ biopsies from the same donors. These data show multiple cell types that are protected from mutagens and display a stereotyped mutation profile, despite their origin from different tissues. Conversely, the same tissue harbors cells with distinct mutation profiles associated to different differentiation states. Analyses of mutation rate in the coding and non-coding portions of the genome identify a cell type bearing a unique mutation pattern characterized by mutation enrichment in active chromatin, regulatory, and transcribed regions. CONCLUSIONS: Our analysis of normal cells from healthy donors identifies a somatic mutation landscape that enhances the risk of tumor transformation in a specific cell population from the kidney proximal tubule. This unique pattern is characterized by high rate of mutation accumulation during adult life and specific targeting of expressed genes and regulatory regions.

Somatic mutation that affects transcription factor binding upstream of CD55 in the temporal cortex of a late-onset Alzheimer disease patient.

Alzheimer's disease (AD) is the most common neurodegenerative disease worldwide. Familial cases suggest genetic components; however, monogenetic causes are few, and the vast majority of incidences have unknown cause. Sequencing efforts have focused on germline mutations, but improved technology has opened up for studies on somatic mutations in affected brain tissue samples. Here we use ultra-deep sequencing on brain and blood from early-onset AD (EOAD) and late-onset AD (LOAD) patients and non-AD individuals (n = 16). In total, 2.86 Mb of genomic regions, previously associated with AD, were targeted included 28 genes and upstream and downstream regulatory regions. Tailored downstream bioinformatics filtering identified 11 somatic single nucleotide variants in the temporal cortex in AD patients and none in the controls. One variant was validated to be present at 0.4% allele frequency in temporal cortex of a LOAD patient. This variant was predicted to affect transcription factor binding sites upstream of the CD55 gene, contributing to AD pathogenesis by affecting the complement system. Our results suggest that future studies targeting larger portions of the genome for somatic mutation analysis are important to obtain an increased understanding for the molecular basis of both EOAD and LOAD.

Accumulation of Progerin Affects the Symmetry of Cell Division and Is Associated with Impaired Wnt Signaling and the Mislocalization of Nuclear Envelope Proteins.

Hutchinson-Gilford progeria syndrome (HGPS) is the result of a defective form of the lamin A protein called progerin. While progerin is known to disrupt the properties of the nuclear lamina, the underlying mechanisms responsible for the pathophysiology of HGPS remain less clear. Previous studies in our laboratory have shown that progerin expression in murine epidermal basal cells results in impaired stratification and halted development of the skin. Stratification and differentiation of the epidermis is regulated by asymmetric stem cell division. Here, we show that expression of progerin impairs the ability of stem cells to maintain tissue homeostasis as a result of altered cell division. Quantification of basal skin cells showed an increase in symmetric cell division that correlated with progerin accumulation in HGPS mice. Investigation of the mechanisms underlying this phenomenon revealed a putative role of Wnt/ß-catenin signaling. Further analysis suggested an alteration in the nuclear translocation of ß-catenin involving the inner and outer nuclear membrane proteins, emerin and nesprin-2. Taken together, our results suggest a direct involvement of progerin in the transmission of Wnt signaling and normal stem cell division. These insights into the molecular mechanisms of progerin may help develop new treatment strategies for HGPS.

Somatic mutagenesis in satellite cells associates with human skeletal muscle aging.

Human aging is associated with a decline in skeletal muscle (SkM) function and a reduction in the number and activity of satellite cells (SCs), the resident stem cells. To study the connection between SC aging and muscle impairment, we analyze the whole genome of single SC clones of the leg muscle vastus lateralis from healthy individuals of different ages (21-78 years). We find an accumulation rate of 13 somatic mutations per genome per year, consistent with proliferation of SCs in the healthy adult muscle. SkM-expressed genes are protected from mutations, but aging results in an increase in mutations in exons and promoters, targeting genes involved in SC activity and muscle function. In agreement with SC mutations affecting the whole tissue, we detect a missense mutation in a SC propagating to the muscle. Our results suggest somatic mutagenesis in SCs as a driving force in the age-related decline of SkM function.

A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation.

AIMS: Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in man, causing substantial morbidity and mortality with a major worldwide public health impact. It is increasingly recognized as a highly heritable condition. This study aimed to determine genetic risk factors for early-onset AF. METHODS AND RESULTS: We sequenced the whole genomes of 8453 Icelanders and imputed genotypes of the 25.5 million sequence variants we discovered into 1799 Icelanders with early-onset AF (diagnosed before 60 years of age) and 337 453 controls. Each sequence variant was tested for association based on multiplicative and recessive inheritance models. We discovered a rare frameshift deletion in the myosin MYL4 gene (c.234delC) that associates with early-onset AF under a recessive mode of inheritance (allelic frequency = 0.58%). We found eight homozygous carriers of the mutation, all of whom had early-onset AF. Six of the homozygotes were diagnosed by the age of 30 and the remaining two in their 50s. Three of the homozygotes had received pacemaker implantations due to sick sinus syndrome, three had suffered an ischemic stroke, and one suffered sudden cardiac death. CONCLUSIONS: Through a population approach we found a loss of function mutation in the myosin gene MYL4 that, in the homozygous state, is completely penetrant for early-onset AF. The finding may provide novel mechanistic insight into the pathophysiology of this complex arrhythmia.

HLA class II sequence variants influence tuberculosis risk in populations of European ancestry.

Mycobacterium tuberculosis infections cause 9 million new tuberculosis cases and 1.5 million deaths annually. To identify variants conferring risk of tuberculosis, we tested 28.3 million variants identified through whole-genome sequencing of 2,636 Icelanders for association with tuberculosis (8,162 cases and 277,643 controls), pulmonary tuberculosis (PTB) and M. tuberculosis infection. We found association of three variants in the region harboring genes encoding the class II human leukocyte antigens (HLAs): rs557011[T] (minor allele frequency (MAF) = 40.2%), associated with M. tuberculosis infection (odds ratio (OR) = 1.14, P = 3.1 × 10(-13)) and PTB (OR = 1.25, P = 5.8 × 10(-12)), and rs9271378[G] (MAF = 32.5%), associated with PTB (OR = 0.78, P = 2.5 × 10(-12))--both located between HLA-DQA1 and HLA-DRB1--and a missense variant encoding p.Ala210Thr in HLA-DQA1 (MAF = 19.1%, rs9272785), associated with M. tuberculosis infection (P = 9.3 × 10(-9), OR = 1.14). We replicated association of these variants with PTB in samples of European ancestry from Russia and Croatia (P < 5.9 × 10(-4)). These findings show that the HLA class II region contributes to genetic risk of tuberculosis, possibly through reduced presentation of protective M. tuberculosis antigens to T cells.

Heterozygote carriers for CNVs in PARK2 are at increased risk of Parkinson's disease.

Together with point mutations, homozygous deletions or duplications in PARK2 are responsible for the majority of autosomal recessive juvenile Parkinsonism. It is debated, however, whether heterozygous carriers of these mutations are at increased risk of Parkinson's disease (PD). Our goal was to determine whether heterozygous carriers of copy number variants (CNVs) affecting exons of the PARK2 gene are at risk of PD that is greater than that of non-carriers. We searched for CNVs affecting exons of PARK2 in a sample of 105 749 genotyped Icelanders. In total, 989 carriers, including 24 diagnosed with PD, were identified. The heterozygous carriers were tested for association in a sample of 1415 PD patients and 40 474 controls ≥65 years of age. PD patients were more often heterozygous carriers of PARK2 CNVs than controls [odds ratio (OR) = 1.69, P = 0.03] and compound heterozygous PD patients for a CNV and a missense mutation were not found. Furthermore, we conducted a meta-analysis of studies reporting on case-control samples screened for heterozygous PARK2 CNVs. Ten studies were included in the final analysis, with 4538 cases and 4213 controls. The pooled OR and P-value for the published and Icelandic results showed significant association between PARK2 CNVs and risk of PD (OR = 2.11, P = 2.54 × 10(-6)). Our analysis shows that heterozygous carriers of CNVs affecting exons of PARK2 have greater risk of PD than non-carriers.

Large-scale whole-genome sequencing of the Icelandic population.

Here we describe the insights gained from sequencing the whole genomes of 2,636 Icelanders to a median depth of 20×. We found 20 million SNPs and 1.5 million insertions-deletions (indels). We describe the density and frequency spectra of sequence variants in relation to their functional annotation, gene position, pathway and conservation score. We demonstrate an excess of homozygosity and rare protein-coding variants in Iceland. We imputed these variants into 104,220 individuals down to a minor allele frequency of 0.1% and found a recessive frameshift mutation in MYL4 that causes early-onset atrial fibrillation, several mutations in ABCB4 that increase risk of liver diseases and an intronic variant in GNAS associating with increased thyroid-stimulating hormone levels when maternally inherited. These data provide a study design that can be used to determine how variation in the sequence of the human genome gives rise to human diversity.

Severe osteoarthritis of the hand associates with common variants within the ALDH1A2 gene and with rare variants at 1p31.

Osteoarthritis is the most common form of arthritis and is a major cause of pain and disability in the elderly. To search for sequence variants that confer risk of osteoarthritis of the hand, we carried out a genome-wide association study (GWAS) in subjects with severe hand osteoarthritis, using variants identified through the whole-genome sequencing of 2,230 Icelanders. We found two significantly associated loci in the Icelandic discovery set: at 15q22 (frequency of 50.7%, odds ratio (OR) = 1.51, P = 3.99 × 10(-10)) in the ALDH1A2 gene and at 1p31 (frequency of 0.02%, OR = 50.6, P = 9.8 × 10(-10)). Among the carriers of the variant at 1p31 is a family with several members in whom the risk allele segregates with osteoarthritis. The variants within the ALDH1A2 gene were confirmed in replication sets from The Netherlands and the UK, yielding an overall association of OR = 1.46 and P = 1.1 × 10(-11) (rs3204689).

Assessment of osteoarthritis candidate genes in a meta-analysis of nine genome-wide association studies.

OBJECTIVE: To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA. METHODS: A total of 199 candidate genes for association with OA were identified using Human Genome Epidemiology (HuGE) Navigator. All of their single-nucleotide polymorphisms (SNPs) with an allele frequency of >5% were assessed by fixed-effects meta-analysis of 9 genome-wide association studies (GWAS) that included 5,636 patients with knee OA and 16,972 control subjects and 4,349 patients with hip OA and 17,836 control subjects of European ancestry. An additional 5,921 individuals were genotyped for significantly associated SNPs in the meta-analysis. After correction for the number of independent tests, P values less than 1.58 × 10(-5) were considered significant. RESULTS: SNPs at only 2 of the 199 candidate genes (COL11A1 and VEGF) were associated with OA in the meta-analysis. Two SNPs in COL11A1 showed association with hip OA in the combined analysis: rs4907986 (P = 1.29 × 10(-5) , odds ratio [OR] 1.12, 95% confidence interval [95% CI] 1.06-1.17) and rs1241164 (P = 1.47 × 10(-5) , OR 0.82, 95% CI 0.74-0.89). The sex-stratified analysis also showed association of COL11A1 SNP rs4908291 in women (P = 1.29 × 10(-5) , OR 0.87, 95% CI 0.82-0.92); this SNP showed linkage disequilibrium with rs4907986. A single SNP of VEGF, rs833058, showed association with hip OA in men (P = 1.35 × 10(-5) , OR 0.85, 95% CI 0.79-0.91). After additional samples were genotyped, association at one of the COL11A1 signals was reinforced, whereas association at VEGF was slightly weakened. CONCLUSION: Two candidate genes, COL11A1 and VEGF, were significantly associated with OA in this focused meta-analysis. The remaining candidate genes were not associated.