Debbie A. Nickerson

Actionable exomic incidental findings in 6503 participants: challenges of variant classification

Recommendations for laboratories to report incidental findings from genomic tests have stimulated interest in such results. In order to investigate the criteria and processes for assigning the pathogenicity of specific variants and to estimate the frequency of such incidental findings in patients of European and African ancestry, we classified potentially actionable pathogenic single-nucleotide variants (SNVs) in all 4300 European- and 2203 African-ancestry participants sequenced by the NHLBI Exome Sequencing Project (ESP). We considered 112 gene-disease pairs selected by an expert panel as associated with medically actionable genetic disorders that may be undiagnosed in adults. The resulting classifications were compared to classifications from other clinical and research genetic testing laboratories, as well as with in silico pathogenicity scores. Among European-ancestry participants, 30 of 4300 (0.7%) had a pathogenic SNV and six (0.1%) had a disruptive variant that was expected to be pathogenic, whereas 52 (1.2%) had likely pathogenic SNVs. For African-ancestry participants, six of 2203 (0.3%) had a pathogenic SNV and six (0.3%) had an expected pathogenic disruptive variant, whereas 13 (0.6%) had likely pathogenic SNVs. Genomic Evolutionary Rate Profiling mammalian conservation score and the Combined Annotation Dependent Depletion summary score of conservation, substitution, regulation, and other evidence were compared across pathogenicity assignments and appear to have utility in variant classification. This work provides a refined estimate of the burden of adult onset, medically actionable incidental findings expected from exome sequencing, highlights challenges in variant classification, and demonstrates the need for a better curated variant interpretation knowledge base.

Recurrent de novo mutations implicate novel genes underlying simplex autism risk

Autism spectrum disorder (ASD) has a strong but complex genetic component. Here we report on the resequencing of 64 candidate neurodevelopmental disorder risk genes in 5,979 individuals: 3,486 probands and 2,493 unaffected siblings. We find a strong burden of de novo point mutations for these genes and specifically implicate nine genes. These include CHD2 and SYNGAP1, genes previously reported in related disorders, and novel genes TRIP12 and PAX5. We also show that mutation carriers generally have lower IQs and enrichment for seizures. These data begin to distinguish genetically distinct subtypes of autism important for etiological classification and future therapeutics.

Detection of structural variants and indels within exome data

We report an algorithm to detect structural variation and indels from 1 base pair (bp) to 1 Mbp within exome sequence data sets. Splitread uses one end–anchored placements to cluster the mappings of subsequences of unanchored ends to identify the size, content and location of variants with high specificity and sensitivity. The algorithm discovers indels, structural variants, de novo events and copy number–polymorphic processed pseudogenes missed by other methods.

Inferring Clonal Composition from Multiple Sections of a Breast Cancer

Cancers arise from successive rounds of mutation and selection, generating clonal populations that vary in size, mutational content and drug responsiveness. Ascertaining the clonal composition of a tumor is therefore important both for prognosis and therapy. Mutation counts and frequencies resulting from next-generation sequencing (NGS) potentially reflect a tumors clonal composition; however, deconvolving NGS data to infer a tumors clonal structure presents a major challenge. We propose a generative model for NGS data derived from multiple subsections of a single tumor, and we describe an expectation-maximization procedure for estimating the clonal genotypes and relative frequencies using this model. We demonstrate, via simulation, the validity of the approach, and then use our algorithm to assess the clonal composition of a primary breast cancer and associated metastatic lymph node. After dividing the tumor into subsections, we perform exome sequencing for each subsection to assess mutational content, followed by deep sequencing to precisely count normal and variant alleles within each subsection. By quantifying the frequencies of 17 somatic variants, we demonstrate that our algorithm predicts clonal relationships that are both phylogenetically and spatially plausible. Applying this method to larger numbers of tumors should cast light on the clonal evolution of cancers in space and time.

