Targeted sequencing of genome wide significant loci associated with bone mineral density (BMD) reveals significant novel and rare variants: the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) targeted sequencing study.

Background: Bone mineral density (BMD) is a heritable phenotype that predicts fracture risk. We performed fine-mapping by targeted sequencing at WLS, MEF2C, ARHGAP1/F2 and JAG1 loci prioritized by eQTL and bioinformatic approaches among 56 BMD loci from our previous GWAS meta-analysis. Targeted sequencing was conducted in 1,291 Caucasians from the Framingham Heart Study ( n = 925) and Cardiovascular Health Study ( n = 366), including 206 women and men with extreme low femoral neck (FN) BMD. A total of 4,964 sequence variants (SNVs) were observed and 80% were rare with MAF <1%. The associations between previously identified SNPs in these loci and BMD, while nominally significant in sequenced participants, were no longer significant after multiple testing corrections. Conditional analyses did not find protein-coding variants that may be responsible for GWAS signals. On the other hand, in the sequenced subjects, we identified novel associations in WLS , ARHGAP1 , and 5' of MEF2C ( P- values < 8x10 - 5 ; false discovery rate (FDR) q-values < 0.01) that were much more strongly associated with BMD compared to the GWAS SNPs. These associated SNVs are less-common; independent from previous GWAS signals in the same loci; and located in gene regulatory elements. Our findings suggest that protein-coding variants in selected GWAS loci did not contribute to GWAS signals. By performing targeted sequencing in GWAS loci, we identified less-common and rare non-coding SNVs associated with BMD independently from GWAS common SNPs, suggesting both common and less-common variants may associate with disease risks and phenotypes in the same loci.

Prevalence of Pathogenic Mutations in Cancer Predisposition Genes among Pancreatic Cancer Patients.

The prevalence of germline pathogenic mutations in a comprehensive panel of cancer predisposition genes is not well-defined for patients with pancreatic ductal adenocarcinoma (PDAC). To estimate the frequency of mutations in a panel of 22 cancer predisposition genes, 96 patients unselected for a family history of cancer who were recruited to the Mayo Clinic Pancreatic Cancer patient registry over a 12-month period were screened by next-generation sequencing. Fourteen pathogenic mutations in 13 patients (13.5%) were identified in eight genes: four in ATM, two in BRCA2, CHEK2, and MSH6, and one in BARD1, BRCA1, FANCM, and NBN. These included nine mutations (9.4%) in established pancreatic cancer genes. Three mutations were found in patients with a first-degree relative with PDAC, and 10 mutations were found in patients with first- or second-degree relatives with breast, pancreas, colorectal, ovarian, or endometrial cancers. These results suggest that a substantial proportion of patients with PDAC carry germline mutations in predisposition genes associated with other cancers and that a better understanding of pancreatic cancer risk will depend on evaluation of families with broad constellations of tumors. These findings highlight the need for recommendations governing germline gene-panel testing of patients with pancreatic cancer.

A 2-step penalized regression method for family-based next-generation sequencing association studies.

Large-scale genetic studies are often composed of related participants, and utilizing familial relationships can be cumbersome and computationally challenging. We present an approach to efficiently handle sequencing data from complex pedigrees that incorporates information from rare variants as well as common variants. Our method employs a 2-step procedure that sequentially regresses out correlation from familial relatedness and then uses the resulting phenotypic residuals in a penalized regression framework to test for associations with variants within genetic units. The operating characteristics of this approach are detailed using simulation data based on a large, multigenerational cohort.

On family-based genome-wide association studies with large pedigrees: observations and recommendations.

Family based association studies are employed less often than case-control designs in the search for disease-predisposing genes. The optimal statistical genetic approach for complex pedigrees is unclear when evaluating both common and rare variants. We examined the empirical power and type I error rates of 2 common approaches, the measured genotype approach and family-based association testing, through simulations from a set of multigenerational pedigrees. Overall, these results suggest that much larger sample sizes will be required for family-based studies and that power was better using MGA compared to FBAT. Taking into account computational time and potential bias, a 2-step strategy is recommended with FBAT followed by MGA.

ABCC9 gene polymorphism is associated with hippocampal sclerosis of aging pathology.

