Candidate gene resequencing to identify rare, pedigree-specific variants influencing healthy aging phenotypes in the long life family study

Background: The Long Life Family Study (LLFS) is an international study to identify the genetic components of various healthy aging phenotypes. We hypothesized that pedigree-specific rare variants at longevity-associated genes could have a similar functional impact on healthy phenotypes. Methods: We performed custom hybridization capture sequencing to identify the functional variants in 464 candidate genes for longevity or the major diseases of aging in 615 pedigrees (4,953 individuals) from the LLFS, using a multiplexed, custom hybridization capture. Variants were analyzed individually or as a group across an entire gene for association to aging phenotypes using family based tests. Results: We found significant associations to three genes and nine single variants. Most notably, we found a novel variant significantly associated with exceptional survival in the 3' UTR OBFC1 in 13 individuals from six pedigrees. OBFC1 (chromosome 10) is involved in telomere maintenance, and falls within a linkage peak recently reported from an analysis of telomere length in LLFS families. Two different algorithms for single gene associations identified three genes with an enrichment of variation that was significantly associated with three phenotypes (GSK3B with the Healthy Aging Index, NOTCH1 with diastolic blood pressure and TP53 with serum HDL). Conclusions: Sequencing analysis of family-based associations for age-related phenotypes can identify rare or novel variants

Surveillance Recommendations for Children with Overgrowth Syndromes and Predisposition to Wilms Tumors and Hepatoblastoma

A number of genetic syndromes have been linked to increased risk for Wilms tumor (WT), hepatoblastoma (HB), and other embryonal tumors. Here, we outline these rare syndromes with at least a 1% risk to develop these tumors and recommend uniform tumor screening recommendations for North America. Specifically, for syndromes with increased risk for WT, we recommend renal ultrasounds every 3 months from birth (or the time of diagnosis) through the seventh birthday. For HB, we recommend screening with full abdominal ultrasound and alpha-fetoprotein serum measurements every 3 months from birth (or the time of diagnosis) through the fourth birthday. We recommend that when possible, these patients be evaluated and monitored by cancer predisposition specialists. At this time, these recommendations are not based on the differential risk between different genetic or epigenetic causes for each syndrome, which some European centers have implemented. This differentiated approach largely represents distinct practice environments between the United States and Europe, and these guidelines are designed to be a broad framework within which physicians and families can work together to implement specific screening. Further study is expected to lead to modifications of these recommendations.This study was supported by NCI K08 CA1939915, Alex's Lemonade Stand Foundation for Childhood Cancer, and St. Baldrick's Foundation (to J.M. Kalish); European Research Council Advanced Researcher Award (to E.R. Maher); and NCI 5P30CA054174-21 (to G.E. Tomlinson)

Deep Sequencing of the Nicastrin Gene in Pooled DNA, the Identification of Genetic Variants That Affect Risk of Alzheimer's Disease

Nicastrin is an obligatory component of the γ-secretase; the enzyme complex that leads to the production of Aβ fragments critically central to the pathogenesis of Alzheimer's disease (AD). Analyses of the effects of common variation in this gene on risk for late onset AD have been inconclusive. We investigated the effect of rare variation in the coding regions of the Nicastrin gene in a cohort of AD patients and matched controls using an innovative pooling approach and next generation sequencing. Five SNPs were identified and validated by individual genotyping from 311 cases and 360 controls. Association analysis identified a non-synonymous rare SNP (N417Y) with a statistically higher frequency in cases compared to controls in the Greek population (OR 3.994, CI 1.105–14.439, p = 0.035). This finding warrants further investigation in a larger cohort and adds weight to the hypothesis that rare variation explains some of genetic heritability still to be identified in Alzheimer's disease

PoPoolation: A Toolbox for Population Genetic Analysis of Next Generation Sequencing Data from Pooled Individuals