Genome-Wide Association Study Identifies Genetic Loci Associated with Iron Deficiency

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged ≥25 y and women ≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF) ≤ 12 µg/L (cases) and iron replete controls (SF>100 µg/L in men, SF>50 µg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10−7 for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P = 7.0×10−9, corrected P = 0.012) was replicated within the VA samples (observed P = 0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.

Association between patterns of nucleotide variation across the three fibrinogen genes and plasma fibrinogen levels: the Coronary Artery Risk Development in Young Adults (CARDIA) study

Summary.  Background: Previous genotype–phenotype association studies of fibrinogen have been limited by incomplete knowledge of genomic sequence variation within and between major ethnic groups in FGB, FGA, and FGG. Methods: We characterized the linkage disequilibrium patterns and haplotype structure across the human fibrinogen gene locus in European‐ and African‐American populations. We analyzed the association between common polymorphisms in the fibrinogen genes and circulating levels of both ‘functional’ fibrinogen (measured by the Clauss clotting rate method) and total fibrinogen (measured by immunonephelometry) in a large, multi‐center, bi‐racial cohort of young US adults. Results: A common haplotype tagged by the A minor allele of the well‐studied FGB−455 G/A promoter polymorphism (FGB 1437) was confirmed to be strongly associated with increased plasma fibrinogen levels. Two non‐coding variants specific to African‐American chromosomes, FGA 3845 A and FGG 5729 G, were each associated with lower plasma fibrinogen levels. In European‐Americans, a common haplotype tagged by FGA Thr312Ala and several other variant alleles across the fibrinogen gene locus was strongly associated with decreased fibrinogen levels as measured by functional assay, but not by immunoassay. Overall, common polymorphisms within the three fibrinogen genes explain < 2% of the variability in plasma fibrinogen concentration. Conclusions: In young adults, fibrinogen multi‐locus genotypes are associated with plasma fibrinogen levels. The specific single nucleotide polymorphism and haplotype patterns for these associations differ according to population and also according to phenotypic assay. It is likely that a substantial proportion of the heritable component of plasma fibrinogen concentration is due to genetic variation outside the three fibrinogen genes.

De Novo and Rare Variants at Multiple Loci Support the Oligogenic Origins of Atrioventricular Septal Heart Defects.

Congenital heart disease (CHD) has a complex genetic etiology, and recent studies suggest that high penetrance de novo mutations may account for only a small fraction of disease. In a multi-institutional cohort surveyed by exome sequencing, combining analysis of 987 individuals (discovery cohort of 59 affected trios and 59 control trios, and a replication cohort of 100 affected singletons and 533 unaffected singletons) we observe variation at novel and known loci related to a specific cardiac malformation the atrioventricular septal defect (AVSD). In a primary analysis, by combining developmental coexpression networks with inheritance modeling, we identify a de novo mutation in the DNA binding domain of NR1D2 (p.R175W). We show that p.R175W changes the transcriptional activity of Nr1d2 using an in vitro transactivation model in HUVEC cells. Finally, we demonstrate previously unrecognized cardiovascular malformations in the Nr1d2tm1-Dgen knockout mouse. In secondary analyses we map genetic variation to protein-interaction networks suggesting a role for two collagen genes in AVSD, which we corroborate by burden testing in a second replication cohort of 100 AVSDs and 533 controls (p = 8.37e-08). Finally, we apply a rare-disease inheritance model to identify variation in genes previously associated with CHD (ZFPM2, NSD1, NOTCH1, VCAN, and MYH6), cardiac malformations in mouse models (ADAM17, CHRD, IFT140, PTPRJ, RYR1 and ATE1), and hypomorphic alleles of genes causing syndromic CHD (EHMT1, SRCAP, BBS2, NOTCH2, and KMT2D) in 14 of 59 trios, greatly exceeding variation in control trios without CHD (p = 9.60e-06). In total, 32% of trios carried at least one putatively disease-associated variant across 19 loci,suggesting that inherited and de novo variation across a heterogeneous group of loci may contribute to disease risk.

FOXE3 mutations predispose to thoracic aortic aneurysms and dissections.

The ascending thoracic aorta is designed to withstand biomechanical forces from pulsatile blood. Thoracic aortic aneurysms and acute aortic dissections (TAADs) occur as a result of genetically triggered defects in aortic structure and a dysfunctional response to these forces. Here, we describe mutations in the forkhead transcription factor FOXE3 that predispose mutation-bearing individuals to TAAD. We performed exome sequencing of a large family with multiple members with TAADs and identified a rare variant in FOXE3 with an altered amino acid in the DNA-binding domain (p.Asp153His) that segregated with disease in this family. Additional pathogenic FOXE3 variants were identified in unrelated TAAD families. In mice, Foxe3 deficiency reduced smooth muscle cell (SMC) density and impaired SMC differentiation in the ascending aorta. Foxe3 expression was induced in aortic SMCs after transverse aortic constriction, and Foxe3 deficiency increased SMC apoptosis and ascending aortic rupture with increased aortic pressure. These phenotypes were rescued by inhibiting p53 activity, either by administration of a p53 inhibitor (pifithrin-α), or by crossing Foxe3-/- mice with p53-/- mice. Our data demonstrate that FOXE3 mutations lead to a reduced number of aortic SMCs during development and increased SMC apoptosis in the ascending aorta in response to increased biomechanical forces, thus defining an additional molecular pathway that leads to familial thoracic aortic disease.

Associations between common fibrinogen gene polymorphisms and cardiovascular disease in older adults - The Cardiovascular Health Study

Elevated plasma fibrinogen is a risk factor for cardiovascular disease (CVD), but associations between fibrinogen single nucleotide polymorphisms (SNPs) and disease risk are inconsistent. We investigated whether common (> or = 5% minor allele frequency) variation in the fibrinogen genes (FGA, FGB, FGG) is associated with fibrinogen concentration, carotid artery intima-medial thickness (IMT) and risk of incident myocardial infarction (MI), ischemic stroke and CVD mortality in European- (EA) and African-descent (AA) adults (> or = 65 years) from the Cardiovascular Health Study. TagSNPs were genotyped in 3,969 EA and 719 AA free of MI or stroke at baseline. Race-specific models included multiple testing correction and adjustment for sex, age and site. Among EA, minor alleles of FGA3807, FGB1437 and FGG902 were associated with higher fibrinogen levels; whereas FGA251, FGA2224, FGA6534 and FGG10034 were associated with lower levels, p<0.004 for each. Strongest associations were seen for FGB1437; each additional copy of the minor allele was associated with 13 mg/dl (95%CI: 9-16) higher fibrinogen level. Similar trends in AA were not significant. Fibrinogen haplotypes were not significantly associated with internal or common carotid IMT. No associations with MI or CVD mortality were seen in EA, though FGB1038 and FGG902 were significantly associated with increased and decreased risk of stroke in men, respectively, as were related haplotypes. FGB1038 was also associated with CVD mortality in AA, HR = 1.9 (95%CI: 1.3-2.7). In conclusion, while fibrinogen genetic variation was strongly associated with fibrinogen levels, there was less evidence of association with the more complex outcomes of IMT and CVD events.

Association Between Celiac Disease and Iron Deficiency in Caucasians, but Not Non-Caucasians

BACKGROUND & AIMS Celiac disease is an increasingly recognized disorder in Caucasian populations of European origin. Little is known about its prevalence in non-Caucasians. Although it is thought to be a cause of iron-deficiency anemia, little is known about the extent to which celiac disease contributes to iron deficiency in Caucasians, and especially non-Caucasians. We analyzed samples collected from participants in the Hemochromatosis and Iron Overload Screening study to identify individuals with iron deficiency and to assess the frequency of celiac disease. METHODS We analyzed serum samples from white men (≥25 y) and women (≥50 y) who participated in the Hemochromatosis and Iron Overload Screening study; cases were defined as individuals with iron deficiency (serum ferritin level, ≤12 μg/L) and controls were those without (serum ferritin level, >100 μg/L in men and >50 μg/L in women). All samples also were analyzed for human recombinant tissue transglutaminase immunoglobulin A; positive results were confirmed by an assay for endomysial antibodies. Patients with positive results from both celiac disease tests were presumed to have untreated celiac disease, and those with a positive result from only 1 test were excluded from analysis. We analyzed HLA genotypes and frequencies of celiac disease between Caucasians and non-Caucasians with iron deficiency. RESULTS Celiac disease occurred in 14 of 567 cases (2.5%) and in only 1 of 1136 controls (0.1%; Fisher exact test, P = 1.92 × 10(-6)). Celiac disease was more common in Caucasian cases (14 of 363; 4%) than non-Caucasian cases (0 of 204; P = .003). Only 1 Caucasian control and no non-Caucasian controls had celiac disease. The odds of celiac disease in individuals with iron deficiency was 28-fold (95% confidence interval, 3.7-212.8) that of controls; 13 of 14 cases with celiac disease carried the DQ2.5 variant of the HLA genotype. CONCLUSIONS Celiac disease is associated with iron deficiency in Caucasians. Celiac disease is rare among non-Caucasians-even among individuals with features of celiac disease, such as iron deficiency. Celiac disease also is rare among individuals without iron deficiency. Men and postmenopausal women with iron deficiency should be tested for celiac disease.