Hippocampal sclerosis of aging (HS-Aging) is a high-morbidity brain disease in the elderly but risk factors are largely unknown. We report the first genome-wide association study (GWAS) with HS-Aging pathology as an endophenotype. In collaboration with the Alzheimer's Disease Genetics Consortium, data were analyzed from large autopsy cohorts: (#1) National Alzheimer's Coordinating Center (NACC); (#2) Rush University Religious Orders Study and Memory and Aging Project; (#3) Group Health Research Institute Adult Changes in Thought study; (#4) University of California at Irvine 90+ Study; and (#5) University of Kentucky Alzheimer's Disease Center. Altogether, 363 HS-Aging cases and 2,303 controls, all pathologically confirmed, provided statistical power to test for risk alleles with large effect size. A two-tier study design included GWAS from cohorts #1-3 (Stage I) to identify promising SNP candidates, followed by focused evaluation of particular SNPs in cohorts #4-5 (Stage II). Polymorphism in the ATP-binding cassette, sub-family C member 9 (ABCC9) gene, also known as sulfonylurea receptor 2, was associated with HS-Aging pathology. In the meta-analyzed Stage I GWAS, ABCC9 polymorphisms yielded the lowest p values, and factoring in the Stage II results, the meta-analyzed risk SNP (rs704178:G) attained genome-wide statistical significance (p = 1.4 × 10(-9)), with odds ratio (OR) of 2.13 (recessive mode of inheritance). For SNPs previously linked to hippocampal sclerosis, meta-analyses of Stage I results show OR = 1.16 for rs5848 (GRN) and OR = 1.22 rs1990622 (TMEM106B), with the risk alleles as previously described. Sulfonylureas, a widely prescribed drug class used to treat diabetes, also modify human ABCC9 protein function. A subsample of patients from the NACC database (n = 624) were identified who were older than age 85 at death with known drug history. Controlling for important confounders such as diabetes itself, exposure to a sulfonylurea drug was associated with risk for HS-Aging pathology (p = 0.03). Thus, we describe a novel and targetable dementia risk factor.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Large common deletions associate with mortality at old age.

Copy-number variants (CNVs) are a source of genetic variation that increasingly are associated with human disease. However, the role of CNVs in human lifespan is to date unknown. To identify CNVs that influence mortality at old age, we analyzed genome-wide CNV data in 5178 participants of Rotterdam Study (RS1) and positive findings were evaluated in 1714 participants of the second cohort of the Rotterdam Study (RS2) and in 4550 participants of Framingham Heart Study (FHS). First, we assessed the total burden of rare (frequency <1%) and common (frequency >1%) CNVs for association with mortality during follow-up. These analyses were repeated by stratifying CNVs by type and size. Secondly, we assessed individual common CNV regions (CNVR) for association with mortality. We observed that the burden of common but not of rare CNVs influences mortality. A higher burden of large (≥ 500 kb) common deletions associated with 4% higher mortality [hazard ratio (HR) per CNV 1.04, 95% confidence interval (CI) 1.02-1.07, P = 5.82 × 10(-5)] in the 11 442 participants of RS1, RS2 and FHS. In the analysis of 312 individual common CNVRs, we identified two regions (11p15.5; 14q21.3) that associated with higher mortality in these cohorts. The 11p15.5 region (combined HR 1.59, 95% CI 1.31-1.93, P = 2.87 × 10(-6)) encompasses 41 genes, of which some have previously been related to longevity, whereas the 14q21.3 region (combined HR 1.57, 95% CI 1.19-2.07, P = 1.53 × 10(-3)) does not encompass any genes. In conclusion, the burden of large common deletions, as well as common CNVs in 11p15.5 and 14q21.3 region, associate with higher mortality.

Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets.

We present the first analysis of the human proteome with regard to interactions between proteins. We also compare the human interactome with the available interaction datasets from yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans) and fly (Drosophila melanogaster). Of >70,000 binary interactions, only 42 were common to human, worm and fly, and only 16 were common to all four datasets. An additional 36 interactions were common to fly and worm but were not observed in humans, although a coimmunoprecipitation assay showed that 9 of the interactions do occur in humans. A re-examination of the connectivity of essential genes in yeast and humans indicated that the available data do not support the presumption that the number of interaction partners can accurately predict whether a gene is essential. Finally, we found that proteins encoded by genes mutated in inherited genetic disorders are likely to interact with proteins known to cause similar disorders, suggesting the existence of disease subnetworks. The human interaction map constructed from our analysis should facilitate an integrative systems biology approach to elucidating the cellular networks that contribute to health and disease states.