Recent statistical analyses suggest that sequencing of pooled samples provides a cost effective approach to determine genome-wide population genetic parameters. Here we introduce PoPoolation, a toolbox specifically designed for the population genetic analysis of sequence data from pooled individuals. PoPoolation calculates estimates of θWatterson, θπ, and Tajima's D that account for the bias introduced by pooling and sequencing errors, as well as divergence between species. Results of genome-wide analyses can be graphically displayed in a sliding window plot. PoPoolation is written in Perl and R and it builds on commonly used data formats. Its source code can be downloaded from http://code.google.com/p/popoolation/. Furthermore, we evaluate the influence of mapping algorithms, sequencing errors, and read coverage on the accuracy of population genetic parameter estimates from pooled data

Profiles of Multidrug Resistance Protein-1 in the Peripheral Blood Mononuclear Cells of Patients with Refractory Epilepsy

BACKGROUND: About one third of patients with epilepsy become refractory to therapy despite receiving adequate medical treatment, possibly from multidrug resistance. P-glycoprotein, encoded by multidrug resistance protein-1 (MDR1) gene, at the blood brain barrier is considered as a major factor mediating drug efflux and contributing to resistance. Given that peripheral blood mononuclear cells (PBMNCs) express MDR1, we investigated a MDR1 status of PBMNCs in various subsets of epilepsy patients and demonstrated their association with clinical characteristics. METHODOLOGY/PRINCIPAL FINDINGS: Clinical and MDR1 data were collected from 140 patients with epilepsy, 30 healthy volunteers, and 20 control patients taking anti-epileptic drugs. PBMNCs were isolated, and basal MDR1 levels and MDR1 conformational change levels were measured by flow cytometry. MDR1 profiles were analyzed according to various clinical parameters, including seizure frequency and number of medications used in epilepsy patients. Epilepsy patients had a higher basal MDR1 level than non-epilepsy groups (p<0.01). Among epilepsy patients, there is a tendency for higher seizure frequency group to have higher basal MDR1 level (p = 0.059). The MDR1 conformational change level was significantly higher in the high-medication-use group than the low-use group (p = 0.028). Basal MDR1 (OR = 1.16 [95% CI: 1.060-1.268]) and conformational change level (OR = 1.11 [95% CI: 1.02-1.20]) were independent predictors for seizure frequency and number of medications, respectively. CONCLUSIONS/SIGNIFICANCE: The MDR1 profile of PBMNCs is associated with seizure frequency and medication conditions in patients with epilepsy

Two Evolutionary Histories in the Genome of Rice: the Roles of Domestication Genes

Genealogical patterns in different genomic regions may be different due to the joint influence of gene flow and selection. The existence of two subspecies of cultivated rice provides a unique opportunity for analyzing these effects during domestication. We chose 66 accessions from the three rice taxa (about 22 each from Oryza sativa indica, O. sativa japonica, and O. rufipogon) for whole-genome sequencing. In the search for the signature of selection, we focus on low diversity regions (LDRs) shared by both cultivars. We found that the genealogical histories of these overlapping LDRs are distinct from the genomic background. While indica and japonica genomes generally appear to be of independent origin, many overlapping LDRs may have originated only once, as a result of selection and subsequent introgression. Interestingly, many such LDRs contain only one candidate gene of rice domestication, and several known domestication genes have indeed been “rediscovered” by this approach. In summary, we identified 13 additional candidate genes of domestication

Rare Variants in APP, PSEN1 and PSEN2 Increase Risk for AD in Late-Onset Alzheimer's Disease Families

