Family-based exome sequencing and case-control analysis implicate CEP41 as an ASD gene.

Autism Spectrum Disorder (ASD) is a complex neurodevelopmental disorder with a strong genetic component. Although next-generation sequencing (NGS) technologies have been successfully applied to gene identification in de novo ASD, the genetic architecture of familial ASD remains largely unexplored. Our approach, which leverages the high specificity and sensitivity of NGS technology, has focused on rare variants in familial autism. We used NGS exome sequencing in 26 families with distantly related affected individuals to identify genes with private gene disrupting and missense variants of interest (VOI). We found that the genes carrying VOIs were enriched for biological processes related to cell projection organization and neuron development, which is consistent with the neurodevelopmental hypothesis of ASD. For a subset of genes carrying VOIs, we then used targeted NGS sequencing and gene-based variant burden case-control analysis to test for association with ASD. Missense variants in one gene, CEP41, associated significantly with ASD (p = 6.185e-05). Homozygous gene-disrupting variants in CEP41 were initially found to be responsible for recessive Joubert syndrome. Using a zebrafish model, we evaluated the mechanism by which the CEP41 variants might contribute to ASD. We found that CEP41 missense variants affect development of the axonal tract, cranial neural crest migration and social behavior phenotype. Our work demonstrates the involvement of CEP41 heterozygous missense variants in ASD and that biological processes involved in cell projection organization and neuron development are enriched in ASD families we have studied.

De novo variants in congenital diaphragmatic hernia identify MYRF as a new syndrome and reveal genetic overlaps with other developmental disorders.

Congenital diaphragmatic hernia (CDH) is a severe birth defect that is often accompanied by other congenital anomalies. Previous exome sequencing studies for CDH have supported a role of de novo damaging variants but did not identify any recurrently mutated genes. To investigate further the genetics of CDH, we analyzed de novo coding variants in 362 proband-parent trios including 271 new trios reported in this study. We identified four unrelated individuals with damaging de novo variants in MYRF (P = 5.3x10(-8)), including one likely gene-disrupting (LGD) and three deleterious missense (D-mis) variants. Eight additional individuals with de novo LGD or missense variants were identified from our other genetic studies or from the literature. Common phenotypes of MYRF de novo variant carriers include CDH, congenital heart disease and genitourinary abnormalities, suggesting that it represents a novel syndrome. MYRF is a membrane associated transcriptional factor highly expressed in developing diaphragm and is depleted of LGD variants in the general population. All de novo missense variants aggregated in two functional protein domains. Analyzing the transcriptome of patient-derived diaphragm fibroblast cells suggest that disease associated variants abolish the transcription factor activity. Furthermore, we showed that the remaining genes with damaging variants in CDH significantly overlap with genes implicated in other developmental disorders. Gene expression patterns and patient phenotypes support pleiotropic effects of damaging variants in these genes on CDH and other developmental disorders. Finally, functional enrichment analysis implicates the disruption of regulation of gene expression, kinase activities, intra-cellular signaling, and cytoskeleton organization as pathogenic mechanisms in CDH.

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.

Exome Sequencing of Phenotypic Extremes Identifies CAV2 and TMC6 as Interacting Modifiers of Chronic Pseudomonas aeruginosa Infection in Cystic Fibrosis.

Discovery of rare or low frequency variants in exome or genome data that are associated with complex traits often will require use of very large sample sizes to achieve adequate statistical power. For a fixed sample size, sequencing of individuals sampled from the tails of a phenotype distribution (i.e., extreme phenotypes design) maximizes power and this approach was recently validated empirically with the discovery of variants in DCTN4 that influence the natural history of P. aeruginosa airway infection in persons with cystic fibrosis (CF; MIM219700). The increasing availability of large exome/genome sequence datasets that serve as proxies for population-based controls affords the opportunity to test an alternative, potentially more powerful and generalizable strategy, in which the frequency of rare variants in a single extreme phenotypic group is compared to a control group (i.e., extreme phenotype vs. control population design). As proof-of-principle, we applied this approach to search for variants associated with risk for age-of-onset of chronic P. aeruginosa airway infection among individuals with CF and identified variants in CAV2 and TMC6 that were significantly associated with group status. These results were validated using a large, prospective, longitudinal CF cohort and confirmed a significant association of a variant in CAV2 with increased age-of-onset of P. aeruginosa airway infection (hazard ratio = 0.48, 95% CI=[0.32, 0.88]) and variants in TMC6 with diminished age-of-onset of P. aeruginosa airway infection (HR = 5.4, 95% CI=[2.2, 13.5]) A strong interaction between CAV2 and TMC6 variants was observed (HR=12.1, 95% CI=[3.8, 39]) for children with the deleterious TMC6 variant and without the CAV2 protective variant. Neither gene showed a significant association using an extreme phenotypes design, and conditions for which the power of an extreme phenotype vs. control population design was greater than that for the extreme phenotypes design were explored.

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 tumor's clonal composition; however, deconvolving NGS data to infer a tumor's 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.

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.

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.

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.

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.