Pathogenic mutations in APP, PSEN1, PSEN2, MAPT and GRN have previously been linked to familial early onset forms of dementia. Mutation screening in these genes has been performed in either very small series or in single families with late onset AD (LOAD). Similarly, studies in single families have reported mutations in MAPT and GRN associated with clinical AD but no systematic screen of a large dataset has been performed to determine how frequently this occurs. We report sequence data for 439 probands from late-onset AD families with a history of four or more affected individuals. Sixty sequenced individuals (13.7%) carried a novel or pathogenic mutation. Eight pathogenic variants, (one each in APP and MAPT, two in PSEN1 and four in GRN) three of which are novel, were found in 14 samples. Thirteen additional variants, present in 23 families, did not segregate with disease, but the frequency of these variants is higher in AD cases than controls, indicating that these variants may also modify risk for disease. The frequency of rare variants in these genes in this series is significantly higher than in the 1,000 genome project (p = 5.09×10−5; OR = 2.21; 95%CI = 1.49–3.28) or an unselected population of 12,481 samples (p = 6.82×10−5; OR = 2.19; 95%CI = 1.347–3.26). Rare coding variants in APP, PSEN1 and PSEN2, increase risk for or cause late onset AD. The presence of variants in these genes in LOAD and early-onset AD demonstrates that factors other than the mutation can impact the age at onset and penetrance of at least some variants associated with AD. MAPT and GRN mutations can be found in clinical series of AD most likely due to misdiagnosis. This study clearly demonstrates that rare variants in these genes could explain an important proportion of genetic heritability of AD, which is not detected by GWAS

Sequencing of high-complexity DNA pools for identification of nucleotide and structural variants in regions associated with complex traits

We have used targeted genomic sequencing of high-complexity DNA pools based on long-range PCR and deep DNA sequencing by the SOLiD technology. The method was used for sequencing of 286 kb from four chromosomal regions with quantitative trait loci (QTL) influencing blood plasma lipid and uric acid levels in DNA pools of 500 individuals from each of five European populations. The method shows very good precision in estimating allele frequencies as compared with individual genotyping of SNPs (r(2) = 0.95, P &lt; 10(-16)). Validation shows that the method is able to identify novel SNPs and estimate their frequency in high-complexity DNA pools. In our five populations, 17% of all SNPs and 61% of structural variants are not available in the public databases. A large fraction of the novel variants show a limited geographic distribution, with 62% of the novel SNPs and 59% of novel structural variants being detected in only one of the populations. The large number of population-specific novel SNPs underscores the need for comprehensive sequencing of local populations in order to identify the causal variants of human traits

A single active catalytic site is sufficient to promote transport in P-glycoprotein

P-glycoprotein (Pgp) is an ABC transporter responsible for the ATP-dependent efflux of chemotherapeutic compounds from multidrug resistant cancer cells. Better understanding of the molecular mechanism of Pgp-mediated transport could promote rational drug design to circumvent multidrug resistance. By measuring drug binding affinity and reactivity to a conformation-sensitive antibody we show here that nucleotide binding drives Pgp from a high to a low substrate-affinity state and this switch coincides with the flip from the inward- to the outward-facing conformation. Furthermore, the outward-facing conformation survives ATP hydrolysis: the post-hydrolytic complex is stabilized by vanadate, and the slow recovery from this state requires two functional catalytic sites. The catalytically inactive double Walker A mutant is stabilized in a high substrate affinity inward-open conformation, but mutants with one intact catalytic center preserve their ability to hydrolyze ATP and to promote drug transport, suggesting that the two catalytic sites are randomly recruited for ATP hydrolysis

Rare event detection using error-corrected DNA and RNA sequencing

Conventional next-generation sequencing techniques (NGS) have allowed for immense genomic characterization for over a decade. Specifically, NGS has been used to analyze the spectrum of clonal mutations in malignancy. Though far more efficient than traditional Sanger methods, NGS struggles with identifying rare clonal and subclonal mutations due to its high error rate of ~0.5–2.0%. Thus, standard NGS has a limit of detection for mutations that are &gt;0.02 variant allele fraction (VAF). While the clinical significance for mutations this rare in patients without known disease remains unclear, patients treated for leukemia have significantly improved outcomes when residual disease is &lt;0.0001 by flow cytometry. In order to mitigate this artefactual background of NGS, numerous methods have been developed. Here we describe a method for Error-corrected DNA and RNA Sequencing (ECS), which involves tagging individual molecules with both a 16 bp random index for error-correction and an 8 bp patient-specific index for multiplexing. Our method can detect and track clonal mutations at variant allele fractions (VAFs) two orders of magnitude lower than the detection limit of NGS and as rare as 0.0001 VAF